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“数值 分 析 " 是 解决 工科 数学 问题 和 工程 实际 问题 的 重要 理论 基 
础 和 实用 工具 ,也 是 工科 各 专业 博士 研究 生 人 学 考试 的 内 容 。 为 了 
帮助 广大 学 生 更 好 地 学 习 和 人 掌握 数值 分 析 课 程 的 理论 精 艇 和 解 题 方 
法 ,我 们 根据 清华 大 学 出 版 社 出 版 的 由 李 类 扬 \ 王 能 超 、 易 大 义 编写 
的 《数值 分 析 } 教 材 的 章节 顺序 ,以 其 内 容 为 基础 编写 了 这 本 辅导 教 
材 。 本 书 共 分 9 章 , 每 章 均 有 重点 .难点 全 析 和 习题 全 解 两 个 板块 。 
书 来 附 了 三 套 自 测试 题 及 答案 ,帮助 掌 生 自我 检测 学 习 效 果 。 

(1) 重 点 .难点 全 析 :精练 地 列 山 了 各 章 的 主要 知识 点 , 理 清 了 各 
知识 点 之 间 的 脉络 联系 , 列 出 了 主要 定理 及 其 相关 推论 ,重要 公式 
等 ,帮助 读者 融会 贯通 ,系统 理解 各 章 的 体系 结构 , 问 定 扎实 的 理论 
基础 。 

《2) 习 题 全 解 ,对 数值 分 析 各 章 的 本题 作 了 详细 解答 ,在 解答 过 
程 中 ,对 重点 和 难点 习题 进行 了 分 析 和 讲解 ,归纳 了 解 题 技巧 。 

我 们 本 着 服务 于 读者 学 习 、 考 试 的 角度 ,通过 对 习题 的 解答 , 揭 
示 了 数值 分 析 的 解 题 方法 ,规律 和 技巧 ,以 帮助 读者 更 好 地 理解 该 课 
程 的 基本 概念 和 理论 ,开拓 解 题 思路 ,所 高 数学 素质 。 本 书 由 杨刚 、 
起 燕 , 王 宇 翔 编写 ,在 编写 过 程 中 ,得 到 了 李强 . 黄 蔬 等 同志 的 支持 和 
帮助 ,在 此 表示 衷心 的 感谢 ! 对 《数值 分 析 }》 的 作者 李 庆 扬 , 王 能 超 ，、 
易 大 义 载 授 表 示 囊 心 的 感谢 ! 

由 于 编者 水 平 有 限 ,本 书 难 免 存在 缺漏 和 不 妥 之 处 ,请 各 位 专家 
和 广大 读者 批评 指正 。 


编 者 
2007 年 3 月 
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第 1 章 绪论 


1.1 重点 ,难点 全 析 


1,1.1 误差 来 源 及 分 类 


数值 分 析 主 要 研究 两 类 误差 . 

1, 省 入 刘 差 

出 于 计算 机 字 长 有 限 , 原 始 数 据 在 计算 机 上 的 表示 以 及 进行 算术 运算 十， 
一 ， 多 , 二) 时 可 能 产生 的 误差 , 称 为 会 和 误差. 

2. 截断 误 益 (方法 误差 ) 

为 了 在 有 限 步 内 得 到 运算 结果 , 用 有 限 的 过 程 职 代 无 帘 的 过 程 时 产生 的 
误差 , 称 之 为 截断 误差 或 方法 误差 . 

用 数值 方法 求解 数学 模型 所 得 到 的 近似 解 包 含 含 人 误差 和 截断 误差 两 
部 分 . 


.1.2 误差 度量 了 
1. 绝对 误差 [误差 ) 
(1) 定义; 近似值 x" 与 准确 值 z 之 差 称 为 碑 对 误差 ,简称 误 莽 .用 公式 表示 


为 etz' ) -2 -= 在 不 引起 混 请 时 可 简 记 e(zr ) 为 e- . 

(2) 绝对 误差 限 (误差 限 ), 绝对 误差 绝对 值 的 一 个 上 界 ,用 符 导 表示 为 
eftz") 或 6. 

2， 相对 广 差 


(1) 定义 :绝对 误差 。" 与 准确 值 x( 非 零 ) 的 比值 乞 一 三 -二 三 称 为 近似值 


中 ”实际 应 用 中 ,误差 指 绝对 误差 限 , 相 对 误差 指 相对 误差 限 . 
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x” 的 相对 误差 ,用 符号 表示 为 e, (x") 或 e? . 由 于 真 值 x 总 是 未 知 的 , 故 
本 
了 


时 - 
当 zx" 较 好 地 近似 z 时 ,两 种 方法 仅 相差 高 阶 无 穷 小 . 
《2) 相对 误差 限 ; 相 对 误差 可 正 可 负 , 它 的 奖 对 值 的 较 小 上 界 叫 居 相 对 误 
差 限 , 记 必 e.Cx') 或 er , 即 用 


计算 ， 
3- 有 效 副 字 
(1) 定义 : 设 准 确 值 x 的 近似 值 x"( 非 零 ) 可 表示 为 
工 ” 一 士 10" X (Cate X10 + 二 a X10 "ny C1.1) 
式 人 ,1 中 心心 一 1 是 0 一 9 中 的 一 个 数字 ,al 关 0,m 为 整数 , 旦 使 


上 二 ”一定 ij 和 也 X ] Dr ”+1 


成 立 的 整数 ” 的 最 大 值 为 册 则 称 z 近似 x 具有 hn 位 有 效 数 字 }aj ,as，*…， as 
分 别 是 z” 近似 工 的 第 1 位 .第 2 位 ……、. 第 # 位 有 效 数 字 - 

(2)? 有 效 数 : 称 近似 值 未 位 也 是 有 效 数 字 的 近似 数 为 有 效 数 . 此 时 ,绝对 误 
差 限 不 超过 末 位 单位 的 一 半 . 进而 可 以 知道 ,对 准确 值 进行 四 舍 五 人 或 对 近似 
数 的 准确 位 进行 四 会 五 人 所 得 近似 值 均 为 有 将 数 . 有 效 数 的 有 效 数 字 位 数 等 
于 去 起 第 一 位 非 零 数 字 到 末 位 数字 的 位 数 ， 

4 三 种 度量 间 的 联系 

《1) 绝对 误 差 与 相对 误差 的 关系 参 钢 相对 误差 的 定义 . 

23 物 对 误差 与 有 效 数 字 的 英雄 参 闪 有效 数 字 的 定 文 - 

(3) 相对 误差 与 有 效 数字 的 关系 参阅 教材 ， 

定理 1.1 设 非 零 近似 数 z” 有 形 如 式 (1. 1) 的 表示 ,车 7' 具有 ?位 有 效 
数字 , 则 相对 误差 限 为 


2a 
定理 1.3 设 非 零 近似 数 ”有形 如 式 (],1) 的 表示 ,车 +" 的 相对 误差 
限 为 


1 _ 
0 
和 24a 十 1 
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则 x" 至 少 有 * 位 有 效 数字 
114.3 初 值 误差 传播 


1， 一 元 画 数 误差 传播 
ELF OF = fr} flr} a f(r" elr') 
etflr ad fir’) letr’) 
EC 


其 中 条 件数 Cf ay 于 {| 


2. 并 元 济 教 误 盖 情 播 . 
Er 


ma Cee) 
和 2 Ex tr 开 
，， | 二 J et: 
ez ， 人 
ER 
下 好] 
.af 
def 一 
Ctf, x) Ty 


frr 
3, 二 元 莫 本 运 舞 的 误差 传播 
ErY 二 T= elrr 十 氏 工 
EEC 十 ECxz ) 


Ea 本 | 工 ? | ECXY } 十 | 证 | el.r2 ) 
(Ee ee te es), 


到 [a 到 玉 0 


grr 直 a max{e CAT Et) Te 0 


Er Cxi Xp ya etar 了 二 人， tr 天 日 


1.1.4 算法 的 数值 稳定 性 


一 个 算 落 如 果 输 人 数据 有 误 莽 , 而 在 计算 过 程 中 售 人 误差 不 增长 , 则 称 此 
算法 是 数值 稳定 的 ,否则 称 此 算法 为 数值 不 稳定 的 ， 
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1.1.5 算法 设计 中 避免 误差 危害 的 蒋 干 原则 


(1) 要 避免 除数 种 对 值 远 远 小 于 被 除数 饱 对 人 革 的 除法 ， 
(2) 要 避免 两 相近 数 相 减 ， 

(3) 要 防止 大 数 “ 吃 挤 ” 小 数 . 

“4) 注意 简化 计算 步骤 ,减少 运算 次 数 ， 


1.2 习题 全 解 


1. 谈 z Ds 的 相对 误差 为 训 , 求 hz 的 误差 . 
解 。 nr: -lnr ~ 一 } :相对 误差 


{etlnr) | er jd 
进而 有 etlntz’ sc 全 
2， 证 了 的 由 对 误差 为 2 各 ,来 如 的 相对 谍 盖 ， 
解 函数 值 的 条 件数 为 
Cr ri) 一 
根据 函数 信 误 差 传 播 笃 式 
ET I a Opry Er = nx 2 = DD. 02n 
3. 于 列 各 救 部 是 经 过 四 会 五 入 得 到 的 近似 数 , 即 误差 限 不 超过 最 后 一 伍 的 
夺 丰 单位, 试 指出 它们 有 由 位 有 北 数 宇 . 
Ti = .1021, zi =0.03], 4: = 385.6, zx/ = 56,430, x7 一 7X1.0 
解 因为 近似 数 的 误差 限 不 超过 最 后 一 位 的 半 个 单位 ,所 以 这 些 数 均 为 有 
效 数 . 而 有 效 数 的 有 效 数字 位 数 为 第 一 位 非 零 数 到 最 后 一 位 数字 的 位 数 , 所 以 
六 有 5 位 有 效 数 字 ， zi? 有 2 位 有 效 数字 ， xf 有 4 位 有 效 数字 ， zy 有 
5 位 有 效 数 字 ， xz? 有 2 位 有 效 数字 
4， 利 用 公式 (2.3)5 给 差 竺 播 公式 ) 计算 下 列 各 近似值 的 误差 限 : 


xz" ACzx* yl | 
Cn 


Cr Yn 


【FTwT 十 < 十 2 rr [和 | 二 | 2 
重 


其 中 X22 ， 嫩 3 均 为 和 3 是 所 给 的 载 . 
解 eta 一 间 5X104eTre = O05X10 er 05X]D 0， 
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stm Y= O05X10 ,etry = 0.5X10! 
Ciyetr?r +r? ry = etxr tr ) 十 ET 一 
ox10 ox107 二 O05 10: 一 
1.05 x 1073 


(i) etr? wi ri as | wi x et y+ a a [etawit Hirr a | etri)} = 
0,031 xX 385. 站 共有 510 1 十 
1.102 1 x 385.6 x 0.5X10 3 十 


1.102 I xX 0,.031 XxX 0,5X 10 = 0,21479 


i 人 | ws | etxs 十 | x | etr?) 
《而 》 (六 )= Cr 
56. 430 XxX 0.5 X10 :++0.031 x 0, 5x10° 
56, 430° 


0.886 54 x 10 


5. 计算 球体 积 要 司 和 对 讲 差 跟 汶 1%, 问 度量 半径 民 时 克 许 的 相 烤 误差 限 
是 多 少 ? 


解 ”计算 球体 积 V 一 xR' 近似 值 的 条 件数 为 


cov' 一 | ) | RArR'Y 
VR) 针 x(R*) 


应 用 误差 传播 公式 有 
EV a CV Je, (R') — 90.(R*) 
令 3e,{R' ) 二 1, 可 得 度量 半径 尺 时 允许 的 相对 误差 限 为 
s.(R') 一 寺 X1% a 0.003 3 
6. 说 下 = 28, 接 道 推 公式 


Y, 一 了 ,一 可 VI, n= 1 2 


计算 到 Yo- 著 取 v783 二 27, 人 兰 有 多 大 误 益 ? 
解 YY 一 Yi 一 一 V78 一 (Y,， 


1 
100 io0 ™ 783) 商 Y "8 一 
1 
YY， 一 高 v783 = (YY, 高 V783 )— D5 vv 了 3 = 
Y，, 一 5 /83 = Y, v783 


— 05 
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鼓 有 Yiw 一 加 一 TBT Vo YY VI Yh, = YY — 27.982 


| et |=| Yh — Yo |=| v8 — 27., 982 | 于 X10 = ein) 


7， 来 广 竹 天 一 58z 十 1 二 0 的 丙 修 根 , 使 它 至 少 具有 4 位 有 效 教 宁 ( w783 
a 27, 982), 
解 ”根据 求 根 会 式 可 知 方 程 的 根 为 
zz 一 28 士 783 


XL 2B wTB3 a 28 27, 982 一 35. 982 


0, O17 B63 


1 1 
Ta 二 28 783 一 se 一 
28+ wTE3 55.982 


故 根 据 误差 传播 公式 可 知 两 根 近 似 秆 均 满 足 题 设 要 求 . 
. Ft! 1 ， 
8. 当 和 充分 大 时 , 怎 料 来 | 、 二 du 
和 解 | 让 二 和 一 arctan( N+ 1) — arctant y) se a 一 让 二 


tana— tang _ 
1 十 tanatanf 


arctant ranta 一 站 一 arctan 
arctan — ETN 一 arctan i 
IT 十 CVD i 二 + N+ NN 
9， 正 方形 的 边 长 太 约 为 100 cm, 应 塌 样 测 醒 方 能 使 其 面积 误差 不 超过 
] ecm’? 
和 解 。 设 正 方形 边 长 为 4a, 面积 为 时 ,于 是 有 
M=a, MM’ =(a):, ethM*}y=1, a = 100 


应 用 误差 传 播 公式 有 
El(M' sa sa } = 200eta’ } 
令 200eta*) 志 1 
得 se(a' EO.X1o: 


故 正 方形 边 长 误差 限 应 不 超过 0. 005 cm, 才能 使 其 面积 误差 不 超过 1 em’. 
10. 设 S 二 方 B1: ,假定 g 是 准确 的 ,而 时 :的 测量 有 土 0. is 的 误差 ,证 明 当 
1 增加 时 ,5S 的 绝对 误差 增加 ,而 相对 误差 却 减少 . 
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证 明 ”由 题 意 知 eli) 一 0.1， 5* 一 主 g() 
绝对 误差 限 
ES | mr i elf) = 0. at" 
相对 误差 限 
Sy KS) Vlg 0.2 
” SS” 1 2 全 
Ett } 


由 于 :增加 时 ,其 近似 值 :* 也 增加 , 故 可 知 题 自 结 论 正确 ， 
11. 序列 1y1 满足 遵 推 关系 
yr = Oy 1 C1, n= 1,2,. 
车 6 二 党 1.414 三 位 有 效 数 生 }) ,计算 到 yo 时 误差 有 台大 ?这 个 计算 过 程 稳 
定 吗 ? 
解 多 一 10 一 1 一 1000o 2 一 17 一 ] = 10y :一 [1 十 101] = 
10010y s = DD 一 [1 二 1m] = 10y,5 一 [1 十 10!: 寺 10:] 一 


| 
= 10"3o 一 > 10 = 10°y;— (0 一 1) = 


10"(E 一 言 ) 上 + 二 


一 n 1 | 。 1 1 1 
yw = 10" (VE— 5 )1 二 10" (4 和 1 一生)+ 语 
绝对 误差 
[ety [= | 一 m0 1.41) 过 


lu 1 2 1 E 
lo x 2 X10 = 2 X10 


计算 到 yu 时 误差 为 万 X 10", 帮 在 绝对 误差 的 意义 下 该 方法 不 稳定 . 


12. 计算 站 二 (2 一 10, 取 Y2 sx 1.4, 利 用 下 列 等 式 计算 , 哪 一 个 得 到 的 寻 
果 最 好 ? 


(3 — 2 v2):, 399 一 90 


1 上 
WE+ 1 (3 + 2D) 
解 ”99 一 ?2 出 现 相近 数 相 减 ,(3 一 2y2) 出 现 较 相近 的 数 相 减 , 故 二 者 不 


8 数 午 分 析 全 析 糊 人 解 


可 能 得 到 最 好 的 结果 . -一 1 一 与 一 | 均 不 出 现 相 近 数 相 碱 ,但 前 者 
{f+ (3 十 2 


的 乘法 运算 次 数 过 ,而 除法 运算 次 数 相同 ， 每 次 秉 除法 运算 必然 引信 新 的 会 人 


误差 ， 故 用 式 一 工 计算 将 得 到 最 好 的 运算 结果 、 
(3 十 2 


13, Fiz) 二 ln(z 一 v2 一 1), 求 1130) 的 值 , 车 开平 市 用 6 位 函 教 表 , 问 来 
对 数 时 诅 善 有 部 大 ?车 改 用 嚼 一 等 价 公 式 
nr— vr DD) = in(rt ve 1) 
计算 , 求 对 数 时 误差 有 和 多大? 
和 解 /30) = ln(30 一 v399), 设 w= 二 v899,y 一 了 (30); 则 jw' 一 29.9833， 
1 


所 以 st ) =. x 10 ,页 


1 2 1 
EC TT Tt 了 TD Fete” 2 a 0, O03 


必用 等 价 公 式 ln ~ 一 1) = 一 ln(r+ vi 一 1) 


则 Fg0) =— Int30 十 本 9) 
,|. .1 
故 Ely*) = | 0 etu”) = 


1 和 -局 
5 5355 本 EC ) 冬日. 834 XX 1 
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第 2 章 插值 法 


1.1 重点 、 难 点 全 析 


1.1.1 播 值 问题 


1]. 定义 

设 函 数 y< 并 z) 在 区 间 [a, 的 上 有 定 祥 ,用 已 知 在 点 a 芝 二 过 癌 二 二 
za 所 由 上 的 值 ov 若 存 在 -简单 画 数 Pix) ,使 得 

Pt 一 1 《2, 1) 

成 立 , 刚 称 Pic) 是 fir) 关于 节点 rz oz 的 -一 个 揪 值 函数 . 

点 .more 称 为 插值 节点 ;[a, 站 称 为 插值 区 间 : 了 (x) 为 被 插值 函数 ; 
式 {2. 41} 即 为 播 值 杀 件 ， 

用 PLr) 的 利 作 略 六 zr) 的 近似 鼻 , 当 工 在 节点 形成 的 区 间 上 有 时 , 称 该 方法 为 
肉 插 法 ; 当 工 不 在 节点 形成 的 区 间 上 但 在 插值 区 间 上 时 , 则 称 该 方法 为 外 插 法 . 

当 播 值 孙 数 Piz) 为 包 项 式 时 , 称 Ptr) 是 六 rz) 的 一 个 播 值 多 项 式 , 相 点 的 
插值 法 称 为 多 项 式 插值 . 


插值 余 项 RC -至 jr) 一 Pix) ,插值 余 项 又 称 为 规 断 误差 


2. 插值 多 项 式 的 在 在 唯一 性 定理 

定理 2.1 满足 插值 条 件 (2. 1) 的 不 起 过 nn 次 的 播 值 儿 项 式 P tx} 是 存 在 唯 
一 的 . 

推论 2.2 若 jz) 是 不 超过 m 次 的 多 项 式 , 则 它 的 关于 十 ] 个 互 异 节点 
zt 的 不 超过 = 次 的 插值 客 项 式 Ptx) 与 被 插值 沙 数 f(z) 恒 等 , 即 
Pr fir). 

3. 误差 估计 

(1) 导数 型 误差 情 计 定理 ， 
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定理 2.3 设 六 x) EC "Da yt， 产 mCt 在 ta 如 内 存在 . 撒 值 余 项 


产 ” (E&Y 


GFT te) 


R(t) = Fx Ptr) = 
def 


式 中 Ptzx) 是 六 rr) 关于 二 殿 节 点 x wil4… ,x 的 插值 名 项 式 ,wws (Tr} 
{ta 且 依赖 于 工 和 插值 节点 . 

《2) 尽 莽 {也 称 洲 商 ) 型 误差 司 计 定理 : 

定理 2.4 没 六 x) 在 a, 办 ,上 有 定义 ,P(x) 同 定理 ?2.3 中 的 说 明 , 则 有 插 
值 余 硕 

RT) = Fry Tas TL rs Th ep Cr) 

(3) 等 距 节 点 的 插值 误差 估计 ,; 当 节 点 zw wr ,tyits 等 距 分 布 时 , 即 有 
如 十 下 0 二 01， 1, 作 灾 换 了 = zi， 则 导数 再 误差 估 
计 简 化 为 

Rizr rth = fir th}— Ptr, rth} oe 


tet itt— pH) 
tn 十 3}1 


如 果 作 变换 一 z, 十 1 , 则 误差 估计 简化 为 
R(T, FA) 一 fr t+) Pt7s 十 1 上) 一 


tt 二 1- tt ny 
Cn 17)! 


er] Fe EE Cro, ra) 


htl Fi {Ey » £ EE (xn 了 


1.1.2 搬 值 多 项 式 的 构造 方法 


1. 佐 乌 插值 条 件 妇 立 线 性 方程 组 
设 相 超过 座 的 播 值 狠 项 式 Ptz) 一 上 十 可 工 十 azr 十 十 cu" 其 系数 
是 如 下 线性 方程 组 的 和解 : 


2. Lagrange 插值 方法 


def 
PA = fr) (TT fro tt fr La (2.2) 
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式 中 :Cr 和: 是 关于 节点 和 ra 的 Lagrange 揪 值 基 随 数 , 它 仅 与 节点 
有 有关. 与 被 插 导 函数 无 闫 .具体 定义 为 


【一 
(一 一 一 


[重地 | 


Wt CD) (k= 0 1 2 sn 


(一 Tt) 
插值 公式 52. 2) 称 为 Lagrange 插 秆 多 项 式 ， 
3. ewton 档 慎 方法 
Ptry = Fro flrs rr zt fry Rs er rT 


Fir Ts ra ET Tr 1) A N(x)y {2.3) 
式 中 三 m， zy 2] 为 束 阶 莽 商 ,插值 公式 (2. 3) 称 为 Newtan 插值 公式 . 
4， 等 距 节 可 播 值 公式 
基于 变换 7 = 十 thth 洁 0) 的 插值 公式 (Newton 癌 前 插值 会 式 ) 


lap, 十 … 十 


I l(tt— + 由、 


nl 
式 中 Atfot& 二 1，2,"…, nn) 为 上 阶 向 前 差分 . 
基于 变换 二 一 加 十 上 人 全 0) 的 插值 公式 (Newton 向 后 插值 公式 >》 


Ne 上 Ti = f+ + 十， 


Nr th) = f+ tA 二 


入 (2.4) 


tf 二 + 1 (ttn 1) 
nl 


式 中 本 ft 庄 二 1， 2 …, 1) 为 站 阶 向 后 差分 . 
1.1.3 世 函 数 、 差 商 和 差分 的 性 质 


]. Lagrange 插值 基 泡 数 
11， 下 一 了 
《Del 一 | 《一 1，… 9) 
志和 关 了 
{2 二 ra] 
C30 rit (tr) zitiad i dr 0 
3. 均 差 
(1) 莽 商 的 函数 值 瑚 示 ，; 


LT (2. 51 
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二 zj) 
fr re 之 [| 
A fr) 
| {i 


寺中 wtf 二 Cr 一 一 ri). 
C2) 对 称 性 ; 差 商 与 节点 的 排列 次 序 无 关 . 
(3) 线性 组 合 性 质 : 
te f+ esa rir ry re = fr Fy i 二 
ca BL Tr es ee 
式 中 心 及 为 常数 ， 
{4) 上 与 导数 的 关系 : 


四 
Ex,, TIT], Eo 一 全 . EE [a.#] 


3 差分 ( 竺 距 节点 } 


CD af CE Df My 人 (js 
J 


pn 


VA -FY DrDT 人 jos 
4 了 


mean 
(2) fn = TAP = ED (js]n 
i 
《3) fF[x,, J -= a (m= 2) 
3 pr Tl mn pe rs 
1 1 
mm 


td) A fs = (BY, rh EE Th 


Flas TE ea] 一 


1.1.4 Hermite 质 值 


1， 插 值 条 忻 
Hi = Fr), His}y= 7 tn) (j=0,1, "m,n) (2. 6) 
2. 在 碍 叭 一 性 
满足 插值 条 伴 式 (2.6) 的 不 焉 过 2x 十 1 次 的 插值 名 项 式 是 存在 唯一 的 . 
3. 大 襄 娄 表示 形式 
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Ha tx) = DF (er) tf Cr pr)] {2.7) 


I 


式 中 mn 12) pln Cr rt) 


Ca 


4. 误 羡 司 计 
定理 2.5 设 六 站 ECmtay af 在 (av 加 内 存在 , 则 有 误 善 
估计 
2 i 
Rs; tT) = fir}— Hetr) 一 i Ta (CX), EE (Cahy 
式 中 上 依 坊 于 及 插值 节点 . 
1.1.$ 分 段 播 值 


1. 分 藤 线 性 播 值 
(1) 插值 形式 :函数 f(x) 关于 节点 ae 一 居 过 下 习 < 一 五 的 耸 毁 线 
性 手 值 函数 


lr) = Et fr) pr) 
TE el 还 kr1 I 
X11 
(2) 谋 差 估计 : 
BM, 
max | 所 一 tr) | 起 站 
a 
def . def def 
式 中 MM 2 max i Fry |, 用 一 下 max hs 
时 Il 


{3) 收 化 性 : 当 六 xz) EE Cla, bl 时 有 limh tr) 二 xz 在 [ay 5 上 一 致 
成 立 . 

2. 分 段 三 次 Hermite 桥 情 . 

(1) 插值 形式 :函数 TD 关于 节点 TT 一 记 的 分 段 三 
次 Hermite 插值 函数 


hn = (1+2 二 一 下 (生殖 -) jen) 1 


TH A hl 


(SR ) df 


Tr ty 


{1 + 2 Ee ) (= ) fe + 


Xr Tt Trttl Tk 
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{ A ) (me yf rn ) 


Tht -Tt 
rr 尼 [x Ti | {A 一 0 1,2, ry nO—1) 


def 
— max i Dr 1 则 


《2) 误差 箔 计 : 当 .Ar) E Ciia, 四 且 有 如 


=! — A 
| RLir) | 一】 Fiz) Oo tr) | 381 


式 中 Amax 入 ， < Ti 一 工 ， 
当 产 z) € Cad] ME max : PCr)y | 时 , 则 
ax | fr nM | 
(3) 收 笋 性 : 当 Fx) € Ci'[a, 及 时 , 则 limli tx) 二 f(z) 在 [ay #5] 上 一 至 
成 立 - 
3, 三 次 样 杂 插值 
《1) 定义 ;着 酒 数 Stx} 所 Ci[a' 台 , 且 在 每 个 小 区 间 [z， x 1](ji 一 9， 
l1，…… ,一 1) 上 是 三 次 包 项 式 , 刚 称 Sr? 是 美 于 节点 & 二 式 过 人 
二 上 的 三 次 样 条 函数 . 若 进 一 步 有 
ST) = Friy C=, 1, m,n) 
则 称 SCz) 为 二 次 样 条 插值 西数 . 
(2) 存在 唯 … 性 :满足 下 列 三 组 附加 条 件 之 一 的 三 次 样 条 插值 函数 是 唯一 
存在 的 . 
第 一 组 (转角 每 件 }; Sr) = 二 (Cr, SCrm) = f(x)s 
第 二 给 ( 弯 第 茉 件 ); Stxe) 一 (Cr) Sz) 二 产 (rn 
第 三 组 (周期 条 件 ): Sa 十 0 一 SC 一 们 ,i 二 0, 1,2. 
(3) 插值 形式 : 


3 四 2 
sz 一周 二 2 二 ‘I i +- (ez 2 ) 2 之 | 
了 人 


(Fa i) 主 一 也 + 在 Lz; 了 it ], 1 = DO, 1, :NO—1 
6 Ah 


式 中 局 二 .| 一 工 ,: 
对 于 第 一 组 边界 条 件 ,参数 1A4 ) 六 。 是 下 列 方程 组 的 解 : 


2 1 [Le xa, x 
I 2 Al fr», ls zs] 


= 6 A | 


站 + Tl Ei 


Hr I Tn zx, 


上 hh 


局 
~ def ftir) fn, x 加 Fr rr 一 产 (Cr 


en -Tn 1 
区 | Ah, 


式 中 i 一 人 一 1 2 


; -3 Fr Tr Fx) ffir, 1 Xa 


六 1 Tn ml 


FL 3 
这 样 方 程 组 右 端 的 差 商 部 分 可 用 带 重 节 点 的 差 商 表 计 算 , 则 有 


To 了 Co 
FF Crs) 
Tn pe) > > f[zr,, os x1] 


a fn fy tn. x a] 
zs fem) > 1 


实 呈 1 Fixe } 


= /Lz, 1 + | 


Fw 六 ~ pr,) fr Xa Xs] 
Tr fix) 
对 于 第 :组 边界 条 件 ,1 和 是 下 烈 方 程 组 的 解 
2 和 Al MT 6FLx 4 I} * 了] | BE 
2 2 由 BM; 6 FL x LE 1] 
tH 2 3 BE 站 站 rs za X14] 
We 2 5 MI; BF[xo 9, rn.2s Ye 1 
Hn 1 2 m -] FL r,s Tl sj — Ast MM 


At 一 Fr M, 一 Fr} 
对 于 第 三 组 边界 条 忻 , {M17-o 是 下 列 方程 组 的 解 (M, 一 Mi 》 
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Lxro, Xi kid 
2 jn fo Fir Tr 1 
此 2 a {Lz x Ts 
2 Pp 天 ss ， Xa EE 
" 。 . 二 : 
尼 2 A FL-cr ? ol Tn ] 
1 呈 1 
2 fx: 9 = fro bz 
Fr " hh 十 Re: 
fy he 
式 中 富 守 局 ' 所 一 二 1 二 


(4) 误差 估计 及 收 训 b 性 : 
定理 2.6 设 站 xz) 亡 世 [a; 有 和 ,Stz) 是 f(x) 美 于 节点 2 二 < 二 < 
之 I 二 为 的 三 次 样 条 插值 函数 ,附加 第 一 组 或 第 三 组 边界 条 件 . 则 有 误差 估计 
max | 一 SI) | OC max | Fx) | Ri {m= 0,1,2) 
式 中 ty 一 人 《1 一 于， ta 一 证 h 二 max A, hs = ZH 一 kr 


战 上 述 定理 自然 得 到 , 当 睛 -9 时 ,Sry，S zy 种 汪 tr) 均 分 别 一 致 收 笋 
到 Fry, F(z) 及 (zr). 


1.2 习题 全 解 


1,， 当 了 工 一 11 一 1]12 时 XY) 一 0. 一 3, 和 ,未 FE) 的 二 次 质 值 多项式， 
解 Lr = fr tn Fri + Flr)l(r) = 
{xT (tr 2) 4 《一 工作 工 十 1 


83 
计 开 十 卫 f 一 站 


2. 给 出 fr) 二 nx 数 慎 下 如 下 ，; 


.4 


0.8 0.6 ,7 心 .只 


—0.9185291 —0.693 147 —0.510 826 C—O0.356 75 一心,223 144 


用 绕 性 后 值 及 二 次 插值 计算 In0. 54 的 近 亿 全 
解 ”依据 插值 误差 估计 式 选 距离 0. 54 较 近 的 点 为 插值 节点 ,并 建 空 差 商 
表 如 下 ， 
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zo =0.5 一 0.693147 
> 1. 823 210 
二 0.6 一 0.510826 > —0.204 115 


2. 027 325 
二 0.4 一 0.916 291 -一 


写 出 Newton 插值 多 项 式 
NCT) =— 0,693 147 十 1.823 2]0tr oe 0,.3) 


N(x) = NCD 0,204 115(r— 0.5{r—0.6) 
计算 近似 值 
NCO. 54) 一 一 0.693 147 十 1.823 210¢0. 54~—0.5) =—0.620 218 6 
Na tO. 54) = Ni (0.54) ~ 0.204 115¢0, 54—0, 53)(0.54— 0.6) =—0.616 839 


3， 给 出 cosT'0 之 工 芝 90" 的 函数 表 , 步 长 上 一 1 一 (高 ) ,车 丁 数 表 具 有 


5 性 有 效 数 字 .研究 用 线性 持 值 求 rosr 近似 情 时 的 蕊 误差 界 ， 
解 ”用 王 数 值 及 近世 值 所 建立 的 线性 揪 慎 和 多项式 为 
当 .c ELzr, zm- 时. 


一 . 了 一 -tl ， ， 人 

Er fn) Et re 
Lr 一 fr fi) -一 
| -1 一 工 ， 

Tn a ee 
式 中 心 E 1 一 1 人 0 1 2 5 400). 
得 到 误差 估计 
| cosz 一 上 7 tr l= | cosre Ltr TEL 一 工人 2 | < 


| eosz 一 上 fr || (xr) — Lr Cr) | 
证 以 看 出 ,如 果 和 忽略 算 术 运 算 过 各 中 引入 的 会 人 误差 ,插值 误差 包含 截断 误差 ， 
cos 一 上 LitxX) | 及 初始 数据 误差 的 传播 两 部 分 ， 
截断 误差 


' Comr Oo Lit) := 让 (~ Sinf (TC TT | 


| tr rr 一 ri | 


max | (r(x 一 cin) | 
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1[i a 2 a 
TLE 0 “1X 
TE Lx, xers) 
合作 误 闫 传播 
[Et |= eA DD | ef ra | 
Ti A To 


! “Er | -+ . I 了 了 本 一 上 工 ; 
maxt| et fF ty | | etf' (rn) 人 本 二 元 ) 
: , ,def 
maxil etf* Ce | etf’ (rT) 11 SS Mn 
利用 有 效 数字 定 义 有 


ea 亲本 xl0 


式 中 mm 为 函数 值 5x) 的 量 级 , 故 有 有 


10% | 六 (Cr 1< 0ri 
进而 Mi 过 magk| X10, 


me 1 rai 
五 X10mn| 一 诗 X10 + 


cosr— Ti ir | .06 x10 二 区 YX LO mr 


TE tIiv Td T= O01, 2. 5 390 
在 区 间 | 0， 也 | 上 的 总 误差 界 为 

| cosr — L? (x) (A 1 0B XxX ID + 二 XX In = 0.50]1 05 Xx 10™ 
4. 设 .i 为 至 看 节点 人 一 Donl 求证: 
C1) Slt x) se RR 0, 1 


7 


Cy Dr — a {ry ed (k= 1 2， 


证 明 i107 令 区 二 式 ,车 插值 节点 为 xz 人 二 01.ryR), 则 ftz) 的 
次 播 值 党 项 式 为 


了 zy 一 Dl tr) 


了 一 中 


i 
Ra) Fr 一 Try = 区 全 ep 


插值 余 项 为 
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及 因 沪 让 记 所 岂 . 
FT =0 R(t)=0 


所 以 人 ty = HR = 0 an) 


a | 于 [3 
Cy) Dr Cr) = bp ED (se zy | 
i 全 一 ; 


一 


5 > jz 《一 ei co |= 


J=0 = 
和 1 一 ; - 


5 > (0 Hn |= 


Fe 了 


[em Dat | 


点 
之 
生 下 . 
> Ek tr) 
em 


5. 设 六 7 ECLa bl 和 有 fa) 王 了 (= 0 求证: 


max | FLEX) |} 所 计 作 一 6? 严 max | F(z) 
可 二 Fr 


证 上 明 令 x=&n 和 z=b， 以 此 为 播 值 节点 . 则 插值 多 项 式 为 


_ 一生 ， 工 一 对 ， 
Litr) fm sp + fib) Fe 习 


局 


应 用 插值 余 项 公式 有 
| ft) ~ Itr) ,= Orr) Eg EE (ay hb) 
六 max | (es) | max {realtr—p) lO 
襄 (6 一 a) Inax | Fr) i 
6. 在 一 所 信 4 上 上 给 出 六 Xx} 一 * 的 等 距 节 志 庆 数 表 , 若 用 分 段 二 次 祷 
值 求 眠 的 近 位 值 , 要 使 截断 误差 不 超过 10， 问 使 用 函数 表 的 砂 长 声 应 取 务 少 ? 
解 。” 若 插值 节点 为 x;_; ， zx; 和 zt 则 分 惨 一 次 播 值 移 项 式 的 捅 值 余 项 为 


Riz) 一 尊严 (Bt Ti TE Ts EE Ci Tri) 


式 中 i 一 点 ， A 二 十 中， 
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i R(t) | 过 二 max | Cec) ror | 
TS 
1 ,2 1，,， el 1 
如 上 k= 
6 3 of 
ep 
他 一 一 js 二 10” 得 上 所 0.006 58. 
gy 和 
插值 点 个 数 
4 一 (一 和 二 1 ?1 <] 了 ef 3 
1 二 3 1 208 1217 TN 
星 奇 数 , 故 实际 可 采用 的 函数 值 表 步 长 
= 0.006 579 


?7. 车 二 2 来 妨 及 站 
解 ”根据 向 前 差分 算 子 和 中 心 莽 分 算 子 的 定义 进行 求解 


1 4 


4 于 
Aly, = (E— Diy, = D(C) DEy, = 2 j= Diy = 
t 


i 


+ 4 
(> (Dm j>。 = (21Hy, y= 2 
ll 了 


六 FE 一 站 一 
FiAy, = Ey) = ys = 2 
8. 如 时 fx) 是 向 次 多 区 式 , 记 坟 F(z) 一 Fr 二 ji 一 Fr) 征明 了 1( 的 
下 阶 差分 各 FCD 直 让 坟 89) 是 如一 下 忧 和 多项式, 并 且 AMmiFftxr) 二 0(f 为 正 
蒜 教 )， 
证 明 对 mm 次 儿 项 式 六 应 用 Taylor 公式 有 


AP) = Crz 十 站 一 rn = 二 人 二 十 村 fr) 


即 afFfz) 为 m 一 1 次 的 多 项 式 、 

A 二 BACTD), 对 黄 一 > 次 事项 式 Aftx) 应 用 上 述 推理 过 程 知 
ACAFLT)) 二 Bz) 是 mm 一 2 次 的 委 项 式 . 

悚 此 过 程 递 推 , 知 Arz)t0 芝 二 过 由 1) 为 由 一 上 次 密 项 式 ， 

所 以 aa" Frz) 为 常数 ,前 aa Pr) 寺 0 为 正 整 数 )， 

9. 证 明 站 六 Bt 让 让 Agi 十 gf 

证 明 Atfigt) = fge 一 fga = fi ge ~ fgul 二 figsn 一 大 有 一 
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Beltfet— ft folgi — Ba) = EA fag 
10. 证 明 fog = fg jog 一 DA 
证 明 ”由 于 
SY pag, 十 Sg af = 一 号 was 十 gi 太太 一 Sana = 


上 = t= 


Dy Cf ge 一 frgey 一 FBt — fg: 


En 
ml +] 

所 以 DP fag frge— fp — gnafs 
Lint =n 


mi| 


11. 证 明 DY ay; = Ay Ay 
T= 


wr | wl 
证 明 Dy ay 一 Dalay,) 一 DY Ap — Ay) = (Ay — Ay:) 
i=n i=? ~ 

Ay A) = Ay 


12, 著 六 让 二 十 a 十 十 gw 1.x" 十 ar 有 于 个 和 不同 实 究 工 re rn， 
or 证 阴 : 


a 此 一 于 一 】 
证 明 ”因为 fix) 有 个 不 同 实 根 ix;%1; 且 
FT) = aT Te 


{I A = ok 


= x En _1 
2 te Cr 上 六 
记 gtx) 二 x ,并 利用 荆 商 的 请 数值 表达 形 还 有 
A 1 CX;) i 
DPT EO BT nel + Xn] 
再 由 差 商 与 异 数 移 关系 知 
中 ， 9, 0 一 2 
zi 一 El 
Fe | =al, k=nl 
! ” (一 1 机 


13, 证 明 n 阶 均 差 有 下 到 性 质 : 
{j)》 著 Fizry 一 crtzryy 册 下 ro is 
{攻关 Fr 一 8 十 呈 C)y 剧 


[rr 
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Fais wr es x = fr Tr es To gin wis mr] 
证 有 明 设 Fiz) 二 gf tz 十 遍 Cr wmi(7T) 一 【一 Te) 一 开工 一 了 
将 差 商 { 均 差 } 用 函 教 值 表示 , 则 有 


Era tr es = p> P(x) -之 af tr 寺庙 (xz) 


TEs1 "Cx ) Ver1 rT Cy 


j= 


D+ 
af Lz Ti + Bbc: Ts rs rnd 

取 庙 一 0, 9 二 c 得 结论 (1 ); 取 a 二 加 二 1 得 结论 (有 

14, fr) 二 十 37 十 11 求 C27 ,21 ,mz ?及 [2'， 28] 

解 ” 人 略为 f(r 二 工 十 十 3rt 十 1; 若 这 二 2'(1 二 0 出 [rn ， 


了 Ta EE 一 二 1 
A221 


HL2", 2 2 前 一 0 
15, 证 明 两 点 三 次 Hermite 揪 情 人 条 项 是 


41 
R; (ry) = er TA rr EE tri, rn) 


并 由 此 采 击 分 段 三 灾 Hermite 插 慎 的 误差 限 ， 
证 明 车 z E Lz， zs;], 且 插值 争 项 式 满足 条 件 
有 = few} Htrer) = fix) 
Hitrn) = fim), Hlrzm) = fF (rn) 
知 播 值 余 项 Rt) 一 六 xr) 一 Ha(lxr) 有 二 重 零 点 rm 和 zi 故 设 
Rix} 一 Rx Ti) ro re )’ 
确定 晴 数 请 (.72: 
当 工 = mm 或 xul 时 ktx) 取 任 何 有 限 值 均 可 
当头 时 EE Cn as) 构造 基于 变量 : 的 国 数 
gtr = FN — Hytty — hE ryt zr) 
显然 有 
RCzt DD, ELT) =0,. grey 二 站 
外 = 0 EH (ren 寺 们 
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在 [aa，xj 和 [rz 上 对 gL.r) 岳 用 Rolle 定理 :存在 六 人 (zx,， XX) 及 
灾 E kz， zn) 使 得 
和) 一 由 (p) = 0 
在 Car h) 《hh 次} 《R17 上 对 g(x) 使 用 Rolle 定理 ,存在 z EE 
frr 让) 让 万 《加 ， 交 )》 和 各 滩 .111 《m+ ar 使 得 
g 《六 = 全 (hs = ET. 1) = 
再 依次 对 中 (和 名 人 使 用 Rolle 定理 , 知 至 少 存在 EE 0z zn) 使 得 
6) 0 
而 0 一 了 0) 一 RCr)41, 将 代入 ,得 到 
Mr) Hf EE€ (ya 


推导 过 程 表明 依 球 于 Th 了 1 贡 世 
综合 以 上 这 程 可 知 
R(xY = ~ ri) Cr re) 
下 面 建立 分 段 二 次 Hermite 插值 的 误差 限 . 记 天 1z 为 关 x) 在 La 到 上 的 
基于 等 谍 节 点 的 分 慌 三 次 Hermite 插 值 隙 数 , z= 二 a 十 台 f 全 二 0, 1, Rn) ,三 一 
刷 一 芷 


nn 四 


在 小 区 间 [x 3 sh 上 有 


FD) | 一 十 | | (rz 一 ar 一 二 


1 max [FD max (ro— x {to re )? 
二 C15 TT 
而 最 值 
一 如 十 电 
max (rx) zr 3 ii maxs: ts— 1)*h! 一 Lp 
| i 16 
进而 得 误差 估计 


re 1 Fei 
| Fir) hr) se sgn Tax Fr) i 


16. 求 一 个 次 数 不 南 于 4 次 的 多 项 式 P(tz}, 使 它 满 是 P(0) 一 P'(0) = 0， 
Pty = P(tl) = 1,P2) 一 1 
解法 一 箱 用 Hermite 插值 可 得 到 次 数 不 商 于 4 的 多项式 


一 一] 二 0 一] m= Dm = 1 
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1 1 
再 ;Cr = Pyotr) Dmp Cr) 
了 一 有 


gr = 1) 
A | i 


mt) = 12 


2 Y: {3— 2r)r: 
Bir) 一 rir 1y 
遍 二 Cx 一 Dr 
所 以 Hitr) = 人 一 2 rr = 
设 已 6z) 二 日 :Cr 十 六 Cr 一 (Cr 一 ,其 中 .六 为 待定 常数 . 令 P(t2) = 
1 得 


于 是 
Pix} 一 二 了 tre 3 


解法 二 ( 带 重 节点 的 Newton 插值 芒 } ”建立 如 下 差 商 表 ， 


ha 二 十 
一 这 
一 过 
[= 上 
| | 
I 一 
中 | 


， 1 


这 样 可 写 出 Newton 搬入 公民 
Plry = d+ ot 0+ 


ltr—0F lr or om1)+ 


OT 


二 ~ 3)? 


17. 设 f(r) = 5 上 取 n 二 10, 按 符 王 节点 求 分 医 线 
性 播 值 汪 数 tT), 计算 各 节点 闻 中 上 筷 处 (zr) 与 Cry 的 值 , 并 居 计 谍 善 . 
租车 一 5iru 一 5, 风 步 长 和 A 二 呈 一 一 了 一 1,7, 一 一 5 十 讨 一 一 5 二 


it0sis10). 在 区 间 2x,. .1] 上 ;分 段 线性 插值 酉 数 为 
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二 一 全 ,1 | EE i 
WO = fo) Tt fn) 
i Ti - 
+ » = 人 DO, 1], "9 
了 十 + 1 十 二] : 
分 段 线性 插值 鲨 数 定义 如 下 : 
RD = = Ef rn! 


1 了 x 1 十 i 


各 节点 问 中 点 处 函数 值 玉 播 值 函 数值 如 下 所 示 : 


fF 庄 DD.3 十 1.3 土 2,8 二 3.3 土 4.5 
Ar) 0. B00 oO DQ. 307 了 0.137 9 nN.075 3 0. DA7 1 


Tn bry 赂 7730 入 dQ. 330 0 0.1500 [9 D, nd 
估计 误差 ,在 区 同 5r.，zrr] 上 
| fi f(x) 1 一 Fr ZrO 
1 max | FE)! max 1 tr ryt iny | 
2 -srs i 
而 
了 二 TT 十 刺 1 
max | ir—a}tr mr.}! es max | sG 一 17 1 二 地 
Pe] rl : 
2 _ 2x 一 Br 一? 
/一 


2 — zx: 
令 扩 (7) 一 2 二 ~ 0 得 六 (>) 的 此 点 0, 士 1, 于 是 


max ¢, Fr = ma FO | | 产 十 11. 天 一 号 | 一 2 
碟 有 辣 论 


| Fen 一 WP 和 让 x2xXL oi, re t,x] 


厂 端 与 i 无关, 故 
| FO 地 rT 35, 5] 
18. 来 fx) 一 在 Le 上 的 分 段 线 性 插值 强 数 J;ir), 并 合计 误差 ， 


和 解 在 区 则 [es 的 上 上, rs 二 二 1) 一 


max pA,, 


ss 1 


86 数 重 分 析 宝 析 格 解 


晃 数 fr) 在 Lz 全 7pP1] 上 的 钱 性 彝 值 喝 疗 为 


下 一 Tl Fri1) 了 一 工 一 王 【Ti 一 工 ) Ei {reo Ti) 
:一 Ti Ti A 上 下， ph; 
分 段 线性 插值 巴 数 
ltr) = I (Cr) 《air 一 工 ) 十 对 《rr 一 Ti 了 EL] 


误差 估计 


| Fzry Bi Cr) | 一 | 二 ff CHT OT | 


让 inax | Fer | max 1 fir){tro rl) | 一 
忆 | 
1 A - 
和 人) 和 eco 
故 
一 rs ce | "eT A 上 
[A | maxr | AT 一 下 ma 二 一 个 
[EE ut 


19. 束 ft) 一 工 在 [ay pl 上 的 分 段 Hermite 插值 .并 信和 计 误 攻 。 
解 在 区 间 Le.p] 上 ,rr = 一 再 由 二 一 max A,. 
在 区 间 ix, rr] 的 Hermite 插值 酉 数 为 


On = ) Hero 全 ) | 


Ti Ti ,it] | 
' 
工 一 二 1] TT 工 . 
Tl (En) + 
At] 一 了， tl 
; z 
$ Ek 四 - 
4 (一 xD 人 ， TE Lr. si] 
Tl 


牛 计 误差 
| fn — Ninery |= [| 


去 max Fr | max (rr){tr— rr)|’ = 
a 了 


rE A 
译 X24X| 香 ] = 证， rE Lr, 7..] 


对 于 人 tx} 有 


~ Ah: 页 
一 | ry 一 
| Fry — Et) (max | Fr 一 上 人) [a Max TH 一 1 
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20. 始 定 载 据 表 加 下 : 


i 25 Bb. 30 39 0D. 45 .3 


yy D, OD 站 D0, 547 了 0, 82 5 D.670 8 0, 728 0 


证 求 三 凑 样 条 插值 SCr) ,并 满足 条 尾 

Ci Sst.25) = 1.000 0 30.53) 一 0.6868: 

(Hy SO. 25) = $00,53) = 0D. 

解 A = 0,.30—0.25= 0.05, = 0.3—0.30= 0.09 
he = D5 0.3 = .05, k= 0.53—0,45= 0.08 


h 用 


J L -= 
4 和 == 
hh 


得 


建立 善 商 ( 均 差 ) 表 


0.25 0.5000 

> 1 0000 | _ 
0. 25 0 5000 0 954 0 > ~ 0200 = fm, nr 如 
0. 30 .347 了 二 ass , > 0.7193 一 并 rz rs) 
0, 39 9.624 5 co 7711 > —0.5440 = ffr:, xs» zx] 
0. 45 0 0708 < 07150 < 04050 = fr, mm] 

- 2 

0.53 0 7280 oo og6 8 025 = fs ws x] 
0.53 0.7280Q 


《1 7) 已 知 一 阶 避 数 边界 条 件 , 弯 答 方程 组 


28 灼 导 分 析 全 析 嬉 解 


2 1 了 了 
1 
5 日 nN, — 0.320 0 
和 2 3 
M, —0,7193 
2。 三 hf | 一 6] 一 0.5440 
5 5 
; — D0. 05 
El ， 于 Mi; [| 
7 7 7 MA — 0,352 5 
L 1 了 j 
解 此 方程 得 
MM, =— 2.027 8, 一 一 1.4643， M, =— 1.0313 
M; = DD. S07 2, M, =— 0. 6539 
三 次 样 条 表达 式 为 


1.878 3.7° — 2.422 ?x + 1,859 1x— 0,157 3,7€ [0.25, 0.30] 
0. B01 Bx — 1, 453 Br’ + 1. 568 5r+o.186 3,.r €E TO.30, 0, 39] 
DO, 622 Br’ — .244 Ox 十 1.4860r 十 0.1970E [0.39, 0,45] 
0, 319 47° — O. 834 Br: 十 1.302 5r+0.2246,.7r EE [0.45, 0.537 
(1 已 各 二 阶 导 数 过 界 革 性,Af = MM 二 0, 普 答 方程 组 


Sr == 


2 2 0 
1 nt -0,7198 
2 2 Z| lm| el_o50 
5 < :21 “? 
MM, 一 0.405 0 
0 3 2 
于 
解 此 方程 得 
MI -1,8809, M, =—0.86] 6, M, =— 1.0304 
二 次 样 杀 表 杰 趟 为 


— 6.269 7x° + 4,702 3 — D0, 205 9r+0.355 5, xtE [0.25， 小 4390] 
1.887 87 — 2.639 3 1.9985 6r+0,.135 3, x € [0.30, 08. 39. 
68 Hr +0,117 Br + 0.921 3x + 0,275 1, x Ee 20.39, 0. 45] 
2.146 Fr — 3.413 2x 十 2.510 3 十 站 0367， rel[o.45,0. 53 
21. 著 扩 人 和 人 [ay 45Str) 是 三 次 样 系 志 数 ， 证明: 


nh | 
{1} | Crcoldz- | [SCCry J dx =| [PCzy Srl dr 


SL) 一 
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2 SCL Fr) 一 SG]dz 
{ ii》 贰 FrzySCoG 一 0 11RD); 式 中 {zx,) ?0 是 梧 慎 节 点 且 有 a 一 


各 则 
[sor Dslde = SD Sm SF Ca) — Sa)] 
十 明 《1) 左 边 一 站 [fCD]az+| [S03 dz 一 
nh 
中 FTIS (ryder < 
| Crde—| [SCT dr+ 
?| ”Cr 全 Cr) Ft) de 
从 而 有 
六 nh 日 
| ECFenp] dr — [Esc ar = | TF) — Stridr+ 


日 
?| SLED SC Nd 


Ci) 左边 亚 > sodcren —S'(r)] = 
5 一 


-| 
PEF eS (全 1 一 
“9 [1 , or 
5| [fF (2 — SC]S"(r)dzr 一 
=n 
Sr tr oS tr 
Sn)LF tr oS (ro 
a 1 十 于 | 
2 (下 [Fx)— SCr}y ddr = 
此 二 11 


(Cor) 为 常数 ,2 三 《5 PE 
SEF Oo StH] 一 
SOF DT Soi+ 


Ss 2 j CFan ~ sco]。 = 一 
号 PEF S60] (ai (a) — Ste) 
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第 3 章 函数 让 近 与 曲线 拟 合 


1.1 重点 .难点 人 笃 析 


1.1.1 蚁 数 通 近 


1l. 概念 
对 函数 类 A 中 给 定 的 函数 了 (xr) ,在 男 -… 类 简单 稀 便 于 计算 的 薄 数 类 BB 中 
求 泗 数 (co ,使 ptz? 与 fix) 的 误差 在 革 种 度 最 意义 下 最 小 ， 
def 


(1) 最 圭一 致 儿 项 式 通 近 问 题 . 设 fr) 后 CT[a.i 求 站 (0r) 二 H, 一 
spanl]; yy ,ll 使 
| 7 一 站 1 一 rmin 工 和 一声 |. 
pn 


def 
式 中 17 一 名 下 -一 max | f(0) 一 ple)! 


<2) 最 佳 平方 通 近 问题 . 设 ftz) E Clay 中 , 求 SO em 于 


span{ gp CT} pr iC CLa, ,使 得 
LAS i= min | 大 一 全 
SE 
def re 

式 中 1 一 5 | [fn S00 Tptr)dz. 

(3) 最 佳 平方 达 角 表 近 问题 . 取 负 = 二 spanll, cosx, sinr: cos2r: sin2r, 
""} 的 最 佳 平方 带 近 向 题 . 

《4) 有 理 逼 近 问 题 . 用 形 如 
ef Pin 7 


i 
a (zi) La 


Dh 
一 让 


的 函数 近似 /4z) € CLa: #8], 使 得 1 了 一 民 .n tr) 上， 最 小 的 癌 题 称 为 最 佳 


Rn Cz) 


第 3 章 儿 数 娩 近 与 曲线 拟 合 31 


有 理 一 臻 逼近 ,使 得 f(z) 一 尺 。(z)1s 最 小 的 问题 为 最 佳 朋 理 平壤 逼近 . 
2， 笋 永 斯 特 接 斯 ( 凤 sierstrass) 定理 误 伯 轴 斯 坦 (Bernsrein) 条 项 式 
(1) 瑶 尔 斯 特 拉 斯 定理 , 设 fz) E Cie, 加, 则 对 任何 e->0, 总 存在 一 个 代 
数 多 项 式 pix) ,使 
| Az pla) | re 
在 La, 5 上 上 一致 成 立 . 
{27 拍 恩 斯 坦 多 项 式 


del 
Bf n= PFE)PAD, Pr) Ciao" 
EE=n 
收 基 C 性 。 若 tz) € C*[o, 1]; 则 limB™ (Ff, x = f(r) 


单位 分 解 2 Pn |= Dp) = (rr4r Ie 
让 三品 是 二 六 


0 过 rl 
稳定 性 | Btf, r) | max | Fr) |e | . 


1.1.2 范 数 . 内 积 


1, 范 数 
(1) 定义 : 设 有 定义 于 线性 空间 S 上 的 实 值 前 数 | .|| 满足 条 件 ， 
461797 上 宇 人 , 当 有 上 且 可 当 x 工 = 人 时 zj = 0;{ 正 定性 } 
(231arl 一 ai zl y{ 奇 次 性 } 
371xr 二 + yy 过 上 cH 十 上 vy;( 三 角 不 等 下) 
则 称 上“ | 为 线性 空间 S$ 二 的 一 种 落 数 ,5 与 上 ，| 称 为 赋 范 线性 空间 . 
《2) R" 上 的 常用 范 数 ， | 
I xl 一 Tmax .7 152 - 范 数 或 最 大 范 数 ) 


zi 二 2 1 1;(1 - 范 数 ) 
EE 


1 rl = (D1)7.02- 范 数 ) 


(3) CLa. 如 上 的 常用 范 数 ， 
| Fe = max | Az) 1sCo2 -基数 ) 


Fj -| | flxr): dr(1 - 苑 数 } 


32 数 利 分析 全 析 精 角 


171 = (Podr) .42 - 范 数 ) 
2. 内 各 
(1) 定义 : 设 久 是 数 域 臣 { 朋 或 Cy 上 的 线性 空间 ,(*,， +} 起 定义 于 XX 站 到 
K 上 的 二 元 烟 数 , 且 满 足 
共 罗 对 称 性 Yu,wE 外 有 (us v0) Cv ii 
线性 组 合 性 Yu vw 世祥, a1 四 民有 Lan 二 Be ww) = eu tv 
ws 
正定 性 YuE 各 有 to 汪 0, 当 且 仅 当 w= 0 时 (tu, 站 一 0 
则 称 (. ,+") 是 苇 上 的 一 种 肉 积 . 定义 了 内 积 的 钱 性 空间 称 为 内 积 空间 . 
C2) 正 变 : 鞠 内 积 空间 中 的 元 素 #, + 满足 C(x, 中 二 0. 则 称 wu, ww 关于 内 积 
(,*) 正 交 . 
(3) 柯 西 - 施 瓦 欧 (Cauchy 一 Schwarz) 不 等 式 ， 
Lu au》 | RE us HC 
《4》 定理 3.1 内 积 空间 关中 的 元 素 列 4 ;人 线性 无 闫 的 充分 必要 条 件 是 
格拉 姆 (Gram) 矩阵 己 非 亲 异 , 即 dettG) 关 0, 这 里 
[| 


def | Eu Wa} Ciar Wa) re (Cua Wy 


Cus Mi as Mn) oe 


(5) 内 积 诱导 范 数 Cu, 4)Y. 
(6) 以 上 关于 肉 积 的 结论 对 于 带 权 内 积 也 成 立 . 


i.1.3 正 交 多 项 式 


1. 定妆 
多 项 式 序列 | pu Cr)}2h 如 果 满 足 
def rs D， jE 
(opis ps | ping rr dr 一 [A [| 


则 称 {p, Cx) i1 是 [a, 如 上 带 权 pz) 的 # 次 正 交 多 项 式 系 . 这 里 gr (zr) 是 首 项 
系数 不 为 惟 的 nn 次 密 项 式 tn 二 90,1, "**》. 

2. 正 变 多 项 式 系 的 性 质 

0]) 线性 无 闫 性 :和 中 2 中 的 任意 有 限 字 序列 均线 性 无 关 ; 进 而. 任何 不 超 
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过 于 次 的 多 项 式 均 可 以 用 { 和 人 -线性 表 出 . 
《2) ga 与 所 有 小 于 并 次 的 多 项 式 正 交 . 
(3) gr 在 fut 有 ) 内 有 点 个 互 不 相同 的 实 零 点 . 
(4) 当 1m 1221 是 首 项 系数 为 1 的 正 交 多 项 式 系 时 有 如 下 递 推 关系 


pr = Crag —phg, n=01,. {3.1) 
式 中 
和 1 


《or fr) (gn 1 gl 
3. 斯 蜜 特 (Sehtnidty 正 变化 方法 构造 正 变 多 项 式 系 


li DD ey,, no 1, 2, 
4 常用 正 交 多 项 式 系 


{1) 勒 让 德 (Legendre) 多 项 式 是 [一 1 + 1] 土 带 权 甫 数 ptxzr) 二 1 的 正 交 多 项 
式 系 . 定义 为 


__ 1 中 :es is， 加 ， 
Ptr) Pn) gn a 1)7?)， A= 1, 2, 
并 具有 下 列 主 要 特征 ; 
了 一 2 二 EL 
DP PY) 0 1 日 


《iD 为 偶数 时 PCz) 为 偶 函 数 ,n 为 奇数 时 已 , (xz) 为 奇 函 数 ; 

CH) Pt) 在 (一 1, 1) 内 有 个 不 同 的 实 零 点 . 

(iv) 首 项 系数 为 1 的 甚 让 德 多项式 了 , (x) 是 所 有 首 项 系数 为 1 的 n 次 多 项 
式 集 人 台中 在 上 外， 1， 意义 下 距离 零 最 近 的 元 素 . 


(2) 切 比 雪 去 (Chebyshev) 多 项 式 是 [一 1, 本 上 带 权 画 数 ez 一 一 1 


一 


揭 正 交 才 硕 式 系 . 定义 为 


T(x) = cos(narceosr), HH = 0 1 
并 且 上 有 具有 下 列 主要 性 质 ， 
(0) T(x) 最 高 次 项 系数 是 2 一 (0n 1). 
ti) CT TY) = x (TT,, TT) 二 六 Cn 1 2 
(iii n 为 偶数 时 工 ,(x) 为 侦 函 数 ,n 为 奇数 了 时 工 ,{.r) 为 奇 消 数 ， 


(iv T(z) 在 {一 1:1) 内 的 有 个 零点 a = CDS 人 


Tr 
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在 及 = cos 所 全 = 0,1，…。) 点 处 交错 取 最 大 值 1, 最 小 值 一 1 
0) (x9 = TirTitz) 是 所 有 首 项 系数 为 1 的 4 次 多 项 式 集合 中 在 
|。 1. 下 距离 零 最 近 的 元 素 . 
1.1.4 最 佳 一 致 逼近 多 项 式 
1. 定 头 


def 
(1) 偏差 AF, PY 一 P|, Fery EE Ca Hl, Pry € H.. 


def 


{2) 最 小 偏差 上. 

3) 最 佳 -一致 逼近 多 项 束 . 如 果 Pr Lr) € 万 ,使 得 Atr P= 二 EB,, 则 称 
Pi tr) 是 x) 所 Ca 在 La, zl 上 的 最 佳 一 致 通 近 多 项 式 或 最 小 坊 差 通 近 
多 项 式 ,简称 最 佳 通 近 多 项 式 ， 

2， 存在 性 及 唯 -- 性 

定理 3.2 若 fx) © Cia;5,; 则 存在 唯一 的 寥 项 式 Pi (x) 所, 使 得 


| A) Co Pi (Cx) tb . = EE,:p 
Mtf; Pi1+= EE, 


i AC Ps}!}, 
PEH, 


3， 女 比 雪夫 定理 
(1 定 吕 : 若 中 (zz 一 Fe 下 Pr 一 天 xz 称 x 是 瑚 zl 的 正 偏 差 
点 ; 普 忆 一 兵 加 一 Ps) 一 下: 称 过 是 Pr 的 负 偏 莽 点 ， 


(2) 定理 3.3 PCr) EE 万 ,是 jz] E CLa, 如 的 最 住 一 致 带 近 多项式 的 充 
分 必要 条 忻 是 Ptzx}) 在 [a, #5] 上 至 少 有 # 十 2 个 轮流 为 正 \ 仙 的 杞 差点 ,中 有 十 
2 个 点 各 小 使 得 
Peri— Fm = C1 Pin Ftryl., sg=+l 
这 桩 的 点 组 称 为 切 比 雪夫 交错 点 组 ， 


un 
4, 在 万 ， 2 Ee Yar', 2 写 RE. 0 所 i 所 nn 一 1 上 的 服侍 一 部 远近 
定理 3.4 设 Fzrl)E ctl,11, T(x) E 日 , 是 f(x) 在 宫 , 中 的 最 


佳 一 至 过 近 多 项 式 .偏差 为 如 


5， 最 佳 一 次 适 近 多 项 式 
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设 扩 mi) 二 [a ,的 上 且 天 (z) 在 (a, 的 内 不 变 号 , 则 六) 的 最 佳 一 致 一 次 


通 近 多 项 式 为 


十 - 
Pitz) 一 Efe + flr)] ta Ea 5) 


m= ED -rz 


bs—a 


变 错 点 组 为 :Ts 和 三 
1.1.5 最 佳 平方 通 近 


1 f(x} Eee[ay tl 在 而 一 Span{ gus Pl ~ 5 全 [ap 的 上 的 里 住 平方 


衣 提 SS=asg 十 qigs 十 十 angpn: 了 票数 {0} 是 法 方程 组 


| Co) {ge pe) an tf, ph) 
def ip Ph) Cg pi) ph) a fg 
ei, Cpr Ph) (gs pn) a Cf, gn) 

的 解 向 量 ,平方 误 善 为 


- def , 2 
| cry Hs 站 一 写作 一 | Fl Oatf, pe) 
式 中 内 积 和 2 - 范 数 分 别 定义 为 
(fg) = | fgpdr, lf = (Ff, 12 


2. 基于 正 交 基 的 里 性 平方 表 近 
设 fir) EE C[a， ,本 数 繁 汪 有 正 交 基 { gp 7-，. 
(1) 表达 形式 。f 在 鲁 中 的 最 佳 平方 逼近 


上 


. a) 
We 


平方 误差 


。 4 , 
lli= IAl— Datpe P= NAN Ym 
i Py | Fx | 2 

tz) 贝 塞 尔 (Bessel) 不 等 式 


Daiinlis I: 
k=0 


3， 广 六 情 黑 叶 (Fourier) 级 载 


(3, 27 


{3.3) 


《3. 4) 
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设 [p;} 训 5 是 区 间 fa， 引 上 上 权 国 数 为 pxry 的 正 交 函数 系 , 且 有 1 入 1 灾 
CLe, 六 ECLa HJ: 称 


ee - > Co ep 
为 子 的 广 祥 傅 星 叶 级 数 . 


4. 勒 让 德 多 项 式 B, 是 堆 函 数 在 严 , 中 的 最 佳 平方 逼近 (ofz) 三 1 一 ] 势 
TS 1). 


1.1.5 最 小 二 乘 井 线 拟 合 


(1) 线性 空间 中 的 最 小 二 渠 曲 线 拟 澡 问题 . 已 知 数 据 计 ， 姑 盖 四 和 在 
更 = spanig， 的， 加) 中 求 一 函数 S' (xr) 使 得 
(二 LS Ci) — fr 


jh 


最 小 . 这 里 iretzr,)}7s 是 权 系 数 , 均 大 于 零 . 
(2) 最 小 二 乘 解 5(z) ~ Pay 的 条 数 满足 的 革 方 程 组 的 形式 交工 (3. 2)， 
只 须 将 其 中 内 积 换 为 如 下 离散 形式 的 内 积 ， 
(fg) 一 SiC) fr gl,) 《3, 6) 


{3) 存在 唯一 性 . 当 { 和 ur 六 ,甘于 点 上 集 iz,17%;, 洽 足 略 东 (Haar) 条件 时 , 则 有 
法 方 积 组 存在 唯一 解 向 量 

陪 尔 条 件 : {gr Ms 忆 Cra 说 的 任意 系数 不 全 为 0 的 线性 组 全 在 点 集 
i un 谤 2) 中 至 名 有 个 不 同 的 罕 虑 (gr ji-, 是 线性 无 关 的 函数 组 )， 

C1) 用 关于 点 集 的 正 交 多 项 式 芍 最 小 二 乘 的 合 . 关于 点 沫 正 交 多 项 式 构 造 
方法 同 式 (3. 1) ,只 须 将 其 中 内 积 换 为 离散 内 积 式 (3. 的 . 用 该 正 交 多 项 式 构造 
最 小 二 乘 拟 侣 解 形式 网 式 (3.+1. 


1.1.7 基于 三 角 函 数 系 的 最 仁平 方 通 近 和 最 小 二 彝 拟 合 
1. 县 性 平方 各 和 近 

设 Flz) 是 以 2x 为 周期 的 平方 可 积 函 数 . 

(1 fir) 在 [0, 2r] 上 的 最 佳 平方 通 近 S$,tr) 为 


Sule) = 六 + DD arcoskr Lbsinkr) 
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式 中 情 里 叶 系 数 a 定义 淤 


a 一 二 | foeoskedr, R01, 2 
Tu 


入 一 {Ansinkrdr, 大 一 
mn 


《2 fA) 的 全 里 叶 级 数 


1, 2 


Fir) ~ 六 a + 之 (ar coskr 十 hsintr) 


(33 当 产 {x) 在 [0, 2 
fT 二 人 二 Fr 一 0 


2xj 上 分 段 连续 时 , 则 有 


] 


2 
C4) 平方 误差 
1f—s. 人 一 站 
(5) 了 下 不 等 玉 


2 


i 


na su 十 Darcoskr + bsinkz) 


| 二 + a+)] 


a Patto) LIA 
T 


2, 本 
(1) f(xr) 为 实 值 时 : 
11， cosr, 


sinr, cos?z, 


i 2 | 2 


1 


23m1 1! 正 交 
fz) 的 最 小 二 乘 三 角 垦 近 为 


St7) 一 5 


式 中 


Dc, 


_ 2 
2m 寸 1] pe 


bh; 一 


2 
2 对 十 1 会 


特 其 地,: 当 二 澡 时 ,Stx) 是 关于 fi 


(2) Frl 为 复 值 时 ， 
已 知 六) 在 | =- 7 ) 


慰 拟 合 


ees 


™ 1 


j= 


sin2r， …， 


Lo, 十 > tarcoskr 十 hy sinkzr } 
一] 


ps, ] sin 


Cosptr: sinmr! 


[| 


i KO len 


2mjik 
2m 十 1” 


5 的 三 角质 值 西数 


下 一 1.2,. .nn 


上 的 值 {7, = (等 


关于 点 梨 
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和 关于 太一 和 | 是正 交 的 . 


Atz) 关于 1 3) ) 的 最 小 二 冬 通 近 为 


ot 
Stz) = DOer, nN 
业 一 全 
1 局 2 
式 中 wi k=0,1,.,N—1 
; 2x. 下] 
特别 地 ,n 一 NN 时 Scix) 是 关于 {3;】 ”的 播 值 , 即 


= Sie, es, j= 0 1 N—1 
式 (3.8) 和 式 (3.9) 分 别 被 称 为 离散 博 里 叶 变 搞 和 反 变 摘 . 
1.1.8 帕 德 (Page] 逼近 


1 定 浆 
度 站 Xx) EC 一 aa a)， NN 二 # 十 所 ;如果 有 理 了 两 数 


ta Pr) 
Rut ad 二 十 gr PP, 
(一 


其 中 Pt 与 @, tx) 无 公园 式 , 目 满足 条 件 
RE 一 FN, k=o lr N 
则 称 Ri tr? 为 晴 数 fx) 在 = 二 个 姓 的 tn mr) 阶 帕 德 遂 近 ， 
2. 套数 inf 的 确定 


《3.7 


{3,8) 


{3.97 


引入 符号 < 二 0 j = 0, 1, 2,… ,十 mm. {7 人 是 方程 组 坊 一 5 


的 解 ,其 中 
Fi -or Cl Ca 五。 
Cn -ere er Ct - be | 
亚 一 p= ce 
er ar 1 
当 j 了 之 0 时 约定 c= 一 0, 二 1, 当 j 守 mm 时 b; 二 人 0 
二 -1 
fa 由 ex 一 DB, 十 上 炳 定 . 
jh 


ttm 


3, 月 理 函 才 的 画 载 恒 可 转化 为 连 分 式 函 载 愤 的 计算 ,这 样 可 以 节省 乘除 法 


需 3 草 ” 消 数 滔 近 与 曲线 氮 合 


39 
的 运算 次 数 . 
4 截断 误差 
人 def 交 。 
Fr Rn (ry = Qa (Cx) 天 抑 me 
当 | 所 1 时 误差 近似 表达 式 为 
def 这 
FT Re a at, po et 
由 = 让 
1.2 习题 全 解 
]1. fr) = sin BR 给 出 LI， 1] 上 的 伯 轩 斯 坦 和 多项式 Btf，I) 及 
Bt ff, ry. 


解 。” 伯 感 斯 坦 包 项 式 为 ”Btf. x) = DE) a 
当 wn 一 1 时 ， 


1 1 
五 区 六 全) 一 fo jaal 一 卫 ) tAD{ ja 一 了) 一 


tsinF XIXrxl—x 


ont ) a + 
OF | 3 14 


Ee 


Osin X31 


Te 


sn I 十 


sin 于 XI1Xxmxl=m53yi 6 3 
2 2 2 
— 0.098r: 一 日 402r2 +1.5r 


2. 当 六 I) 一 工时 ;求证 BB (ff, Ir) 二 渡 . 
证 明 车 fry = 二; 则 
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2 9 nn 
Btf, Xx) =>) 去 亿 je 一 x)"* 二 >) 二 人 和 02 


二 一 中 二 一 1 


3 肯 nin Dn kD x 
此 一 】 n 


1 


La Dr t = 上 一 1 
[> {nO 17(n 一 2)r(n 一 目 二 Dm on 
t= 


{一 1171 
『 3- DY 了 
| 二 {ne 1)tnoe 2) 人 1— x+ Dd | 
| m! 


[ri itl x] r= 
3, 证 明 函 数 1， 工 ，"…， x" 线性 无 关 . 
证 阴 只 须 证 明 如 下 等 式 : 
如 十 电工 十 GT -十 ar 一 D WT7TE 民 
只 有 零 解 @ = (ay al ,a7 二 从 为 此 ,分 别 取 tk 二 0,1. 
式 商 端 在 [0，11 上 作 带 投了 肖 数 ef) = ] 的 内 积 , 得 


[1 1 革 1 

: 2 3 nn 十 ] er 1 r 

i 
2 3 4 再 二 了 us |= 1。 
1 1 1 ,,, 上 L | 

nil n+2 2 十 3 29 上] 


，8) ,对 上 上 


该 廊 程 组 的 节 数 年 阵 为 希 尔 犀 特 窜 阵 ,对 称 正 定 非 麻 异 . 故 只 有 零 解 a = 0 现 


此 函数 lx x" 线性 无 关 ， 


4. 计算 下 列表 数 天 2 在 EL 1 上 的 下 与 上 天 | 


四 


fi 一人 一 2 ftr) 一 


| 一 


|; 


{HFCrY 2 一 "mn 为 正 整 教 ; 
4 = rt 1 "er. 
解 (rr 一 3 一 1 0reEts 1). 故 fizr) 单 增 . 
了 na | 一] 一 maxtl F007 |, i fly 1: -= 1 
:= | Fr) dr =|.0- mdr 一 本 
+ 


1 fF" 一 站 rooa -- i rd] 万 
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1 
= 村 | = max{ 1 ro /( 志 ) ;| f(D 1 二 地 


上 业 一 上 , 
| 六 站 < sn ly .| 1! 
1] 一) | = 
CHF EY = rm mt rr E Los 


信 六 (z) 一 0. 得 驻 点 -外 f(=) 一 (和 )】( 轴 和 5)】' 于 是 


全 
四 i . #4 i Hn 
f! = max| FOO f(D, f(s) - 
! " ' ! ml a minl! 
FF | | Fr rr |: [中 mi 


iFi = [Fe z | = [| 【1 qz | 一 


CF Ee = 10tr+ lMe itr DD itoe)™ 
0 在 [9.1] 上 单 增 


Tm 2m 12n)! FP 


f= max Fr (= maxi| ft0Y |, All) i 
1 


| Fl; = | i fr) :dx 一 [et 1 "ed = 


1 
-re 十 | otr+l Ye dr=3— 15 


Tp 1 : 
i fF: = | De :rdr, = 了 3 一 二 】 
5. 证 明 ‖ 了 一 8 蕉 7 一 TH 
证 明 一 7 一 有 过 一 
所 以 7 


6, 对 了 x) gr) 二 Cfa., ,定义 
(1) tf, g} = [fn nade 
《2) CF k=- | feng rdet fa ga). 


问 它 们 是 香 构 成 内 各， 
证 明 ty 取 Ac 一 1 天 DRzrlE Cia; 办 ,由 
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人 7， 门 = jf DF dr =0 
知 正定 性 条 拍 不 满足 ,也 即 相 构成 办 舱 . 
(2) 正定 性 《六 万 一 「 LF Odrt Fra SD 


车 tf,， 让 = 癌 则 有 | [fFar =0 和 和 fla} = 0. 


[Wd = 0=>f (x) = 0 Om, | > 
a . Ta) =0 
正定 性 得 证 . 
对 称 性 (fr 一 【和 六 7 
齐 次 性 laf 8} = atf, 8) 
是 显 的 成 立 的 (ea 为 实数 ). 


设 ,gE Cia ti, 则 


{fig:h) | [f+ gh dri[Fit) + eta) lhta) 一 


| Fh dri+ ftayhta) +| ER'dri+ gtalhta) := 


Cf, hy ig: AY 
综 台 以 上 + 点 知 定 义 (2) 构成 C'[a, 直上 的 内 积 . 
了. 令 Ti7 Cr 一 TC27 一 1, 了 EE [0 1 , 试 证 1T Cr 是 在 1D， 1 上 带 权 


二 一- 的 正 灾 多 项 或 ,并 求 TT? (x2), TY (7), TT Cx) Ty (x). 


ptrx} 一 


了 工 一 工 - 


1 1 1 
证 明 [ 了 丁丁 Crypl ry dr =| Ti27r— TT, tr 1 — dr 


TK 一 
5 二 0 


. 1 
故 1T Cr :i 是 在 [6.1] 上 带 权 pir) 一 一 -一 一- 的 正 交 和 多项式. 
， 出 2 I 一 zr 


了 = l= T(r} = 1 
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ry rm Cr) := Srl 
2t27--1) 1=8ri— Br+l1 


Tatr) = 27 ~ lw tr) 
Tr = 4 rT (r= dr 1 — (2r—1)= 
327* — 48x? 十 18r 一 1 
8. 对 权 函 烤 交 一 区 间 [ 一 1]，1.J, 该 求 首 项 系数 为 1 的 正 变 多 项 址 
nat) A Or 1, 2 3. 


解 在 区 间 [ 一 1,1」] 上 定义 内 积 (ff, g) = [ rigtriot rdr 


一 1 


HT 


~ Cp po) J 0 _o 


Fr} 一 1， a Cp pe) B73 


Cm py) 0 tp pi) 
mT pp) 6s 0 8 B73 一 


| 


mr rp 


a = {Crp 0 一 60， 此 = to Pe) 13673525 17 
tp gp) 136:525 (gm gi 了 16715 了 7 


I) 二 下 


9 该 证 zy 一 sinf Cn 十 1)atrecoszr] = 六 
访 证 明 ws (7) A 给 出 的 第 二 类 切 比 雪 去 多 项 式 闵 


tt tE， 是 区 间 [ 一 1，1] 上 带 李 pa 一 1 一 x 的 正 交 多 项 起. 
[ sine + 6sinCm + 1)898 


1 B= NS 
证 朋 | ustIIun CTI oT} dr es 
4 

def 

= lj, 1m) 


当 n = 加 时 ， 


Tn my = 中 1 cot t+ Dag Ey 


当 n 了 mm 了 时 ， 


i ， eostm+ la [~ costmt 100 了 了 
Tin ny | sin 上 了 ed 下 m1 sinta 178 ， 1 


] costn + Dirostm + 1 df =: 
ju: 
sinmt+ De 


nl S 
I Sosa 178d | 
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并 十 1 - _ 下 
[= HT 
(2 ) | sincn + DosinCm 一 1)6d8 = (Ei) lr. m) 


于 是 Ion mw)| 1 一 (21 ) j= 9, 进而 Ka m) = 0. 正 交 性 得 证 ， 


m1 
10. 证 明 场 比 雪 志 多 项 式 T(r) 满 是 微分 方程 
(AT) Tr T= 
证 明 切 比 雪 专 多 项 式 为 了 fr) 一 cosfrmareebsryy :了 | 过 1] 
nt 


wr 


Ti tr} = sin(narccos.r} 


a nm: 四 
TTY -2 Costnarceos.r) 1 5 一 SLC aceDg 
一 ”了 


nT 


i | 二 + 
1 一 23 T, tr} 1 一 六 rl, tr) 


两 边 同 莱 以 ] 一 xz .并 移 项 得 
TT T(r) 0 
11. 要 证 fr 在 [a, 上 8 上 连续 , 求 fir) 的 香 次 最 住 一 致 通 近 多 项 或 ， 
解 由 闭 反 间 这 续 函数 性 质 邯 ,存在 .mr E ai 使 得 max fC = 
Ar 和 min fr) 一 flrs) 成 立 ， | 


取 忆 二 二 5fCn) 一 (xz)]; 则 有 


max zr— Pl= ma Fn} -Pl | tr pl} 
而 左 


Fe 一 已 = 一 (1 一 下 一 天 2 了 天 


用 上 了 一 下 = Ar 
进而 有 - 和 是 卫 通 近 了 (zx) 的 两 个 交错 的 篇 车 点 .机 雪 比 雪夫 定理 知 忆 就 是 
AT 的 零 次 明 佳 一 禾 通 近 多 项 式 ， 

12. 选取 常数 41' 主 max rar | 达到 概 小 ,对 问 这 个 解 是 香 哈 一 ? 

解 x 一 ar 是 [一 ]，11 上 的 奇 郁 数 , 故 
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ax | —ar!:= max | 于 一 az | 一 fr 一 cr) 一口 放 - 
娄 使 上 式 达 到 极 小 ; 妈 求 5 在 [一 1. 11 上 的 三 次 最 佳 一 至 通 近 多 项 式 . 由 
切 比 雪 去 多项式 的 性 质 知 


Ti—ar Tt ry 3 【43 一 Sr 


时 上 并 一 gr 一生 / 达 到 最 小 , 帮 a 一 六. 
由 最 佳 .- 致 逼近 多 项 式 的 叭 一 性 知 a 也 是 唯一 的 ， 
13. 求 f(z) = sinx 在 | 0 要 | 上 的 最 佳 -次 逼近 多 项 式 , 并 信 计 误 关 


解 二 次 导 晒 数 Cr) = 二 (sinr)” 二 -sinzt 当 x [9， 六] 时 不 变 号 .页 
六 xz) 的 最 佳 -- 次 鼻 近 多 项 式 为 


Pt) 


f(t 
FO0 + Fr) 2 0+ x 
他 这 T {x )= 


于 js) 三 (= 52x ) 


A 和 A000) 


L 
5 骨 


Fr ,BI cos == 二 , 求 得 
x nt 


Ti =: ATCOS 二 2 B80 9 rl 一 sinf arceos 全) =- 人 ,了 ?1 18 
™ Tx 
将 以 上 数据 代 人 Pitx) 得 
Ptxry = 三， 十 日 105 26 
x 


误差 为 


| sinxr 一 PCr ll 一 |sin0 一 Pt = 0.105 26 
14, 求 f(tx) = 在 0. 11 上 的 最 性 一 次 逼近 多 项 趟 ， 
解 F(t) 一 在 0, 1] 上 不 变 号 , 它 的 好 仁 一 次 逼近 多 项 式 


yo Or Fo fi0) D+ ra 
P, (ry 5 = 地 {x 5 ) 


1—0 


] - . ~- ] 
Lilt rt ce—1}{r~ 二) 
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xz: 满足 方程 fre) = A 
即 be 一 上 上 一 上 


解 之 得 re = nfe 一 1),f(r) 一 er 一 e 一 1. 故 有 


PoD 一 全 十 te Dz sine » |]= 


(e Dz 二 了 fs te— llnte-- 1)] 


15. 隶 fz) 二 区 十 87z? 一 1 在 区 间 [0, ]] 上 的 三 次 最 佳 一 笋 讶 过 多 项 式 . 


解 令 :=2(z 一 全 ), 则 re [1,1 
f(D0 一 (等) +3( 等 1) 一 1= 


To 和 加 | def 
Fa 10r 2 T2219) gitt) 


16g() 是 首 项 系数 为 1 的 站 次 多 项 式 , 记 它 的 三 次 最 佳 -- 孝 带 近 路 项 式 为 


六 0 这样 
16 人 一 Pt i, = | [lg — Ptryd—o).. 


这 到 了 最 小 . 利用 切 比 雪夫 和 争 珊 式 的 性 质 知 
16g( — Pe 一 了 ti = Ler 一 8 十] 


1 


书 ( = 16g(0 一 于 (9 二 109 二 25 十 32 一 宣 


ki) 的 最 佳 一 臻 记过 多项式 为 让 Pt0 ,进而 (x) 的 最 佳 一 至 逼近 多 项 
式 为 


工 一 = 二 | rl1y: 六 3 gn 全 |= 
Bh: 1) 6 02 一 十 2552 人 一 1 十 2 一 ]) 8 
5 129 
Br 十 地 585 


16. f(x) 一 工 | ,在 [一 1,1] 上 求 关于 四 二 spani1， 17 ，r!}) 的 时 仁平 方 通 
近 争 项 起 ， 
解 内 积 (f 四 一 | fgcndr 记名 i 


2 2 
(om i) = 2, (om 二， (or 和) 一 
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47 
一 一 了 = 二 
(pp 二， (wp; Fe) = 7 to ps) 9 
(fm ll, fm), (frm 
,oh ， ， Eg + pe 本 
求解 法 方程 组 
2 2 
了 村 三 ， 1 
BI 
2 2 2 _ | 至 
3 
2 2 2| 2 |3 
5 ? 9 ' 
和 15 105 105 
:二 一 了 5， 一 一 AR 1, 25，a。 一 一 了 Se 
得 a. 128 日 .1]17 187 5 人 4 ] .64D 625, a 1 呈 
— 0. 820 312 5. 


故 F(xr) 的 最 佳 平方 逼近 密 项 式 为 
a tag Tag = 0.117 187 了 二 1.640 625 x — 0. 820 312 Sx 


17, 求 涡 雪 f(r) 在 指定 区 间 上 甘于 外 一 spanll， x+ 的 最 性 平方 通 近 多 
项 式 : 
(1) fex) 一 二 ， r€ i, 3]; (23 tr) = er Ero 


(3) flr) = cosxx, x ECO, 1]; CAFEr) 一 nz 和 12]. 


解 (1) 内 积 (f, g) = | fengcndr 由 


CD) 一 2， xy 一 4， (rzr) 一 剖 


Cf 1) = ln3, Cf, rry 一 2 


4 An In31 
解法 方程 组 | 。 26 | 1=[ ] 有 
4 一 I 2 
3 
a 一 ln3 6 ~ 1.141 0 a = 3— 3n3 0.2958 


这 样 求 得 f(x) 的 最 佳 平方 通 近 多 项 起 为 
pr) 一 如 XI 十 ai XT=1.1410—0.295 8 


【23 内 积 ( 六 8) 一 | 
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gr) 一 10 一 故 有 


， ，_ 1 1 
站 
tf) el (fm = 1 
1 
mua | 
i 1 | La ] J 
2 3 


da; = ,187 8, cw = 1.6244 
喜 求 得 六 x) 的 最 佳 平 方志 近 者 项 式 为 
tx) = nar 二 0,187 8+1.624 4x 


(3) 对 f(z) 一 cosxx, x E [0, 1] 作 线性 变换 < = 方 十 记 4， 即 


2 


于 ] 
Fix) = cosmz 一 cos{ x) 2 Ee), (EeE [1,1] 


利用 勤 让 德 正 变 密 项 式 如 人 二 1， pitt) 一 1 为 基建 立 g (2) 的 一 次 最 佳 平 
方 允 近 名 项 式 


5 P Sp D+ EA p 


0 ip, pi pn 一 


ir) 一 一 过, 
抱 


他 PCD 十 2 


六 x) 的 最 佳 平方 通 近 为 


F211) = 【2 工 一 下 2,431 T+ 1,?159 
4 
(4) 内 积 { fy) 二 [ fingtrdr: 则 gtr) 二 1.gi(r) 一 z, 栈 有 
ln-llig li 


(fig) 2ln2 一 1 fp) er 22 一 二 


[1 | ， 2In2 一 1 3 
方程 | . |= 
解法 方程 组 四 加 | ] | 
2 
a =—0.637 1, a = 0.6822 


故 求 得 六 zy 的 最 佳 平 方 遇 近 多 项 式 为 


第 3 和 草 函数 盟 近 与 曲线 拟 合 和 


Pr = a ar 0.682 sr 0.6371 
18. fr? 一 sin 本 <, 在 [一 1， 1J 上 按 勒 让 德 多 项 式 展开 求 三 次 最 性 笠 方 通 
近 多 项 式 . 
解 记 ip, i 总 , 为 勒 让 德 正 奖 多项式 


T 
(Cp, Ps = | Pcidr 一 区 十 


1 
8 | 
Cp =0, (pf =, (hp, f=0 


.1 = 1 10) 


x) 的 三 be 页 式 为 
a Cp , 188( — 10) ,00) 
站 放 ot pr 站 十 1 psi {x) 


120(31 — 2x:} 420{n -— 10) ; 
十 I A 


上 


1 


Ci TT 
1. 53533 2.r -0.302 2r’ 
19， 观 . 铀 蝎 蛋 的 直线 运动 ,得 出 以 下 教 据 : 


地 运动 方程 . 
解 经 措 图 发 现 !+ 和 :* 近 和 伺服 从 线性 规律 . 故 作 线性 赚 型 ，= a 十 bt, 令 
更 一 spani1. t+ 计算 离 蕉 内 积 有 


(1, 1) 一 PD =6. (i.1)= Yu = 14.7 
一 人 


(Ds = 80 GH- Yt; — 1078 
求解 法 方程 组 得 


T 台 | [2 = 00, 
L14.7 353.63. Lp| Liorsl 


#5 048， = 22.253 26 
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运动 方程 为 1 一 一 了 .855 O48 二 22,253 ?76 


平方 误 盖 = Ys 2,1 Xx 107 


j= 
j= 


20. 已 知 实 验 教 据 如 下 : 


用 最 小 二 来 法 求 形 如 w= a 十 Mr! 的 经 验 公 或, 并 计算 均 方 误 差 ， 
解 下 = spani1l,x .计算 离散 内 积 


+ 
(1 = Dl=5, (x) zy 5 327 
了 上 一生 癌 


(ez 一 一 7277699 


jn 


+ 
(=o D2, Cr yo Dry, = 369 321.5 
41= 


解法 方程 组 
「 5 5 327 2] 271.4 
L5 327 ?7 277 5 = [00 321, | 


a 0 972 579, = 0.030 035 


得 


] 


均 方 谋 差 12II(0zi 一 四 卫 上 se0.1226 
21. 坦 某 化 学 反应 中 ,上 由 实 验 得 分 和 远 物 浓 关 与 时 间 关 系 如 下 |; 


时 间 并 


滞 度 了? 
(YX 104 3 


站 


用 可 小 二 蓝 法 求 了 一 fn), 
解 ”观察 所 给 数据 特点 ,建立 拟 合 模型 y 一 ae (a,p 0). 该 模型 关于 参 
数 非 线性 ,两 边 取 对 数 得 


jny = ]na 一 一 记 lna 二 六 


名 一 span|1, 一 地 ,由 于 + 一 0 时 条 件 自然 满足 ， 内 积 (f，g) = 
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11 


> FoO8t0 3 经 计算 有 


四 


.1 一 1，f, 一 十 ) = 一 0.603 975 


(一 了 二， 一)= 0.062 321， (C1, iny) 一 - 87.674 095 
{ 一 证， Inw) = 5.032 489 
解法 方程 组 

FF —0.603 9751 FA Fr—87.6740957 

L 一 各 603975 0,02 321 | |, 一 | 489 
得 A 一- 7.558 T87496 163， ua er 有 215103X10: 
得 到 拟 合 模型 

21 03 HY X 10r 

拟 合 平方 误差 


= DTy) ,3 910: 
1=] 
22. 给 出 一 张 记 录 i 大 1 一 人 3，2， 1 人 0, 12,3), 甩 FFT 算法 求 i 的 
离散 谱 1ere 


解 区 二。 区 一双 一 一 1 af 一 全 这 
计算 过 程 如 下 : 
kh 加 ] 2 3 t 6 了 
和 3 2 上 1 2 3 
六 | 4 4 4 2uw， 站 4 一 3 
A | 8 4 4 8 of 1 -2 
Be = A | 16 4 二 2 0 43 0 4-2Y 0 41+2 


[的 丙 基 语 是 (2.2 在. 0 2 0 于 在 0 2 并 ) 


， 本 3 三 十 Br 本 
23. 用 辑 转 相 除 法 将 Ry: Cr) 去 干果 下 化 为 连 分 式 ， 
解 
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ul 
[3 


3z 6r 0, 12x+ 18 4 12 
xX 二 Br 一 6 二 br 十 日 了 十 4.5 一 


Ras {tr) = 


24, 求 站 5) = sinr 在 工 关 人 处 的 (3,3) 阶 帆 德 通 近 Rs (7r), 


、 ;十 qT 十 HT 十 a 
Rs x 1 2 工 小 
解说 Rstz) nT Th 


由 fry = sinr 一 工 一 十 于 一 到 十 …, 得 
5 了 1 


31 
二 lO rs 0 CE 0 
{1-1 是 下 列 廊 程 组 的 解 
tl rs Cs 内: ey 
一 |c Ca oe 由 | 一 | 
3 1 让" 
1 
1 0 一 到 
6 Fh; fi9 
一 一 = | .1 
即 9 一 ?| = 1256 
_l 2 Ll* Lo 
6 120 
了 1 
得 A # 0, FF; 
解 ， ”20 
确定 fa,1-， 
| 


3 二 一 一 高 


ds 二 十 r= 0， 


故 有 
?os 
. : 3 
Hr ATA dr 的 0 ?7 


开本 二 -一 一 一 一 一 一 一 一 一 一 
1+b thr + hr’ 1+ or: HO — 电工 


25. 求 六 Tz) 定 在 x 二 0 处 前 (2.1) 阶 粕 捕 带 近 Rs, (x). 


解 ” 设 Rtr) = 和 ftir) = 二 ef 的 Taylor 公式 
] 
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由 一 一 cz 名 ,得 站 一 一 


知 二 1 二 lr ee 了 “ 


四 | 一 


a 二 1 二 


a 二 二 中 机 十 cc 一 二 


rr 


54 数值 分 析 全 析 精 解 


第 4 章 ”数值 积分 与 数值 微分 


1.1 重点 、 难 点 全 析 


1.1.1 数值 求 积 公式 
1. 裁 值 求 积 的 基本 思想 
积分 TCP) 一 | emar 的 数值 求 积 公 式 可 表示 为 


而 
1) = | fdra AFCr) C4. 1) 
“ = 
数值 求 积 公式 54. 1) 的 余 项 为 
ROP = IN — YA) (4. 2) 
t=" 


式 中 民 : 称 为 求 积 节点 ,A 称 为 求 积 系 数 . 亦 称 伴随 节点 x 的 权 , 六 翁 与 节点 
的 选取 有 关 ,而 不 依赖 于 被 积 另 数 F(x) 的 具体 形式 ， 

2， 播 值 型 求 各 公式 

当 Fr) 用 Lagrange 插值 多 项 式 上 ,tz) 近似 时 ,构造 出 的 求 积 公 式 


1 = DAF) C4. 3) 
称 为 是 插值 型 的 . 式 中 求 积 系数 4 通过 插值 基 消 数 避 tr? 积分 得 出 : 即 


Tri) tx) 


叶 T 二 0 1 4. 4} 


让 让 
六 = | (ed | 


OT 3 (re) 
式 中 心计 C) 一 《一 你) 一 
对 于 搬 值 型 求 积 公式 (4.3) ,其 余 项 
bk TH] > 
RADI-D | dr (4. 5) 


J nt 1) 
式 中 上 与 变量 > 有 关 . 
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3 求 和 公式 的 称 定 性 

定理 4.1 车 求 积 公式 半 , 1 中 系数 生 汪 0 二 1、 87) 则 此 求 积 公 
式 是 稳定 的 . 

4， 求 各 公式 的 代数 精度 

求 积 公式 (4,1) 苦 对 于 次 数 不 超 过 m 次 的 多 项 式 均 能 准确 成 立 , 但 对 于 
好 十 次 的 多 项 式 就 不 能 准确 成 立 , 则 称 该 求 积 公式 具有 严 次 代数 精度 . 

定理 4.2 形 如 式 (4.3) 的 求 积 公式 至 少 有 次 代数 精度 的 充 要 条 件 是 它 
是 插值 型 的 ， 


1.1.2 生 顿 - 柯 特 斯 公式 


设 将 积分 区 间 7a .条 划 为 等 分 , 步 长 入 一 < 二 < 选取 等 距 节 点 +， 一 a 十 
地 张 sl ay 则 可 构造 出 牛顿 - 柯 特 斯 (Newton-Cotes) 求 积 公式 
一 Ch a) Yee pen) C4. 6) 


式 中 CC” 称 为 柯 特 斯 系数 ， 
nn 一 时 . 求 积 公式 为 梯形 公式 


[fn 十 Fb)] (04. 7) 


[fd 
一 2 时 , 求 积 会 式 为 辛普森 公式 
pp-—u) ， 
本 | pe 和 
了 二 时, 求 积 公式 为 柯 特 斯 公式 
| Flrdr ss 7 fr,) Aa 12Ftry 十 3270r) 7A) 


(过 2) 


[fendr 人 )+ AD (4. 8) 


本, 99) 


式 中 ,ri 二 gg 二 万， A 


当 呈 六 如 时 , 柯 特 斯 系数 和 7” 诈 现 负 值 , 此 时 计算 不 稳定 ,故地 袜 日 时 的 牛 
顿 - 构 特 斯 公式 是 不 使 用 的 ， 

2.， 代数 精度 

定理 4.3 当 n 为 闸 数 时 .牛顿 - 柯 特 斯 公式 (4.5) 至 少 有 nn 次 代数 精度 ; 当 


1 
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# 为 偶数 时 . 则 至 少 有 n+ 1 次 代数 精度 ， 
3， 见 种 低 阶 永和 公式 的 余 项 


梯形 公式 : 
KR: =— a, PE Ca. bi 
辛普森 公式 ， 
p—arts 
Rs =~ 180 { 7 (a, 人) 
柯 特 斯 公式 ， 
42 一 2 
R (Fe) poem ye te 


1.1.3 和 做 化 求 积 公式 


{4. 107 


{4. 127 


把 积分 区 间 [Lutb: 分 成 4 等 分 ,在 每 个 子 区 间 上 应用 低 喻 求 积 公式 的 方法 


称 为 复 化 求 积 方 法 . 
人 
1 T= $ [r+ + 0] 
其 余 项 为 


RN =~— TR, jE ab) 


th 
~ 12 
复 化 辛普森 公式 (一 人 -2 ) 


Ts S$, = ft4D fn. 2 fe + 0) ] 
其 余 项 为 
Rf) = 一 和) (PD Ea 
1.1.4 龙 贝 格 求 积 公 式 
龙 贝 格 求 积 算法 的 计算 步骤 可 用 下 表 表 示 , 其 详细 过 程 见 教材 - 


4.137 


Ka. 1 2 


{£4.15) 


4. 16) 


be 
~ 
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区 间 等 分 数 梯形 序列 辛 浦 生 序 列 | 柯 特 斯 序列 龙 贝 格 序列 
如 2 Tut St- ut Rr1 


0 2:=1 TL | 


1.1.5 高 斯 求 积 公式 


1 一般 理 论 
定 尽 二 1 [a, 5 上 上 带 权 ot.r} 的 积分 了 一 上】 reeeceyar, 它 的 求 积 公式 为 


， 。 
全 roomd a DY ALF) (4.17) 
= 


J 


若 求 积 公式 5d4. 17) 具有 2n 十 1 次 代数 精度 ,网 称 其 节点 zz (二 0 1. 下 ,1) 为 
高 斯 点 ,相应 的 公式 (4.17) 称 为 高 斯 求 积 公 式 . 
定理 条 节 播 值 型 求 积 公式 [17) 的 节点 6 受训 之 旋 是 
高 斯 点 的 充 要 本 件 是 ,以 这 些 节 点 为 零点 的 密 项 式 
WT = Er Tr 
委任 何 次 数 不 超 过 = 的 多 项 式 忆 671 带 权 pt.x) 正 变 . 即 
{powm np dr = 0 C4.18) 


根据 定理 二 4 在 [es 让 上 带 朴 ptx) 的 正 交 孝顺 式 的 零点 即 为 高 斯 求 税 公 
式 的 节点 ,其 系数 


站 2 . 
六 二 | OA) dy C4. 19) 
{TT Te Sr be) 


其 余 项 为 


rl rh 
R(t) 一 i Cryotridr, WE ta. db) (4. 20 
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定理 4.5 高 斯 求 积 公式 (4, 17) 的 求 积 系数 训导 二 0， 1 al 全 是 目 
的 , 压 高 斯 栈 积 公式 荐 稳定 的 . 
定理 4.6 若 六 x) EE Cia 万 1, 则 高 斯 求 积 公式 (4.17) 是 收效 的 , 即 
tim DAF Cs) = | ffirpC) dr (4. 21) 
常用 的 沿 斯 型 求 积 会 式 有 :高 斯 - 勒 让 德 求 积 公式 ,高 斯 - 切 比 雪夫 求 积 公 
式 ,高 斯 - 拉 盖 尔 求 积 会 式 及 高 斯 - 埃 和 尔 米 特 求 积 公式 等 .可 参考 教材 . 


].1.6 数值 微分 


数值 求 导 数 常 用 泰利 展开 法 , 恤 值 孙 数 法 .数值 积分 法 和 三 次 样 条 法 等 ,也 
可 使 用 理 吾 森 外 推 法 . 考 虚 到 会 八 误差 的 因素 ,数值 微分 公式 中 的 步 长 一 般 不 
宜 取 得 过 小 ,具体 求 导 公式 可 熏 考 教材 . 


1.2 习题 全 解 


1]， 确定 下 列 求 积 公 式 中 的 待定 季 数 ,使 其 代数 精确 度 尽 是 商 ,并 指明 所 机 
造 出 的 求 积 公 式 所 具有 的 代 米 精度 ， 


| re A A + A fs 


oh 
2) fndr se A fA AO TAA): 


rl 了 一 " 本 
(rod 到 了 十 2 十 27， 
二 


Gy | Fondr a 1] Fh Pt0) 一 产 (6 


解 0 将 并 xz) = 1 x 分别 代 入 公式 两 端 并 令 其 左右 相等 ,得 
AAA := 2 
| ,Ra 一 站 


LA 1 十 大 六 邓 并 
解 得 入， 一 A = 全， 和 A 一 加 ,所 求 公式 至 少 具有 2 次 代数 精确 度 ， 又 由 于 


让 
本 
| .= dr 二 本 (一 科 ! 十 司 * 所 
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站 
| rdr 尖 在 (一 册 十 在 。 外 
a 中 电 


fl)de fh) 7(0) | 生产 5) 具有 3 次 代数 精确 度 - 


{20 仿 I = ] Tr ,信人 公式 两 端 令 其 左右 相等 ,得 
4 一 扫 
一 


蔓 | 


， ,9 
和 一 二 A 十 A 一 了 


解 得 A = A Eh A C=- 


令 Ar = ,得 
Ps & 
=| rr 
| " 
令 Fr) = x! ,得 
16 


人 
Sh 全 = dr ALt 而) 十 页: = 可 各 


页 求 积 公 式 具 有 3 次 代数 精度 ， 
(3) 当 .Ar) -= 1 时 , 易 知 有 


上 
三 Frydr = 工 [Lr( 一 有 十 270c 上 37 


令 求 积 公 式 对 Flr) 二 x. 淮 确 成 立 . 即 
i—] Far | dr 


十 2ri 十 3 二 2 


则 可 解 得 
[: —0.2899 1 = 0.689 9 
irs 0. 326 [rs =—0.1266 
将 Fr = x 代入 已 确定 的 求 积 公式 , 则 


| eez Zz 亏 [7C 一 1 2ftry -ht 3H.rs)] 
破 求 积 公式 具有 2 次 代数 精确 度 ,所 尿 节 点 为 z, 一 一 小 289 9, .ro 一 .526 6 或 
| 一 0589 8，r， =— 0. 126 6. 

{4} 求 积 公式 只 含有 一 个 待定 参数 mw 当 产 r) 一 1， 时 ,有 


网 


“ 下 a A . 
lde= Ti+ lt .0-0 | ldr= 1+ 
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| xd 一 生 [0+ 问 十 咕 :0] 一 1] 
故 令 FIz) 一 天 时 求 积 公式 准 夏 成立, 即 
[ar 一 全 [0 十 有 + 大 :[2x 0 一 2p] 


= 
解 得 a 


将 大 7 = ,zt 代 人 人 上述 育 定 的 求 积 公 式 , 有 
由 2 
fear 二 全 [0 十 和 如] 十 等 [0 一 3 和 


看 2 
iar Es 到 fo 十 岂 志 十 生 10 一 463 


故 所 求 求 积 公 式 具 有 3 次 代数 精确 度 . 
2. 分 虽 用 梯形 公式 和 幸 普 森 公 式 计算 下 列 积分 . 


1 . 四 3 二 
| Idx A 8 | Ue DF, 0 
(3) | Var 4 《4 下 VA sinip dp, n = 6. 
解 ”1 复 化 梯形 公式 一 8a 一 0 一 17x) = 了 4= 吉 


T= 和 [ft fn) J] 0.111 #0 
复 化 学 普 琳 公 MR 


一 三 A pea) + 4 DD ps 2 fm tft] 111 57 
A 
1 


(tle) 
z 


(2) 一 10 二 0 一 1 fC) 
复 化 樟 形 公式 为 
Ty = f+ 2 fem) + 7] 一 1.391 48 
复 化 辛普森 公式 为 | 
Sm = [f(a) + fm te De + f(D)] = 1.454 71 


C3) N= 1h = = fr or 
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复 化 梯形 公式 为 
T, 一 [f(a) 十 2 fe 十 FY)] 二 17, 227 74 
复 化 辛普森 会 式 为 


3 4 
Si Lf HID fm) +t oD Fr) + /6)] = 17. 322 2 
k= Fe 


让 Er = vA sn 


4)n Ba Dp 


6 
复 化 梯形 公式 为 
T; = Ei 十 2 Fa + hy] 一 1.D35 62 
复 化 辛普森 公式 为 


Si = EO tA rms) 2 ft) + /CB)] = 1.085 77 
3. 直接 验证 各 特 斯 公式 [4.9) 具有 5 次 代数 精确 度 . 
证 阴 柯 特 斯 公式 为 
pa 
rendr = 0 
令 rT) 二 1, 则 


[Ff Cm) 32F CT) + 12F Cr ) + 32F (73) + 7 FC,)] 


让 
| 六 fr 一 出 一 应 


[7 for, 十 3270m + 2 fr + 32F(7 + ?Fr = ba 


令 Aer) = zx: 则 


nh hb 1 
| Flridr 一 | Xdx 一 上 全 
看 


TTF) + 2) + 12 Fr) 1 3270c) 二 ?Fr 一 去 必 — a) 


令 Frl = 二 x, 得 


[renee | dz 1 《8 一 全 3 
a 3 


-nn 


i + 32F¢r) — 12 fr) + 27 Cr) + 7 F711)] = 言 (5 四 


令 Fr) = 所 , 则 
| renaz 一 | xsdr ~ 3 Ch — a!) 


62 


歉 值 分 析 全 析 精 解 


误差 


两边 


pC— 
on 


全 Fir} 二 x, 则 


中 


QT Fer) + 32F Cr 12f ry + 92 Cr TH) = hat) 


人 5 1  - _ 
| reom 一 | :dr 一 六 —a1) 


2 BY Fm) TI2 fr +t 27) + TF rz = 于 他 —a) 


令 站 77 二. 则 


上 rdr = | rd par) 
fr 了 一 | 7 人 -~ 
7 fe ?十 327r) 12Ftr) 二 + 327trs) ?fr = 二 ay} 


他 六 io 二, 则 
在 一 
90 


因此 ,该 柯 特 斯 公式 具有 5 次 代数 精度 . 


4, 用 半 普 说 公式 家 积 分 | edz 并 售 计 误 基 ， 


ETF fm) fF BOF Cr HE I2f Cr) TF B27Crs) 4 ?Fr1)] 


8 
| rd 区 


解 3:= re "de 二 el] 一 0.632 33 


一 也 一 531 
2 2 】 的 


1 1 _ 
1 ~ a ~ = DD. DD 35, FE ,1? 
5, 推导 下 列 三 种 蝶 形 求 积 公 玄 ， 


nm 


| RiF)Y = - 


{1} 


bh 
fondr ss 十 一 ay Flay 十 LR:s 


J 


(2) | firydr = (ha) ftby.. Dba) 


/Ddr = Hof (se) re). 


解 (1)7 左 算 形 公式 ,将 天 六 在 “处 展开 ,得 
一 一 
在 [au, 户 ] 上 积分 ,得 
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[fondr = | radar rec-odz= 
(ph—a)ftay + | re 一 adz 
由 于 x 一 5 在 La; 如 上 不 变 号 , 故 由 积分 第 二 中 信 定 理 : 有 ?YE (ay 5) ,使 
[erar = Bf +f P| de 


故 有 | erar = Ba fn a FE Ca: 
《2) 右 纸 形 公 式 , 阿 (1) ,将 rr) 在 点 处 展开 并 积分 ,得 
hb 
| firydr = (Half tb) 一 二 广 (用 志 一 oj ， nCa, b) 
(3) 中 短 形 公式 ,将 f(x) 在 “于 处 展开 ,得 


ft) (一 < 人 人)+ 


rn (rt). 二 E Ca. 有 
两 边 积分 并 利用 积分 中 值 定 至 ,得 


[rnd = (2)6 ay 7 (oe)| (re 人 二] dr 十 
二 | refz- 上 二 2) qr = 


他 一 oj(2 二 全 二 二 产 全 (一 经 2) dr 局 


Bf (ee) 站 六 tp 一 沾 1 ， IE Ca, 


6, 车 用 复 化 实 形 公式 计算 积分 『 一 | eae, 问 区 闻 [oJ] 应 分 多 少 革 分 才 


能 使 截断 误 兰 不 超过 x 10-*? 著 改 用 复 化 说 普 森 公式 ,要 达到 同样 精度 .区 朵 


应 分 多 光 等 分 ? 
解 于 Fa 一 er 一 Pr 一 ea 一 1. 故 对 复 化 梯形 公式 ， 
要 求 
， 由 一 2 1 fl -1 : 
| RIi{f) i= |— |< Ta (地 ) eXI0', y€ {0,1 


6 和 数值 分 析 全 析 梧 解 


即 六 实 语 X107.n 守 212.85. 取 nw 一 213, 即 将 区 间 [0,1] 分 为 213 等 分 时 可 以 


满足 误差 要 求 . 
用 复 化 辛 普 琳 公式, 要求 
， | 让 一 af 页 和 Fa 1 1 
人 - Teo CE) 7 人 | < mr (3) ‘< 


1 

之 

即 由 执 全 x 10m 371. 取 2 一 上 即将 区 间 等 分 为 8 等 分 时 可 以 满足 误 养 
要 求 . 

7 如 果 (x) > 0, 证 明 用 梯形 公式 计算 积分 了 一 | f(x)dx 所 得 结 果 比 淮 


确 慎 了 大 ,并 说 明 其 几何 意 头 . 
证 明 ”由 梯形 公式 的 余 项 
(ay 
jpD, FE bash) 
知 ,车 产 ( 了 0 且 请 2 网 民 CR < 之 0, 从 而 


5 
Jf evar 一 T+RCP 二 工 


X10 FE (0,1 


RfY =— 


即 计 算 值 比 准确 值 太 . 
其 几何 意义 为 ,六 tx) > 0, 融 六 rz) 为 下 凸 函 数 , 梯 形 面积 大 于 曲 边 梯形 
面积 . 如 图 4. 1 所 示 ， 


| y=) 


8. 用 龙 贝 烙 杂 积分 方法 计算 直列 积分 .使 误差 不 超过 10. 


Sr 9 
4) | erdes (| zsinrdr 09) | z VITZd. 
|， ， 
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解 ”(1) 计算 如 下 表 所 示 ，; 


Th? 
0 D771 了 43 3 
I NM. FT28 059 9 0.713 51l2 1 
2 D.716 9B2 8 0. 713 287 0 D. 713 272 0 
3 D.714 200 2 OQ. 713 272 6 0,.713 271 7 D713 271 了 


因此 ,Ts>= 0.713 271 7. 
《2) 计算 如 下 表 所 示 : 


了 


人 


3.451 313 x 10°5 


38.628 283 x 10 一 二. 446 923 X 10"21 


因此 ,T xs 一 4.446 923 X 103 a 0. 
{3} 计算 如 下 表 所 示 : 


六 | Te Ty Tt Ti Tit! 


0 4. 230 249 5 


1 .171 $69 9| 10.151 743 4 
2 | 10. 443 796 8 | 10. 201 272 831 10. 204 574 4 
3 ] 10, 266 367 2 | 10. 207 224 0 | 10. 207 fn20 了 | 10. 207 669 1 


4 | 10. 222 270 2 | 10, 207 57F 2 | 10. 207 504 3 | 10. 207 593 9 § 10. 207 593 € 


5 | 10.211 260 7 | 10. 207 590 9 | 10. 207 592 2 | ]0. 207 592 2 | 10. 207 592 2 | 310, 207 592 2 


因此 .IT 10. 207 592 2. 
9 用 及 二 2.3 的 高 斯 - 勒 让 德 公式 计算 积分 | ersinzdx. 
] 
解 ” 居 为 +E [3j; 令 1 二 一 2, 则 :EE [一 1,1 
故 [ e" sinrtdr 一 人 er sinttt 2)dr 


6 数值 分 析 全 析 精 解 


3 _ _ 
| sinrdx 2 有 555 555 6 XX CFE 0.774 396 7 了 7) 工 Fn. 774 596 7) 1 


0D. BBS BBS 9 x FUNDY = 10.948 4 
当 n= 3 可， 


El 


je sinrdr 2 1. 347 B54 8 XL FAC— 0. 861 136 3) + fi0. B61 136 371 
] 


0.652 145 2 xX LF{—0. 339 981 0) 十 FAO, 339 981 0 = 
]D. 950 14# 
10, 地 球 卫 里 轨道 是 一 个 椭 国 , 补 轩 周 长 的 计算 公式 是 


3 一 a fl (£) sing dy 


这 里 “是 精 国 的 长 半 轴 以 是 地 球 中 心 与 加 道中 心 ( 术 轩 中 心 ) 的 距离 , 记 刻 为 近 
地 点 距离 .万 为 迁 地 点 距离 . 尺 二 537] km 为 地 球 半 径 ; 则 


2R7+7 TT 由 _ Hh 
站 一 里 人 5 


4 4 


a 


载 国 第 一 颗 人 造 节 球 卫 时 近地点 距离 为 卢 一 439 Km 远地点 焉 离 为 万 一 
2 384 km 该 上 东 卫星 轨道 的 周 长 . 


解 a = LRH+AY = 7 782.5 


36 5 3 


[ 
一 


六 二 


a 


48 3 


en 


1.361 B46 2 


[3] 


361 646 3 .S64 了 


了 1.364 B46 3 


.364 BAG 3 


后 和 


网 ”7 一 人 |= 10” 必 广 X10%*, 故 积分 已 有 ?位 有 效 数字 , 取 
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了 一 1.564 646， S= dul = 48 707.44 km 
11. 证 明 竺 式 


J 1 

- 亚 T 开 
HSIn = 了 十 = 
nn 31p sln 


试 依据 nsin Cn 一 3,， 6, 12) 的 植 ,用 外 推算 法 未 的 近 拟 值 . 


证 明 车 ftn) =， nsin 至 ,此 函数 的 素 勒 展开 式 为 


着 记 TY 一 wain 于 we r 其 误差 为 Of ( 工 ) ) 


由 外 推 法 ,TS 一 计 04I 和 8 -了 ) ~ x 其 误差 为 O( (三 ). To = 


站 067& 一 TE) = ,其 误差 为 0( (三 ) )， 
据 上 述 公式 列表 计算 下 表 所 示 ， 


15 


.598 076 


3.000 ONn 3.133 975 


105 B29 3. 141 7035 3, E41 580 


鼓 x A 3. 141 58. 
12. 用 下 列 方法 计算 积分 | 冬 , 并 比较 结果 ， 
(1) 龙 贝 格 方 法 ; 
(2) 三 点 及 五 点 高 斯 公式 ; 
(34 特 积 分 区 间 相 擎 分 ,用 划 世 两 点 高 斯 公 臣 ， 


解 (1) 计算 见 下 宕 ， 


68 数值 分 析 全 析 精 解 


k Ti TH Ta! TT) Ty) 
| 上 | 

0 1.333 333 

1 1, 166 66? 1,111 111 

2 1.116 B67 1. 100 000 1.199 239 

3 1.103 211 1. 098 725 1.098 640 1.098 531 


4 1.099 ?68 1. 098 620 | 1.098 613 1 .098 813 1. 098 6l13 


联 1 二 1. 098 613, 
(2) 积分 区 间 -1, 3] ,来 用 高 斯 公式 时 , 需 先 变换 到 [一 1, 1] 上 . 喜 作 变换 


= Lp of 
YY Fa 十 折 } 十 信人 qayt 二 1 十 2 
| 
则 当 y EE [1 时 ,1 E€ [一 1， 1 dy — de, | gf dt 
[用 十 了 
三 点 高 斯 公式 


| 站 一 | ps 一 0.555 555 6 (F077 007 + FTO TT S987) 
D. S88 B83 9x pr 一 1.098 039 
五 点 高 斯 公式 
,由 -| 7 和 0 236 826 9 (0 J 3+ Te 17 6)+ 
D 478 628 "(ss M+ TT)+ 
0. 568 888 9 x zi = 1.098 609 


(3) 将 区 间 [1, 3] 四 等 分 ,在 每 个 小 区 隔 上 用 两 点 高 斯 公式 ,得 


1-5 t 
二 dy_ [| od 
i | 由 1 


1 1 
-1 A++- 4 

- TY 1 |= § 405 4 

non ns roel rv 0. 405 405 4 
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『 1 十 1 | 
nD. 3x 1 1 |=0.287 6712 
一 二 3.5++0.5xXx— 
[3.5+0.5x { 局) 霹 
»=| 宇 - ”osd 
2 永 -1 44.3 二 0. 3 
TCD 
0.5X 1 1 |= 0.2234 1405 
. 5 一 二 4.5+4+0.5Xx—— 。 
[5+0.5 ( 方 ) ; 万 
7 站 dh osd 
ay 1 5.5+ 人 O53 ~ 
1 1 
0.5Xx 1 --， ] i= 0.182 3204 
5,.5 二 0.5X 1 一 一 = .03 二 .5X - 
| ( 雄 ) 三 
艳 了 一 上 十 天 十 下 十 万 1.098 538 
13. 用 三 点 公式 和 积分 方 活 求 fr) 一 林产 二 在 z 一 1.0. 4.1] 和 1.2 处 
了 


的 尾数 值 ,并 估计 误差 、Fz) 的 值 由 下 表 给 击 : 


解 “由 带 余 项 的 三 点 求 导 公 起 可 知 


Ft) = 去 3Fte) tf fr)] + 人 


fr) 去 [ Flr) + fr)] FS) 


fr fe) 
志和 (xnyrz) ai 一 日 1，2 
取 上 表 中 心 一 1.0, zm 二 11, ro 一 1.2, 分 别 将 有 关 数 值 代 人 以 上 三 趟 , 即 可 
得 导数 近似 值 . 


由 于 六 人) | 和 max | F(z) |= max 
oT. ] 


村 


从 而 可 求 得 误差 上 限 与 导数 值 如 下 表 所 示 ，: 


4 | 
tr | 


4 _ 
2 一 0.75 


70 数值 分 析 全 析 辆 解 


1.2 


—D. Er 


.O02 5 


— .187 B28 了 


一 由 215 959 4 


数值 积分 法 . 设 gtx) = 了 (xz), 由 
fx) = Fr 十 | gen)dr 
了 


对 积分 采用 樟 形 公式 ,得 


fn) Fe 十 ptre) 十 os 一 


《sr 一 
tg 《全 ]， 巴 所 ZE Ter) 


令 上 万 二 0. 1, 得 


PT) PTL) A LF) 一 f(r )] 
| 2 , 
ELT) 二 各 (Te ) A 下 [Fa 一 Flr; 
同样 对 Ti 一 证 | gender 
el 


el 


有 fril? = Hr) +t yr) + px JJ? 1 一 


(Tol 一 3 


13 Pm), hE CFL Tl 


眠 而 有 
Tr) (rs) a LF) ~ Fee)] 


代 和 人 数值 , 解 方程 ,部 得 ptc)tk 二 0，]，2) 如 下 表 所 示 : 


l. 


2 


| — 0.187 


一 0D.213 959 4 | —0.187 828 7 


0, 000 828 7 


D. 001 O40 对 
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第 5 章 解 线性 方程 组 的 直接 方法 


1.1 重点 .难点 全 析 


1.1.1 高 斯 消去 法 


对 方程 组 4x 一 由 其 中 生 非 奇异 , 高 斯 消去 法 是 将 方程 组 约 化 为 一 个 于 三 
和 形 方程 组 ,然后 回 代 求解 .而 高 斯 - 若 当 消去 法 是 将 方程 组 的 系数 矩阵 约 化 
为 单位 阵 , 则 由 应 的 右 端 列 向 量 就 约 化 为 方 囊 组 的 解 .在 消 元 过 程 中 .为 避免 
主 元 素 相 对 其 他 元 素 过 小 .以 致 引起 舍 太 误差 过 太 的 情形 .可 将 向 列 中 绝对 值 
最 太 的 元 素 搞 到 主 元 素 位 置 上 , 即 为 列 主 元 素 消去 法 . 

利用 消去 法 计 求 扼 阵 的 道 . 行 列 式 及 鲜 杀 数 息 阵 相同 的 线性 方程 组 系 . 


1.1.2 和 抢 阵 三 角 分 解法 


高 斯 消 元 过 程 实现 了 汪 的 一 个 三 角 因 子 分 解 各 = LU, 其 中 卫 为 单位 下 二 角 
阵 ,5 为 上 三 角 阵 . 此 分 解 称 为 4 的 杜 利 特 尔 CDPootirtle) 分 解 . 

定理 $5.1 各 有 唯一 的 Doolittle 分 解 的 充 要 条 件 是 4 的 前 一 1 个 顺序 主子 
式 均 不 为 零 ， 

定理 5.2 若 4 对 称 , 电 4 的 所 有 上 顺序 主子 式 均 不 零 , 则 4 有 瞧 -的 分 和 解 式 

A = LDL' 

其 中 工 为 单位 下 三 角 阵 , 卫 为 对 角 阵 ， 

定理 5,3 若 和 4 对 称 正 定 , 则 存在 一 个 实 非 奇异 下 一 角 阵 工 .使 4 一 于 5. 
当 限 定 上 的 对 骨 元 为 正 时 ,这 种 分 解 是 唯一 的 . 

对 定理 5.3 的 分 解 忒 , 因 分 解 过 程 中 有 并 方 运算 , 故 称 平方 根 分 解 或 
Chelesky 分 解 . 

点 分解 后 , 解 方程 组 Ax = 一 5 等 价 于 解 两 个 三 角形 方 穆 组 
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[y= 


Ly==b 
[一 政 当 4 对 称 正定 时 ,| 
=y Lr-=y 


Ux = 


11.3 向 量 和 系 阵 的 范 数 

1. 定义 ”在 R" 上 的 实 值 函数 ‖ .|i 称 为 向 量 范 数 , 如 果 对 R" 中 的 任意 向 
景 x 和 yy, 它 满足 ; 

C1) 非 负 性 :xl 宇 0, 量 但 当 x = 二 0 时 |x| = 0: 

C2) 齐 次 性 ;对 任 一 数 kER, 有 | Rx = 二 |k&i xl 

(5) 成 立 三 角 不 等 式 : 上 x+y 达 ixl| + 1iyl. 

2. 常用 的 向 景 范 数 是 p - 范 数 . 


co - 范 数 (最 大 范 数 )， | 

1 - 范 数 ， zl = Ta 

2 - 范 数 ， Ix: = (D7) 

pp - 范 数 ， lx = (Dn). 1 pete 

3. 定义。 定义 在 R* 上 的 实 值 两 数 | 。| 称 为 矩阵 范 数 ,如 果 对 于 R”" 
中 的 任意 矩阵 4 和 B, 它 满足 : 

(7 141 尖 0, 且 仅 当 4 一 0 时 1 41 = 60， 

C2) 对 任 一 数 & ER. 有 1 态 目 二 | 站 | 入， 


C3) | 起 十 BB 所 上 | 二 B1; 
4) ABI < 过 14i181， 


4- 常用 的 矩阵 范 数 

co - 范 数 ( 行 范 教 ) 141- = me ,| 
1 - 范 数 ( 列 范 数 ): 1Al = max D) | as | 
2 - 范 数 : 41s -Ar 
FE - 范 数 ， IAlr = (Ba) 


a jl 


其 中 hnw(4 7 4) 表示 半 正 定 矩阵 4 4 的 最 太 特 征 值 . 矩阵 的 前 宇 种 范 数 分 别 与 


第 5 意 ”和 解 上 线性 方程 组 的 直接 方法 73 


向 量 的 = - 范 数 ,1 - 范 数 以 及 2 - 范 数 相 容 . 
1.1.4 误差 分 析 


设 1。 | 为 非 奇 民 算 阵 4 的 菜 种 算 子 范 数 ， 称 数 Condt4)。 一 
141.14a721u= 一 12 或 < 天 ) 为 纸 阵 各 的 条 性 数 . 当 和 的 条 件数 Cond( 间 )， 
准时 , 称 方 程 组 4x 二 是 "病态 "的 ,否则 称 为 " 良 才 ”的 ." 病 柱 ” 方 程 组 很 难 
用 常规 方法 求 得 比较 准确 的 解 ,但 对 其 中 部 分 方程 组 ,可 通过 余 量 校正 选 代 求 
解 方程 组 . 


1,1.5 上 矩阵 的 正 交 三 角 化 及 应 用 


(1) 设 疝 量 w ER 有 旦 ww 一 1, 则 称 矩阵 肝 (w) = 了 一 2ww7 为 初等 反射 阵 
(或 称 豪 斯 杜 尔 德 (Householder) 变换 ),， 它 在 计算 上 的 作用 是 实现 约 化 身 量 与 
钴 阵 , 如 对 向 量 x 关 8 则 可 选择 日 使 Hx 一 士 x 二 se. 

2) 平面 旋转 变换 ,又 称 (iivens 变换 ,可 通过 选择 平面 旋转 矩阵, 特 向 量 中 
指定 元 素 变 为 零 . 

‘3) 矩阵 的 QR 分 解 , 使 用 多 次 初等 反射 变换 或 平面 旋转 变换 ,可 实现 4 的 
分 解 丰 地 .其 中 心 为 正 交 算 阵 ,及 为 上 三 角 阵 . 当 限 岩 尽 对 角 元 为 正 时 .此 分 
解 唯 一 . 

44) 利用 4 的 正 交 约 化 可 求解 超 定 方程 组 . 


1.2 习题 全 解 


三 1 了 
1. 设 点 是 对 称 阵 且 ai| 关 0, 经 过 高 斯 消去 法 一 步 后 , 生 多 化 为 | 。 上 
证 明 A 是 对 称 和 矩阵 . 
证 明 ”由 消 元 公式 及 4 的 对 称 性 得 


。 a Qj 。 
2 1 一 1 = - 
dy = a a = di a j= 


冬 t1 Ql 
故 4 是 对 称 卸 阵 . 
2. 向 各 一 fa 是 对 称 正 室 上 元 和 阵 ,经 和 过 高 斯 消去 法 一 而后 点 鲍 化 为 


1 ®] 
[ 0 4: 


了 4 数值 分 析 全 析 糖 解 


其 中 4 一 Ca) 证明。 
(1) 入 的 对 外 元 妹 a 人 一 1，2， Pi 
{2) 4 是 对 称 正 定 姓 阵 . 
证 明 (1) 办 各 对 称 正 定 , 故 
Bi = (Ae;, 6) = eiAe: > HO i= 11,2, ,7 


其 中 8 = 0, …, 90， 1, 0 0)7 为 第 i 个 单位 向 量 . 
{2) 由 起 的 对 称 性 及 消 元 公式 得 


137 由 Hil 


妾 5 ty a i a ap ; i+J 二 2 nn 
故 有 4; 也 对 称 ， 
1 er 
及 | 上 = 工 :如 
0 六; . 
1 
2 1 
如 
其 中 也 一 ， 11 
Rn ] 
Ql 


显然 工 ; 非 奇 异 ,Yx 天 0, 有 
Ex 天 QiCx LALIX) = (Lix, arTx) 0 《由 和 4 的 正定 性 ) 
故 LtAL,” 为 正定 矩阵 . 


11 口 
又 上 ALT = 人 ， | mn > 0, 故 4; 为 止 定 答 阵 . 
3. 设 亏 为 指 禄 为 下 的 初等 下 三 角 阵 ( 除 第 点 多 对 衣 元 以 下 元 素 外 ,L， 和 单 
在 阵 于 相同) ,好 
下 
1 
1 


LL; 一 . 
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求证 当 记 了 二 时 全 一 瑟 上 L441, 也 是 一 个 指标 为 下 的 初等 下 三 角 阵 ,其 中 局 为 


初等 置 措 阵 . 
证 有 明 因为 
要 
| 
1 = 上 
Pa I 
7 上 1j 


上 且 4 为 初等 冒 搞 阵 , 当 ?> 上 时 ,将 的 第 ; 行 和 第 j 行 互 换 
门 . 一 


! 
Ppt].# 
1.L = : 0 -= 1 i 
: : 日 1 
1 1 5 + 
L 1 


将 上; 的 第 i 列 与 第 j 列 互 搞 
Li = = 
所 以 六 为 指标 为 & 的 初等 下 三 角 阵 。 
4， 斌 推导 些 阵 上 问 的 Crout 分 解 各 一 工 D 的 计算 公式 .其 中 工 为 下 三 角 痊 阵 ， 
可 为 单位 上 三 角 瞪 阵 . 
解 年 阵 有 4 的 Crout 分解 生 二 LU 


1 [3 "1 ma 1 .203 a 
[5 la: fa 1 ae 


ba, Hr a | 1 | 
由 和 手 阵 嫌 法 知 a = i 2. 
故 上 


又 由 au 一 020i 故 2 一 2 了 一 2 3 a 
11 
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假说 工 的 前 上 一 1 列 和 的 前 上 一 工行 已 知 . 出 


上 上 -i 


Qn Daun = Dts FH i= 二 lm 


故 fi 二 a 
同 理 有 


内 .| 
au 一 ts =: Dus haus, j= 一] 十 2 


Eo 
| 

2 一 > teu 
rl 


故 ui 一 f ， 了 一 点 十 T 下 二 号 
3 


综合 上 面 推 导 结 果 得 Crout 分 解 公式 为 


ty dant i= 1 2,,n 


Hl) 一 了， 2 

1 nn J 
| 

A = a Duas, 了 一 下 点 十 1， 
rl 


A—] 
2 

Up dn 
r= 


Hy 


5. 设 [x 旦 ,其 中 在 为 三 前 蝶 阵 ， 


(1) 就 U 为 上 及 下 三 角 丐 阵 推 导 一 般 的 求解 公式 ,并 写 出 算法 ， 


(2) 计算 解 三 角形 方程 组 [Ux :od 的 来 除法 次 数 ; 
(3) 设 下 为 非 奇 异 阵 , 斌 推导 求 品 ， 的 计算 公 武 . 
解 1 设 U 为 上 三 角 阵 


HH ln "1 i 
Ha Ha Ls i 
Hm i -上 


了 一 此 十 下 天 十 2 2 克 一 2 3，，- 
Ea 
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二 一 > Hy ke 
:和 , 
= 1 一 一] 一 
Hn 


Xr 


当 口 为 下 三 角 阵 时 ， 
13 TI | a 

Hl [a 2 加 

了 


Th | 


2 te 


ol 
:一 ] 
A 一 > 上 
部 of 1—l ， : 
得 zl ,x 二 一 
| 地 语 


《2) 除法 次 数 为 #. 先 法 次 数 为 
1 二 2 十 TD 一 


故 总 的 乘除 法 次 数 为 十 全 一 一 于 


43) 设 忆 为 上 三 记 王 二 S 则 SS 也 是 上 三 角 阵 ,由 


[ | "Kis 0] 
“fo 1 
Sn 
得 + 2, ne 
3 可 法 证 时 
员 mn in 1 2, 
当 V 为 下 王 角 阵 时 ,由 
Hl B11 了 
[ul ee Nal Ns2 1 
| Ha Sn Sm 
得 wss 1, i= 1.,2. 亚 
-1 
到。 
S 一 一 2 j= 
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8. 证明 :(1) 如 到 并 是 对 称 自 定 阵 , 别 态 '， 也 是 对 称 正 定 隆 ; 


(2) 和 如果 间 是 对 称 正定 阵 , 则 各 可 唯一 地 写成 各 一 1 .其 中 虐 是 县 有 还 对 


前 元 的 下 三 前 阵 - 


证 明 ” (1) 因为 坟 是 对 称 正定 阵 . 故 它 的 特征 值 1, 金太 于 0,&' 的 特征 亿 


hi 1 蕊 全 大 于 巡 


1 
故 是 对 称 正定 短 阵 . 


(2) 因为 生 对 称 正 定 , 故 和 的 所 有 硕 序 主子 式 均 不 为 零 . 从 而 4 有 了 唯一 的 


Doolittie 分 解 太 二 LU, 又 


1 名 En 
TL Wilt :1 
| Uyy Ha 
U = ] a | = PU, 
Hm : 
1 
其 中 P 为 对 角 阵 ,UU, 为 单位 上 三 角 阵 ,于 是 
= LD = LDU, 
于 站 二 A7 二? 
由 分 解 的 唯 -- 性 即 得 
Ur =L 
从 而 有 是 一 下 DT 
又 由 4 的 对 称 正定 性 知 
d= Di, ad = >0. 一 了 
:—] 
ra, 
dd 
邦 DD = 一 
可。 
| va va 
Vd a 一 DD 


第 5 章 解 线 性 方程 组 的 直 授 方法 79 


训 
A EDL = LDTDIET 一 《ED ) cp = 了 7 
即 矩阵 及 可 了 唯一 地 分 解 为 4 = [LT ,其 中 1 一 了 二 为 对 角 元 为 正 的 下 三 角形 
抵 阵 . 
?用 高 斯 - 若 当 方法 求 鼻 的 逆 钴 阵 ， 其 中 
2 1 一 3 一 1 


风 3 1 0 
1 2 4 一? 
1 .051 5 
F213 -1:1 9000 
和 DLL 0 70901900 

-112 4 一 2i0010 
Lo-l 5i0001 
P 5 一 1: 5:0 0 0 3 
01 3 —8i0 1 一 3 

2 3 1 o1 中 
0 1 -1 -100 一 2 
一 ] 5 ou D 1 

1 3 -8i0 10 一 3 
0 0 ~3 1945 —2 7 
no -4 ~3il 10 1 
- 4 ， 2 
100 -了 :0 于 -也 -| 
010 io -1 1 t 
5 
lb 00 -等 1 计 -村 -等 | 


80 数 慢 分 析 全 析 精 解 
广 ; 2 了 
ee 
so 
90010- 栈 订 一 霄 一 
| 刘 二 站 让 | 
故 
| 车 站 一 种 一 站 
33 6 4 13 
| 
3 3 _3 _8 
条 7 
.31 4 5 
| 1 17J 
r 0.047 0589 0,588 2353 一 9.2705882 ~0,941 1765 
0.388 235 3 ~—0.352 9412 0482 3529 0.7647059 
—0.223 5294 0.2941176 ~ 0.035 2941 —0.470 588 2 
[0.035 294 1 —0.0588235 0D.0470589 0,2941176 
8. 用 追赶 法 解 三 对 角 方 程 组 Ar 一 日 ,其 中 
2 1 0 0 0 
-1 2 -1 0 0 0 
4-|1 01 >-1 川 sl 
0 6-1 2 一 1 0 
0 0 0—1 ? 站 
解 设 有 分 解 
2 -1 0 0 0 a 0 0 0 0 
-1 2-1 0 0 1 ww 0 0 0 
0 -1 221 oo -1 % 0 0 
0 0~-1 2 -1 0 0 -1 wm 0 
0 0 0~-1 2 on 0 0 1g 
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9 1 Ro 0 
9 0 ] 0 
9 0 0 1 
Dn ,| 
由 公式 
自 一 mm 一 的 商 
ee i = 2 3 4.5 


c= oh, i 2 3 
其 中 性, ae 分别 星系 数 征 阵 的 主 对 角 线 元 素 及 其 下 边 和 上 边 的 次 对 角 线 元 
未. 故 有 


9 ,二 2 二 = 了 了， ,= 
人 1 + 安 : 2 在 3 了 rl T : 3 
一 -2 2 一 二 
癌 罗 、* 总 村 4 1 再 5 
由 
P 也 姜 0 站 妆 
3 
一 】 六 0 D a 4 ri 
: 0 
4 了 
全 一 1 一 0 0 
守 3 | 一 | 
0 0 一 1 地 0 | 
yd Lo 
在 
[oo 0 0 1 于 | 
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1 站 [0 总 
0 1 -全 0 0 
3 

_ 2 
0 0 ] 3 
a 0 0 1 一 生 
上 2 


解 得 x z 一 本 ， 一 二， 于 
1 


9. 用 政 进 的 平方 报 法 解 方 程 组 


T 2 —1 1 [Ti] 
一 】 一 2 3 322 [一 
] 3 1 lz; 
解 设 
| 2 一 1 1 1 了 Ta 
| 1 一 2 = |b 1 
i 1] 3 1: a1 14 L 
由 矩阵 乘法 得 
_ ll 
d! 一 上 at Eh ta 
一 互 一 了 
d; 一 2! dz 5 
2 
已 


clon 
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| 
ce 
| 
1 
上 
| 
cs| 一 


1 1 2 J 2 
1 | ra ! 1 
| 一 2 7 | i 
得 1 7 二 2 所 
?| Le 27| 国 23 
1 了] | | 3 
» 
访 x 二 二 9.355 555 6 rs 一 序 = 0.777 ?7? B.C= 总 :Titit 
即 x = 1, 11 111 1.D,777 777 RB,2,555 555 6)T 


]U, 下 述 手 阵 能 否 分 解 为 LET 其 中 于 为 单位 下 三 角 阵 , 开 为 上 三 南 几 1 车 能 
分 解 , 那 入 分解 是 否 唯 一? 


| 2 3 1 11 1 2 61 
A= |2 4 | B= | 2 | C=- |2 3 15| 
i 6 7 3 3 1 6 15 46] 


解 ” 坟 中 加 ;二 人 , 放 不 能 分 解 .但 dettA&) = 一 10 尖 D, 故 基 将 和 中 第 一 行 
与 第 三 行 交 拘 ,出 可 以 分 解 , 且 分 解 上 唯 一 . 
县 中 ,A 二 A 二 9 得 它 仍 果 以 分 解 为 


1 i 3 1 
加 = : 1 . 轩 —1 
ts nD 0 zj 


其 中 心 为 一 任意 常数 . 旦 五 奇异 , 故 分 解 不 唯一 ， 
对 Ca 关 D = 1 2 3 坝 妇 可 分 解 旦 分 解 唯 一 . 


Ei 


站 ,有 .8 了 
11. 设 各 一 | 】 。。 ji 计算 妨 的 行 范 数 . 列 范 数 ,2- 范 雪 以 及 下 - 范 热 ， 
.1 人 J 
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上 上 a1. = ms le ! 


Ti 


:5 
上 六， 二 aX >， | a, | 一 中 8 
J 


| al: = PE 
. 0.6 0.13 0.8 0.5 0.37 0.33 
4- [os oo] [or od- [os oa 
0.3 0.3 1 0.3 


.33 0.34 


Xm (和 一 0. 85 3 和 4 


故 让 各 一 Vi A 二 站.827 85 
12. 求证 :17 x. 芝 1 站 zi 委 王 | zi; 
1 ，， 
lAlre Ali: 41r. 
证 有 明 (1} 由 定义 知 
lx max lt i ln l= 


EE nx :: 


即 jx x nx 
{2) 由 范 数 定义 ,有 
放 太 1 二 A 
Da 1 tt Da = DD = jiAl 
?一 1 i=1 7=1 i=1 
| 起 :| := 一 cr > 这 
二 [LA CATA) 十 (AT 六 十 人 CaT4)] 一 可 a | 
故 


二 1 A A 
13. 证 PE R"" 且 非 奇异 ,反讽 x 为 R* 上 一 向 量 范 数 , 定 义 上 x-|。 
Pr | , 试 证 明 | xl ,是 民 " 上 向 重 的 一 种 范 数 ， 

证 明 (1) 因 P 非 奇异 . 故 对 性 意 的 x 关 0, 有 Px 关 鲁 帮 | x 
| Px |‖ 实 0;, 当 且 仅 当 x = 二 0 时 ,有 |x| ,= :Px = 成 立 ， 

(2) 对 性 意 a ER'., 有 


i 
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Torl,= iPrl| =ile| Pri 一 ae， xd, 
《3) x+ ys = Pexiyyil = PrPl+lPpyi = 
Hxil ,+t yl; 
故 1 站 zl 是 R" 上 向 量 的 一 种 范 数 . 
14， 设 各 所 Rom” 为 对 称 正定 姑 阵 ,定义 |x| 一 (Ax x) 了 , 斌 证 明 ixla 
为 R" 上 的 一 种 向 重 范 数 . 
证 明 (1) 因 生 正定 对 称 . 故 当 兰 一 和 时 ,xs =D, 而 当 半 短小 时 ， 
下 ={r Ar)t 0. 
C2) 对 任何 实数 a, 有 
or Ala 一 ai vxAr =|e| xl 


| ax i = 和 


(3) 因 点 正定, 故 和 可 崔 一 地 写成 上 4 一世 7 : 则 
jx = CA = x = LT Lx = HETx|, 
鼓 对 任意 向 量 x 和 和 3, 总 有 
Txt pli = itr+ryl: = |LTx+Liyil, < 


让 二 rz 十 Lyl: Co dxlat+ | rl, 
故 x。 一 (xTAx)3 是 -种 向 量 范 数 . 
15, 设 | 和 让， 141, 为 有 ”上 任意 两 种 矩阵 算 耶 范 载 , 证 明 冯 在 常数 
CC 守 0, 司 对 一 若 六 扎 R”"" 满足 
Ci 141 过 本 Cs i. 
证 了 明 因为 


由 向 量 范 数 的 等 价 性 知 ,存在 常数 cf ,cz >> 0, 使 对 任意 x, 有 


Car 之 让 area 
他 
了 了 4， | Axi 
A) i ， _- 1 i 了 中 
故 cixi Trl < er TT 
er _ Cr i 
令 一 二 .六 二 GC. 网 有 
Cu [时 | 
上 ax | ar， -， |Ax| 
， :eC, 
Tx Wx ~ xl 


8 


数值 分 析 全 析 籍 解 
、 | Ax|, | Ax |, 、 | ar | ， 
和 
即 中 1， 
16. 说 入 为 非 村 和 异 算 阵 : 求 证， 
] ,Ar | 
ATi 9 Ty 
和 -] [i 仿 辆 fr | 
证 明 4A]. 一 max 人 Xb A 
工区 和 0 "x | :oo 上 大 看 | Ay | 2 
70 TAF Te 
站: 
] ,| A | 
故 = 一 mnH 一 - 
ai en yi 


17, 下 阵 第 一 行 束 以 一 数 成 为 ” 


24 A 
< 
1 1 


证 明 当 4 一 二 字 时 ,Cond(A)。 有 最 小 值 ， 
证 明 设 4 关 0, 则 


3141，14! 亏 
让 六 “= 
2 ， | 壹 之 
3 
i 一 : 
A 
六 各: 一 
一 了 9 
El 
， 1 
|- 二 一 一 2 
发 1 Ti 
6 | | 十 3， | 之 
Condta) = IA -下 41、 一 


页 当 | 4 1 一 也 时 , 即 和 一 = 也 时 ,Condt4) .有 最 小 值 .此 时 


min CondtA).. 一 了 
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100 
: 设 问 二 i g8 oo | 计算 4 的 素 件 教 CondiA)，。 (t= 二 2, oo0). 
.i 98 99 | 
解 4 ! ] 上 A 一 199， | 和: 一 159 
L 99 一 100 
Condf 是 一 于 下 昌 让 一 39 #01 
_ 19 801 19 6027 
总站 一 | : 
19 802 19 405 ， 
、 _ Mn A . 
Condt 4 = |: 一 
故 end 一 和 FT 


19. 证 明 ; 若 点 是 正 交 阵 . 则 Cond(A&)， 一 1 
证 明 国正 交 , 页 "有 二 A47 一 开 4 一 上 4， 从而 有 
441: 一 vp A = vB =.1 
ailis = iA Tl: = volh4T) = vo = 1 
故 Condta): = 让 aa 1 1141: = 
20, 设 有 4, BER ,下 下， 为 R""” 上 大 阵 的 算 于 范 数 ,证 明 ; 
CondtABY < CondtAyCondtB) 
证 有 明 Cond(AB8) 一 | (CAB | 上 ABi | 上 aEis lilaAl lal = 
上 Cond(A}YCondtB) 
21. 证 4x 二 下, 其 中 肪 EER"* 为 非 奇 异 阵 , 证 明 ， 
(1) 六 "入 为 对 称 正 定 起 阵 ; 
{2) Cond{& > = [CondtA): ]. 
证 明 (1) 144)T 一 ATCAT)T 二 各 入, 孝 及 7 有 为 对 称 和 矩阵 . 
及 二 韭 冀 异 , 故 对 WxY 关 00. 有 Ax 天 0 从 而 有 
A 
即 474 为 对 称 正定 矩阵 . 
[2 Cond{tATA); 一 | CaTAY 11 4。 一 
VA [ (CAT AY 0 TT VAL CAT AITO4TA 一 
Vim [LAT A TD] 0 
VCATATTT VAba tA 
[ viara Ty] a 
a Cond( dy), 
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第 6 章 解 线性 方程 组 的 迭代 方法 


1,1 重点 .难点 全 析 


1.1.1 选 代 法 基本 概念 


对 线性 方程 组 Ax 一 上 其 中 由 一 (ai) ER 为 非 奇 异 息 阵 .可 变换 成 
X= Br 二 + (6.1) 
的 形式 ,从 而 构造 出 一 阶 定常 迭代 法 公式 
x = Bret ff ke 0, 1 2, ee 《6, 2) 
其 中 一 ”为 任 取 的 初始 向 量 ,B 为 常数 抵 阵 , 称 为 造 代 抑 阵 .了 为 常数 向 量 . 
车 limx* 存在 ( 记 为 x" ), 则 称 此 迁 代 法 收 俩 , 且 x" 就 是 方程 组 的 解 ,否则 
称 此 选 代 法 发 散 ， 
定理 & 区 ( 选 代 法 基本 定理 ) ”对 任意 初始 向 量 癌 " ,一 阶 定常 迭代 法 (6. 2) 
坡 合 的 充 要 条 件 是 掀 代 和 矩阵 B 的 谱 半 径 ptBY 一 |， 
定理 乱 .2{ 奢 代 法 收 争 的 充分 条 件 ) ”车 交代 拓 阵 8 的 某 种 算 子 范 数 
181 = 二 g 过 1, 则 一 阶 定常 闭 代 法 (6.27 对 任意 初始 向 量 x 都 收 敦 . 
定义 看 】 称 RiB) 一 一 Jnpt 有 ) 为 千代 法 10. 2) 的 渐 近 收 化 速度 ,简称 选 代 
读 收 化 速度 ， 


1.1.2 雅 可 比 选民 法 


苍 了 00 一 11 2， -#7); 则 可 了 
B=D (Lt, f=D'b 6.3) 
其 中 口 为 4 的 对 角 线 元 素 构 成 的 对 角 算 了 泗 , 工 . 旧 分 别 为 4 的 对 前 屿 俯 下 及 以 上 
部 分 元 素 的 反 号 数 构成 的 矩阵 . 这 样 构造 出 的 迭代 法 称 为 雅 可 比 选 代 法 ,其 用 
分 量 表 示 的 计算 公式 为 
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a 

Lk 

吕 一 > 入 
j= 


e471 i 
1 一 
[re 


定理 .3 雅 可 比 选 代 法 收 误 的 充 要 条 件 是 ptB8,) < 1 
定理 44 若 好 1 一 1 或 玫 可 -< 或 生 为 严 档 对 角 占 优 算 阵 , 或 
入 为 并 对 前 占 优 且 不 可 约 矩 阵 : 则 稚 可 比 选 代 法 收 笋 . 


1.1.3 高 斯 - 窄 德 尔 选 代 法 


类 似 于 雅 可 比 适 代 法 , 取 
B= DPD LUV, f= Dp—L)'p 《6. 37 
这 样 构造 出 的 造 代 法 称 为 商 斯 - 塞 德尔 类 代 法 ,其 用 分 量 表示 的 计算 公式 为 
五 一 Tae y》， Qj 
j=1 si"rl 


zl 一 : 二], 2, m,n (6. 8) 
Ee 


定理 56.5 高斯- 塞 德尔 远 代 法 收 或 的 充 要 条 件 是 p(B.) < 1， 
定理 6.6 若 1‖ Bi 二 1 或 上 有 1: 天 1 或 4 为 严格 对 前 占 优 矩阵 ,或 
入 为 弱 对 角 占 忧 且 不 可 约 矩 阵 , 则 高 斯 - 塞 德 尔 选 代 法 收 敏 . 


1.1.4 和 解 大 型 儿 琉 线性 方程 组 的 膛 次 超 松 弛 选 代 法 (SOR 法 ) 


站 《6. 4) 


A 


对 高 斯 - 蹇 德尔 千代 法 作 收 正 , 取 送 代 矩阵 为 
b= Dw DituD, f=wpD—uy pbp (06.7) 
其 中 必 为 松弛 因子 . 这 样 构造 的 选 代 法 称 为 逐次 超 松 弛 和 闪 代 装 (CSOR 法 ). 该 方 
法 的 松弛 内 子 m 可 根据 收 化 速 康 选 取 最 佳 的 这 sw SOR 法 用 分 量 表示 的 计算 公 
起 为 
zh 一 Ga 十 一 oa asst) 


i 1 一 这 


《6. 各》 
定理 6.7 SOR 只 代 法 收 敏 的 充 要 条 件 是 ptL.) < 1， 
定理 6.8 车 ‖E.1 1 或 上 1: 之 1, 或 坟 对 称 正定 上 且 0 加 二 2， 
或 雪 严 格 对 角 占 优 旦 0 之 w 志 1, 或 和 故 弱 对 角 鼎 谍 且 不 可 约 且 0 过世 之 1 出 SOR 
先 代 法 收 襄 . 
定理 6.9 SOR 谈 代 法 收 襄 的 必要 条 件 是 0 < 过 赤 过 和 2 
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1.1.5 分 块 选 代 法 
将 系数 年 阵 4 及 右 端 向 量 扣 分 块 , 可 得 拨 雅 可 比 选 代 法 


a 

A bh Ax i=1,2, ,gqg (6. 09) 
rm1 
| 


块 SOR 千代 法 CRSGORY> 


i 日 
A 
sl J—i 


(6, 10) 
和 舞步 求解 时 都 过 解 低 阶 方 鹅 组 凡 冰 得 小 块 x 让. 
定理 看 10 若 和 为 对 称 正定 年 阵 ,月 < 一 芭 拓 2. 则 BBSDR 迁 代 法 收效 ， 


1.2 习题 全 解 


Bx ! 2rs + 二 = 一 12 
1. 3a A + 4Ts + 2 = 20 


2 一 3 人 -10rs 一 了 


{1) 考察 用 难 可 比 秋 民法、 南 斯 - 究 苇 尔 迁 代 法 解 此 方程 组 的 收 部 性; 

(2) 用 雅 可 比 先 民法 及 高 斯 - 守 德 汞 迁居 法 解 此 方程 组 , 受 求 当 xt 一 
XI 时 迁 先 终止 ， 

解 (1) 因为 拭 阵 盘 严 格 对 角 圳 优 , 帮 雅 可 比 选 代 法 与 高 斯 - 塞 德尔 选 代 
法 均 收 禹 . 

《2) 雅 可 比 法 的 过 代 阔 式 为 


T= 2 全 名 ] 全 后 ]2 
J 了 i 

tl 计 可 ce 二 5 
四 1 3 ， 3 
Et LE! 
[ 5 1 ti6 


取 x" 三 {1,1, 1) 和 迷 代 到 17 次 达到 精度 要 求 
XY 一 《一生 000018 6，2.999 991 5. 2. 000 001 2)7 
高 斯 - 塞 德尔 法 的 进 代 格式 为 
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i 1 12 

[= 一 
LT] 1 - 

| 和 1: 一 二 2 1 本 开 十 5 
， 1 3 ， 3 

1 CE: i 

531 To ti6 


最 x 一 1 1) : 造 代 到 & 次 达到 精度 要 求 
xX = OD O18 6, 2.999 99] 5，2.0009 001 277 
2. 证 方程 蛆 
[el 十 日 .dr 0. 4 = 1 入 十 2 一 2.7: := 1 
C1 3d 十 x 二 1 
dr Br | = 3 四 二 2m -ms 一 1 
试 者 察 解 此 方程 组 的 雅 可 比 千代 法 及 高 斯 - 塞 德尔 移民 法 的 收效 性 . 
解 (1) 雅 林 比 法 的 挝 代 夭 阵 


0 一 总 .4 一 看 .条 
吾 ; 一 及 (一 丰 ) 一 0 总 oa 


一 总 4 一 站 0 
上 好 一 号 ;= (A—0.8( +0. Ba—D. 32) 
eB 二 1.093 21, 故 雅 洒 比 竟 代 法 不 收 襄 - 
高 斯 - 塞 德尔 法 选民 和 矩阵 


0 一 00.4 ~~0.4 
B, = tD— LU = 。 0, 16 -ea 


D 0,032 0.672 
otB y= FRI..=0.8<1 
故 高 斯 - 塞 德尔 迭代 法 收 敏 . 
《2) 雅 可 比 法 的 适 代 矩阵 


0 一 2 2] 
B:=D:(L-U) -| 1 旧 一 1 
2 .2 ol 


i B= ntBi y=0<1 
故 牙 可比 选 代 法 站 父 . 
高 斯 - 塞 德 尔 法 的 先 代 算 阵 
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一 二 2 
B.= Dp— LU = i 2 :| 
| 2 

} 


| 好 一 吾 |= AA 2 , otB. 

克 高 斯 - 塞 德尔 选 代 法 不 收 笋 ， 

3. 求证 lim4, 一 入 的 充 要 条 件 是 对 任何 向 和 本 x ,都 月 limAx = Ax. 

证 明 ”必要 条 件 。 由 lim 一 不. 知 lima 人 = aa: 从 而 有 | 4 一生 一 0 
(一 cc) 页 对 Yx, 有 

上 一 起 ， 过 | 二 一 总 1 二 一 全 【上 一 

即 Ax Ar Lim ex 一 4r- 

充分 条 性 对 人 rERR 有 汪 r 一 Frf 下 -cc 了 服 和 2 10， 1， 


0 0 


x = a a a A Rr 20) 
A a dir oan)! 
艳 站 人 


即 各 -A 1im4。 一 上 入. 

4. 设 Ax 二 是 ,其 中 及 对 称 正定 , 问 解 此 方程 组 的 雅 可 比 先 党 法 是 否 一 定 收 
鼓 ? 试 考察 习题 2(1) 方程 组 . 

解 涉 一 定 : 因 短 阵 的 谱 半 径 ef 及) 无 法 确定 是 否 小 于 1， 

对 局 题 251) 各 对 称 , 有 并 0 一 | 一 有 2967>> 
0, 故 各 正定 , 扯 其 雅 可 比 选 代 法 不 收 敏 . 

ddr] -ra 一 1 
5 用 SOR 亨 法 解 方程 组 4 一 十 dr; 一 了 一 4 【分 别 取 松弛 固守 ww' 一 
一 了 + +r = 一 3 

1.03, w=1. w= 1.1) oY 
03， 人 友 二 1 也 一 1) 畏 靖 解 工 = (二 1. >) .lx x | x 
“时 选 代 续 止 ,并 了 且 对 痊 一 个 昼 情 确定 选民 次 教 . 

解 ”SOR 方法 的 先 代 格式 为 


| = ri wl 


| 
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当权 妈 一 1.03, 韦 什 #2 二 (0.0, 0 时, 进 代 5 次 达到 精度 要 求 , < 一 
(£0, 500 DNd 4, 1.000 O01 6, 一口.499 999 7)7. 
当 了 避 也 二 1. 初 值 x 二 (0, 0, 0) 时 ,和 进 代 5 次 达到 精度 要 求 ,x 一 
(0.500 003 8，]1.000 001 3, — 0. 499 999 5)7, 
当 取 名 二 1.1, 初 值 x = C0, 0， 0) 时 ,迭代 5 次 达到 精度 要 求 ,x 一 
Cn, 500 N03 6, 9.999 998 5, — 0.500 000 09)7. 


上 e107 时 选民 图 目 ， 


解 ”SOR 选 代 格式 为 


Bx] 十 .23 十 x 一 一 上 
6. 用 SDR 方法 解 方程 组 ( 取 志 一 0.974 一 i 十 47; 十 2 汪 20， 要 束 当 
BT 3 FF lOr = 3 


TD = Fl 12 工作 : I : 了 ) 
本 ws tT 一 x -x ) 
Tl 二 亩 一 5 {til + 高" zx) ) 
取 初 始 值 x 一 (1 ,1.1)' 计算 如 下 表 所 示 

上 地 人 2 Te 
t 1 0 J 一 0 
1 | — 2. 600 OM 站 3.563 A000 1, Bon 550 全 
2 — 4.027 4990 3 140 65 2 2,022 B22 4 
3 一 4.057 281 4 二 .0 加重 吕 1 2,010 121 9 
4 | 一 .0D4 235 4 32.393 572 3 了 DOD O42 了 
5 -3,098 119 3 2, 999 ?6E 2 1.399 01 3 
6 ,000 G54 2 3.000 233 让 | .39 和 91 生生 
站 一 4.000 042 4 3.000 031 4 2.000 012 2 
吕 — 4.000 O17 7 2, OA Mad 了 2 O00 O02 7 


国 x 一 一 0.000 037 7 之 ]971, 大 取 x 六 


2.999 993 7, 2.000 O02 7)7, 


《一 4,000 O17 7， 


4 数 审 分 析 全 析 精 解 


7， 设 有 方 种 组 入 Fr -= ,其 中 各 为 对 称 正 定 阵 , 渤 线 公 式 
Et = x 二 reth oe Axr':}, 上 直 二 人 1, 2 


该 证 明 当 日 二 也 二 襄 时 上 述 选 民法 下 证 ( 其 中 习性 a 开 (A) 所 国 ， 
证 明 。” 选 代 格 式 茧 改写 为 
1 
故 渤 民 答 阵 8B 二 了 一 ww 入, 其 特征 秆 p 二 1 一 wa tA). 
由 [| 二 1 .1 一 wh(A4) | 二 1 得 


2 
0D < Ww A 
大雪 和 二 名 二 二 时 ,有 0 二 二 ,区 而 有 1 二 io0B) <1, 纤 出 的 选 代 
所 
格式 收 禹 ， 
8. 证明 矩 阵 


对 于 一 总 二 ec1 是 正定 的 ,而 雅 可 比 选 代 只 对 一 十 < a < 二 是 收 敏 的 ， 


证 明 当 一 广 之 4 <1 时 , 自 
det | 
出 本 

故 东 是 正定 的 . 又 雅 可 比 法 儒 代 拓 阵 


DO -qa —a 
B; = 区 得 
a ~—u 日 


二 tt —B)= la a al=i da +20 = A 一 a) (A Le) 


du 
|= la DD, detltA)= (1—a)(l+2a) 0 


喝 如 
玖 ptBi) 一 | 2a |, 即 只 有 一 二 之 4 之 让 时 , 雅 可 比 先 代 法 收 做 


9, 给 定 迹 代 过 程 x 二 Gx 十 训 ; 其 中 寂 蕊 恨 全 二 让 12、 , 试 证 
明 :如 果 避 的 特征 慎 4.(G) 二 0 一]，2,，…, 7), 则 此 迹 忧 过 程 最 多 竟 代 次 妆 
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数 于 方程 组 的 解 . 
证 明 ”出 答 降 6 的 车 当 标 准 形 了 .存在 可 逆 阵 卫 合 


G= Pi:JP 
由 于 如 的 特征 值 全 为 堆 , 故 托 - 定 有 如 下 形式 ， 
mF 1 9 1 
了 站 | 
了 一 ， 了 一 i 
了 中 


方程 组 x*， Gx 一 # 等 价 于 {fT 一 Gx 一 gg; 由 于 A10) 一 站 , 故 A 一 他 ) 一 1 一 
At 二 1 关 人 D0, 从 而 I 一 上 非 奇 异 , 即 tf 一 GYx 二 g 有 唯 -…- 解 x'. 于 是 


x Gr 1g 
与 所 述 送 代 格式 相 减 ,有 
Xx = x x 
故 x rr 
又 GC =P IP=0 
故 Xx 一 x 一， 二 x "=x 


因此 ,人 性 取 .x ,至 多 揭 代 次 即 可 收 化 到 方程 组 的 解 . 

10, 访 和 上 各 为 严格 对 衣 占 就 阵 , 且 站 < 四 扫 ]. 求 证 解 Axr 一 占 的 SOR 淡 代 法 
收 玖 . 

证 明 国 和 4 严格 对 角 占 优 , 故 as 天 DG 一 1, 2 nm) , 且 各 非 河 异 . 

利用 SOR 人 选 代 法 

LL: = Dw 1 — wD | wd) 

其 中 各 一 用 一 工 一 号 ,而 鼎 . 一 工 , - 品 分 别 为 志 的 对 角 . 严 格 下 三 前 与 严格 上 
和 角 ， 只 需 证 明 0 才 本 硫 ] 时 ,olIL) < 上 即 有 结论 成 并. 

用 反 证 法 没 L. 有 -个 特征 什 4 满足 14133 1. 则 有 

dettaf* IL.}=0 
具 面 有 
det{ (D— m7 | Du Din |)= 0 


即 dottD ~ wl: set| (1— Tl “Dp— i -vo | 日 


因 4 严格 半角 占 优 , 故 de 万 一 到 由 
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(0 一半 G 一 ie))9 一 号 一 到 [出 


1 w|i | 六 (170 le | 六 
下 有 Er ial 上 2 


口 六 -1 本 
要 ; 本 TE 的 ] 1 
丸和 | fT dg | 
7 一 ] j=1 A il 


i= 1, 2 


如 


a 


这 表 戎 亿 在 0 声名 过 1 时 也 严格 对 角 占 扰 , 吉 detC 关 0. 这 与 det( 对 一 一 
0. 子 盾 , 故 扔 设 不 成 立 ,从 而 | | 二 1, 即 pL )》 过 1,SOR 揭 代 法 收 化 . 
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第 7 章 ” 非 线 性 方程 求 根 


1.1 重点 、 难 点 全 祈 


1.1.1 问题 简介 


求 单 变 县 非 线 性 方程 
Cr 一 站 (7.1) 
的 根基 指 求 x* (实数 或 复数 } ,使 得 rr ) = 0. 称 了 为 方程 (7.1) 的 根 .也 称 
zx” 为 函 熬 ftir? 的 零点 . 车 fix) 可 以 分 解 为 
fri = rr "g(r 
其 中 mr 为 正 吾 数 , 且 glx') 关 人 . 则 ”是 方程 (7.1) 的 根 , 当下 一 工 时, 称 了 
为 单 根 ! 当 m1 时 . 称 xx” 为 ms 重 根 . 蔡 g(x) 充分 光 汕 ,rz” 昆 方程 (7. 1) 的 闸 
重 根 , 则 有 
Fr 

车 Ft) 在 [a, 5_ 上 连 总 且 Ca} (二 全 ; 则 方程 (7.1) 在 fz, 如 内 至 少 有 
一 个 实 根 .区 .a. 的 为 方程 (7. 1} 的 有 根 区 间 . 可 通过 函数 作 图 法 或 乏 侈 搜索 法 
求 得 有 根 区 间 . 

1.1.2 方程 求 根 的 几 种 常用 方法 

1 站 分 法 

设 赤 z 在 [at 上 寺 连 续 ,， Fear cD 则 六 re) = 在 rap 内 有 根 x*. 
再 设 f(x) 一 0 在 (e, 4b) 内 仅 有 -个 根 . 令 a = 6 页 一 呈 计 算 居 = 二 (eu 十 


名) 和 六 zc) 车 0) = 二 0; 则 zx” 二 xo; 结束 计算 ; 若 Fian) Firo) 人 0 网 令 
= Y= ,得 新 的 有 朴 区 间 ! ai 而 ji 车 了 a0) fo) 之 站 则 邻 二 gn， 机 二 


Xu'+ 得 新 的 有 根 区 间 ial ， 记 浅 [a 3 页 | 六 [Tan Bo ds BL 一 二 六 (6 一 uri， 此 今 
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XT] 二 a 十 如) 计算 fir;), 同 上 法 得 出 新 的 有 根 区 间 -2: ,办 .如 此 反复 进 


行 , 吕 得 :. 有 根 区 间 套 
人 [ay 和 [ay 人 [ay 


Ea Lr bh ot 1 a Dp, 1 一 上 1) 一 … 一 去 ( 


ar 鼓 


mkas 一 ao) = 0， lim = im (as +h) or 


因此 ,x 一 去 (as 一 名) 可 作为 A(z) 一 0 的 近似 根 , 且 有 误差 估计 


2. 迁 代 法 
将 方程 式 (7.1) 等 价 变形 为 
I= ptr) 《7, 3) 
其 要 求 x” 满足 fc) 二 全 则 -glr' ;上 反之 亦 然 . 称 z 为 函数 zfr) 的 
一 个 不 动 点 . 求 f(z} 的 根 等 价 于 求 Ptz) 的 不 动 点 ,由 式 (7.3) 产生 的 不 动 点 
迹 代 关系 式 ( 也 称 简单 迭代 法 } 为 
stl 一 他 ra 二 一 站 ] 2 1 
酒 数 ptr} 称 为 造 代 范 数 . 如 时 对 全 何 和 二 aa, 们 ,由 式 人 7. 4 产后 的 序列 1 
有 极限 


timr = 
则 称 先 拓片 程 (7,4) 收藏 . 
定理 7. 15 不 动 点 存在 性 定理 ) 恋 lry E CTw, 太 , 满 足以 下 两 个 条 件 ; 
1 六 任 意 x EE La 如 用 gtr) 有 是 
2 存在 正常 数 工 过 1. 使 对 任意 x, vw & za, 们 ,都 有 有 
机 {7.5) 
则 gee) 在 :ae. 的 上 存在 唯一 的 不 动 点 了 . 
定理 7, 妆 不 动 点 选 代 法 的 全 局 收 化 性 定 蛙 ) 设 ytx) E CLa, 加 满足 定理 
?. 1 中 的 两 个 条 件 , 则 对 任意 七 【ay ;由 式 (7.41 得 到 的 选 代 序 列 ! 收藏 
到 gtr) 的 不 动 点 ,并 有 误差 知 计 式 


| ri Oo re (7.67 
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上 
利 1 Ter | elm .| £7. 7) 


定理 7,3t 不 动 点 移民 法 的 局 部 收 化 性 定理 )” 设 为 gtr) 有 不 动 点 ， 
多 1 在 rr 的 某 个 邻 域 连 续 , 旦 | (Cr) 1 二 于 则 选 代 法 (7.4? 局 部 收 襄 ， 

收效 阶 的 概念 。 设 选 代 过 程 rr = glx? 收藏 于 方程 <= pr) 的 根 工 和， 
如 果 造 代 误 差 @ 二 和 一 工 ” 当下 一 时 成 立 下 列 朵 近 美 系 式 


[i 


CC 《 漠 数 个 半生) {7.8» 


别称 该 选 代 过 程 是 p 阶 收 得 的 ,特别 地 ,p = 1 时 称 线性 收 谨 ,p> 1 时 称 超 线 
性 收 但 ,p 二 2 时 称 平方 收 租 . 
定理 7.4( 收 竹 阶 定理 ? ”对 于 达 代 过 程 和- = gtxw) ,如 果 g'"”(x) 在 所 求 
根 ” 的 邻近 连续 ,并 且 
的 


) (7.9) 
gt ) 天 ! 
贿 访 选 代 过 程 在 点 =” 的 鱼 近 是 所 阶 收 喜 鬼 :并 有 
DJ Dre 7 ) (7. 10》 


斯 带 芬 森 (Steffensen) 选民 法 “ 当 不 动 点 适 代 法 (7.4) 只 有 线性 收效 阶 ， 
苏 至 于 不 收 人 第 时 ,可 用 斯 带 劳 森 先 人 民法 进行 加速 . 具体 公式 为 
二 
{7.11» 


Ee 
ba 
| 
~ 
Ee 
江 
| 
c 
上 
ra 
~ ” 


7 


此 法 也 可 写成 如 下 不 动 点 选 代 式 : 
i 


[Letey 1 (7, 12) 
PO) tr | 


定理 7. 5( 斯 蒂 芬 群 迭 代 收 第 性 定理 ) 设 .r” 为 式 (7.12) 中 和 肖 代 函数 风 x) 
的 不 动 点 . 则 x” 为 gz) 的 不 动 点 ; 设 P tz) 存在 g(r) 了 1: 则 x* 是 g(r} 的 
不 动 成 , 且 斯 蒂 苏 森 壬 代 法 (7. 11) 是 2 阶 收 得 的 ， 

3. 牛顿 球 代 法 

牛顿 迁 代 法 是 -种 特殊 的 不 动 点 选 代 法 ,其 计算 公式 为 


Firey 
Far) 


二 


下 二 D1 2 (7. 131 


Tl 2 
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其 迭代 捕 数 为 
fr) 
Fo) 


后 顿 选 代 法 的 收 区 速度 当 六 r = 一 D tr ) 关 0,(x") 关 D9 时 , 容 


易 证 明 ,ortr ) 二 Derta ) = 三 关 0. 由 定理 7,4 知 , 牛 顿 选 代 法 是 平方 


Pr) = 


收 般 的 , 且 


[| 


重 根 情形 的 牛 椒 选 代 法 ” 当 r"” 是 fir)y 二 0 的 m 重 根 Lm 这 2) 时 .选民 消 
数 r(z) 一 < 一 外 生 在 x 处 的 导数 gC") 1 一 上 关 0, 且 gC) ll 
所 以 牛顿 选 代 法 求 重 根 只 是 线性 收 筷 . 车 x 的 重 数 mr 知道, 则 入 代 式 


re = x m 0 《7. 151 
二 


求 重 根 二 阶 收 敏 . 当 由 末 知 时 ,xz* 一 定 是 请 数 m(r) = -2 的 单 重 夫 点 ,此 时 
先 代 式 


wi EC 二 Fr)f (Cr) 
A CE) [CF Er — For) fF Crs) 
R=0, 1 2. [7. 16) 
也 是 二 阶 收 雍 的 . 
简化 牛顿 法 ”如 下 选 代 法 : 
x 一 pat k=0 1 ae 
称 为 简化 和 牛顿 法 或 平行 弦 法 . 
牛顿 下 出 法 ”为 防止 先 代 不 收藏 ,可 采用 牛顿 下 出 靶 . 
4. 台 蕉 法 


将 牛顿 选 代 法 567- 137 中 的 产 () 用 fo) 在 rc， 处 的 一 阶 差 商 来 代 
替 . 即 可 得 强 规 法 


flr) 一 
Tl Ts FF 7 i e717) 


定理 7.6 假设 fir) 在 其 零点 x "的 贸 域 4: | fr 一 r” | 所 8 内 具有 二 阶 连 
续 导 数 , 且 对 任意 xzE& 有 ftir 和 关 0: 又 初 值 ,7 E 4, 则 当 邻 域 各 充分 小 
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时 , 弦 戴 法 7. 17) 将 接 阶 户 二 1 填写 v1. 618 收 各 到 了- 这 里 是 方程 4: 一 4 


2 

一 ] 一 总 的 正 根 . 

5. 抛物 线 法 

芒 截 法 可 以 理解 为 用 过 (ze ， Fa 《xe， fxr2) 两 点 的 直线 方程 的 
根 近 位 代 蔡 产 z) = 5 的 根 . 者 已 知 六 rr 一 0 的 三 个 近似 根 .rm zz 用 
过 (人 Fe Cre Core za rn 1)) 的 抛物 线 方程 的 根 近 似 代 普 
了 tz) 二 浊 的 根 , 所 得 的 壬 代 法 称 为 抛物 线 法 ,也 称 密 勒 (Miiller) 法 . 

当 Fr 在 z' 的 种 近 有 三 阶 连 匀 导 数 , 产 Cx*) 冯 0. 则 抛物 线 法 局 部 收 襄 ， 
且 收 敲 阶 为 产 一 1 839 = 1. 84. 


1.2 习题 全 解 


1[， 用 二 分 法 求 廊 程 一 x 一 1] = 二 但 前 征 要 ,要求 误差 小 于 0,05. 
解 ” 今 让 让 = 二 了 一 x 一 世 则 有 (1 == 一 1 0 人?) 二 1 癌 0;, 坡 [1 ,25 


为 函数 A(x) 的 有 根 区 间 ,又 产 (x) 一 27 一 1, 故 当 0 < 之 了 之 十 时 .C7) 单 威 ， 


当 c 这 方 时 f(x) 单 增 . 而 /( 去 ) 1: 由 单调 性 知 f(r) 一 
0 的 瞧 一 正 恨 .x”E (2). 根据 二 分 法 的 误差 估计 式 !7.2)? 知 .要 求 误差 小 于 


5 . 
d 人 


0.05, 只 雷 吉 二 0.05, 解 得 £4 1 > 5.3232, 故 至 少 应 二 分 5 次 . 具体 计算 结果 
见 下 潜 : 

点 ER I be {| Th mt 的 符 导 

0 } 2 1.5 一 

1 1 .5 2 1.735 十 

1.5 1.75 1. 525 十 

3 1.5 1. 625 1, 562 5 一 

入 1.56B2 3 1, 625 1. 53593 75 一 

加 1. 593 75 1.€25 ]. 6 S75 一 


Bx a 一 1. 609 375. 
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]. 二 附近 的 一 个 根 ，, 设 竺 方程 朴 写 成 下 


2. 为 求 方程 阅 一 归 一 1] 一 站 在 加 
列 等 价 形 式 , 并 建立 相应 的 选 伐 公式 : 


(1) 一 1 小. 迁 代 公式 ml 十 寺 ， 
-TT Tg 
ws 


一 ] 十 x ,更 懂 公 碟 Ti-| 一 


(2) 7 
a .+ _l J ; 1 
37 二 一 一 . 搂 代 公式 ill 二 一 

Xl1 vm 1 


斌 分析 每 种 选民 公式 的 收 贷 性 ,并 选取 一 种 公式 来 出 具有 四 位 有 北 数 字 的 这 你 


根 . 
解 厂 z 一 1.5 的 邻 域 -1.3,，1.6 来 考察 . 
(1) rE [3, 1.6] 时 ,glx) = 1 7 去 E [1.3. 1,67, | ¥ (x) |= 


一 一 工 < 上, 孝 办 代 式 .a,, 一 1+ 评 在 [1.3. 1.67 上 整体 收 策 . 


1. 3 


(2) 当 rerl.3，1.6] 时， 
et 一 人 上 人 3，16 


a 


= |= 3 | 
3 11 


(1 地 


: . 
| Ft 二 可 


才 zi 二 3 工商 在 [1.3, 1.6] 上 整体 收 敏 . 
-1 |.. 1 
acr— 1 ?| 人 


= 


《3) g(r) = 


一 


一 一 一 发散 ， 
ws 一] 
沽 为 选 代 公式 (2) 的 世 较 小 , 故 到 公式 (2》 进行 选 代 . 雪 求 结果 具有 由 位 有 


痪 数字 ,外需 
上， 国 ， 
| Er 工 = TI 了 = X10 
即 
1 ! 上 < 十 1 


取 这 : 二 1.5 计算 结果 岂 下 表 : 
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A 


A 


1. 481 248 D3 
.44720573 


1, 468 B17 314 


由 二 和 一 如 < 二 x 110 , 故 可 取 汪 


= Ti4 — 


1, 4656. 


1. 465 B7 


1.467 OFF 973 


1.456 243 O10 


76 B20 


3, 比较 求 e' 十 10r 一 2 二 0 的 根 到 三 位 小 数 所 需 的 计算 景 : 
{1} 在 区 间 [0., 1- 内 用 二 分 法 : 


(2) 用 选 代 法 %-， = 了 ce , 取 初 值 了 一 0， 
解 (1) 因 z EL 70 一 一 DOA 一 10.7183 之 0, 静 0 有 
2 人 1 用 二 分 法 计算 结果 见 下 表 ， 

中 EE bs rT 4) 的 符号 eT 
n 0 1 [Es 十 站 .5 
1 0 ,5 0. 25 十 | 0. 25 
2 总 DQ, 25 .125 十 0. 123 
3 人 所. 1253 .DB2 3 DV. O62 5 
4 0D. O62 5 0, 125 .093 75 十 0.031 25 
.2 3 ND. N03 75 QO. O78 123 一 1 .013 B25 
各 站 ATH 125 .3 ?5 QBS S37 5 一 | ODOT S12 5 
了 .5 947 3 四 ,93 5 NNRB9 B43 73 一 他 .由 有 全 电站 用 25 
8 .89 B43 了 5 站 -站 及 3 了 5 DO. D091 796 B75 一 站 .DO1 953 125 
9 D089 BA3 73 | .O91 795 B75 1 0.0D90 820 312 + D. O00 97B 562 
ln | 0.089 843 73 | ng B20 942 | .N090 342 31 一 | .DODD 488 2%1 
11 | ND. 090 332 431 | D0.Na0 B20 312 | 0.090 5376 171 十 : ,OOD 244 14 
12 [00.000 332 031 | 0.090 576 171 | 0.090 454 101 一 人 1 122 07 
13 | 0.090 434 101 ] 0, 090 3S76 171 | .090513 13§ 一 0. Dno Del 035 
14 | 0.090 515 136 | OO.090 376 17110.090 515 B53 十 0. O00 30 517 


104 数值 分 析 全 本 楼 名 


二 分 法 迭代 14 次 后 满足 精度 要 求 , | zu 一 r” 宅 Fis 一 0.000 030 517 < 也 XX 
10" ,x se xu 具有 二 位 有 效 数 字 . 

(2)7 当 rE [Lo 0.5] 时 ,Pr) € Ln, 0.5], gm 一 市 -= | 去 工 一 
0. 825， 故 送 代 式 tt 二 证 (2 一 em ) 在 [0，0.53 上 将 体 收 敏 . 取 rw 一 0, 选 代 
计算 结果 见 下 表 : 


中 LI 中 A 


OD. om S12 816 


D.DBS 182 908 .O90 326 468 


9.090 639 135 .O90 321 831 


选 代 6 次 后 满足 精度 要 求 x 一 zx" | 护卫 7 | 一 | 坟 0.000 007 20 二 


二 x 107, 故 zs x 精确 类 三 位 小 数 ， 

. 给 定 画 数 (zx), 设 对 一 茹 了 ,f(r) 育 在 且 0 忆 肋 王 天 (zz) 开 角 ,证 明 对 
于 范围 0< 二 六 内 的 任意 定 教 1, 选 伐 过 程 .c 却 工 一 AFCT4) 均 收 项 于 六) 
一 曲 的 根 二， 


证 明 由 于 te) > 0,f(z) 为 单 增 两 数 , 故 方程 fx) = 0 的 根 x+” 是 叭 
一 存在 的 (假定 方 种 有 根 +*). 连 代 本 数 gt) = 一 工 一 Ar，| g(r? | 一 
1 一 f(z) | 由 0 之 加 各 六 (之 导 及 0 入 得 ， 0 
M2 一 1 |) | 过 工 一 
maxf | 1 一 Aa ，:1 一 Wf} -<< 1. 由 此 可 得 

| oOo Lr | 0 ) 
fimx, 一 

5， 用 斯 蒂 芬 森 选 抠 法 计算 第 ? 题 中 (2):13) 的 近似 枢 , 精 确 到 10“ 

和 解 ” 记 第 ? 题 中 (2) 的 选 代 函数 go (xc) = (1 十 对 这 ,03) 的 迭代 函数 为 
,利用 选 代 式 (7, 11) ;计算 结果 见 下 表 ， 


PitT) 一 1 
wr—l 
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中 加 速 gzt7) 的 结果 必 上 加 速 ps tz) 的 铺 果 x 
| 1.5 加 1.5 
1 | 1. 465 558 485 1 1. 467 342 286 
1,. 65 S71 233 2 1.d85 B76 OSS 
3 1. 465 371 232 3 1. 465 571 232 
| 1 L485 53571 232 


6 访 旬 一 工 一 疡 (zz 一 GEr) PCr) 试 确定 画 吉 pir) 和 g(r ,使 求 
解 .AKCz) 一 品 且 以 gz) 为 选 伐 函 数 的 选 伐 法 至 少 三 阶 收 开 . 
解 过 求 :or = glz) 二 阶 收 襄 到 fir) “0 的 根 x" .应 有 Flr "=r, 
gp Cr" ) 二 0. gtr") 一 0 于 有 是 由 
pl I Fr ) 
Fr ptr df tr 
Fr = 2p FG pelr Fr) 2 Fr I=0 


i LF 
得 pir ) Fry gtr") 2 Tr 
故 取 
= yD 
piry FE’ Ei 和 过] RE 


即 送 代 至 少 必 阶 收 襄 . 

7. 用 下 列 方法 求 作 T) 二 沾 一 3 一 1] 二 0 在 zx 一 2 附近 的 很 ,和 根 的 准确 
秆 ”一 1.879 385 24… ,和 要求 计算 结果 准确 到 四 位 有 效 数字 ， 

(1) 用 牛顿 法 ， 

《2) 用 弦 截 法 , 取 7 二 2, x) 二 1.9; 

3) 用 搜 炳 线 法 ; 取 谨 = 二 1 二 4,ry 一 2， 

解 = FO 0 flr) = m3= rr-1) 0. 
f= br 0, 财 Yrie!ll. 2 ， 

C1) 有 联 以 一 2, 用 牛顿 选 代 法 


了 一 3 一 ] 2 
3 下 一 3 3 一) 


Xa 
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计算 得 ri = 1. 888 889,x: 二 1.879 452，| ze 一 ”| 之 


Te = 1.879 452. 
(2) 取 = 2, x 1.9, 利 用 东 堆 法 


(一 Fe 
FE 


EE TE 


得 ,r: 二 1.98] G94, x 二 BRO 84l, ri = 1:879 489; [| re :< 汪 x 


]073, 故 取 == x -= 1.879 489. 
(9) rm .7 二 3, 二 2. 抛物 线 法 的 近代 式 为 
2 fur) 


wi sign{r) wa 二 Fr 并 人 ri 


区 | Tk 


w= fe Tes] Flr: Tas Kes {pO— ol) 
达 代 结果 为 :rs 二 1.953 968,， z， 二 1.878 015， x 二 上 .879 387 已 达 四 位 有 效 
数字 . 


8. 分 唱 用 二 分 法 和 和 牛顿 法 求 tanr 一 中 的 最 小 正 契 . 


解 显然 x” 一 0 满足 x 一 tant 二 0. 另外 当 | x: 较 小 时 ,tanz < 十 六 站 
十 … 十 元 +…, 故 当 工 毛 fo， 下 时,tanr 2 因此, 方 种 天 一 tmz 一 站 的 


最 小 正 根 应 在 区 间 ( 等 ， 二) 内 


记 fw) = rotanr, tt { 竺 ， 袜 x) ,容易 算得 A4) = 2.842 > 0, 


0 [4.6] 基 x) 二 0 的 有 根 区 间 . 
对 于 二 分 法 ,计算 结果 见 下 表 : 


Eq a bs Ik 六) 的 符号 
0 4.0 由. 五 4.3 4 

上 二 .3 4.6 4. 45 十 

2 .45 4.6 4, 525 一 

3 十。 和 4. 525 14.487 3 十 

4 497 5 1. 525 4, 306 23 

5 44. 二 人 了 5 二 5 和 23 4. 496 B73 一 
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续 表 

9 a ba Tk Fr 的 符号 
6 4, 487 3 4. 406 B75 4.492 187 3 十 

了 4d. 492 187 5 4. 496 375 4. 494 531 25 一 

只 4.492 187 5 二 ,4 531 23 4. 493 358 375 十 

9 恒 493 359 375 ,94 531 25 4. 493 4d43 313 一 


| -3 
此 时 | mx T1107. 


车 用 牛 帧 选民 法 求解 ,由 于 六 (x) = 1 一 seer 一 一 《tant)2 之 0,7(r)= 


小 ， ， 
一 2tanzr on 世人 D, 邦 取 x = 二 44. 1, 选 代 计算 于 果 见 下 表 : 
中 x 点 Ik 
1 543 7 了 32 122 生 4 和 93 412 197 
2 .SDB 1435 388 ,03 409 458 
3 4. 434 1]71 63 时 和 4 和 93 408 45 和 8 


所 以 了 一 tanr = 0 的 最 小 正 根 为 x a 4 493 41, 
9. 研究 求 Ya 的 牛顿 公式 


一 计 ( 二 + 衬 )， 2 
证 明 对 一 切 二 一 1， 2, 下, 之 Ya 上 且 序 列 x, ,xs，… 是 递减 的 - 
证 法 一 用 数列 的 办 法 . 因 z, 0 出 王 = 了 (2 十 到- ) 知 mm >>0, 且 


石和 = 1,2,3,…. 又 由 


-1 


El 三 局 


= 言 + 溃 六 二 + 下 一 1， YA21 
[| 所 单 减 有 下 界 Ya. 根据 单调 有 界 原 理 知 , {x } 有 极限 . 易 
证 其 极 根 为 Va， 


证 法 二 设 大 了 一 下 一 GE 0 十 叶 ) 网 00) = 一 e017( 二 =) 汪 


Tk 
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0D, 许 此 x) 一 日 在 [0, 十 55) 肉 有 根 . 由 后 (x) 一 2 人 0 知 , 六 xz 在 工人 0 时 
单调 增加 , 故 fx} 一 0 在 人 ,十 ce 内 有 了 唯一 实 根 va， 又 产 (r) 一 2, 故 当 了 六 
Ya 时 对 tz) 应 用 牛顿 选 代 法 ,得 


Fire) ral | 如 
PE 3 (x1 隐 )]. ££=0,1,2, 


| I 


收 敏 , 当 0 一 辣 过 wa 时 ， 
工 | ; (x + 二) ; | 去 人 ] ES 


En 


此 时 ,从 Tl 起 , {zl 单 浅 有 下 界 ya, 且 极限 仍 为 ya. 


10. 对 于 Fizy 一 0 的 牛 疾 公式 za 一 一 {7 "证 明 
3 


开工 1 
Ri = 
{rs A 了 


收效 到 一 区 "这 旺 >- 为 Flz) 一 0 的 根 ， 


证 明 根据 牛顿 选 代 法 有 
1 下 一 1 23， 


了 TL 一 工 4 -] Fry 
一 Fe) 
页 了 1 FPCz 
AR 上 | 一 1 Fm) ? = 
Crp 一 2 大 Ce fx) 


1 一 站 2 {Fr 一 


Ero} Fry PCr 


FSDTF (Cr) (rom) 
Cr ly CB) (rer rx) 


其 中 所 介 于 i 与 工 之 间 , 生 ; 介 于 忒 ;与 x" 之 间 , 根 据 式 (7.14) 得 


: EY | = FS VT rs) : XI 一 
lim {I Or) lim Fr fF C8) | re sO— "Yl 
加 
2 fF(lr') 


11， 用 牛顿 法 和 和 求 重 要 前 牛 额 选民 法 (7.15) 和 (7. 18){ 书 中 是 式 (4. 13} 和 
式 (4.14)) 计算 方程 Fr) 一 【sinz 一 持 ) 一 0 的 一 个 运 仆 要 , 准 确 到 10 7 , 初 
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怒 值 一 分 
解 fC = (sinz 一 子 ) 的 根 z” 为 2 重要 
ft) 一 2(sinr 一 子 ) (cosr 一 方 ) 
用 牛顿 法 选 代 公式 为 
(snm — Fn) 
(os 了 
a 直 二 Bl: 


2cosr, 一 上 
令 .二 到 ,出 2 一 1.785 398, x3 一 1.844 562，…， 和 迷 代 到 > 一 1.895 494， 


| 了 ”一 1.895 49 |< 107. 
用 求 重 根 的 选 代 公 式 47. 15) , 迁 代 公式 为 


Sin 一 方 二 
Ts EE - ， 让 二 0, ,2 
cosm 一 可 
取 一 到 ， 则 r， 一 2.000 000， rs = 1.900 996, r; 一 1.895 512，z， 一 


I.895 494，25 一 1, 895 494. 四 次 迁 代 让 到 了 上 而 rw 的 结果 . 
车 用 公式 (7,16) , 则 有 


Frm) {ry 
TO TT Fem Fr) 


独 了 Fi 产 (z) 及 产 (z7 一 2 { cosz 一 二) 一 2sinzfsinz 一 诗 *) 代 人 上 述 远 代 
公式 ,得 


(sinr 一 二 ) {cosz, 一 村 ) 


(cosx, 一 地 } 十 sinxre ( sine 一 本 ) 


XH] A 


取 z， 一 站 .得 x 一 1.801 749 ,wz 一 1.889 630, rs 一 1,895 474， 一 
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1. 895 494. x: = 1.895 494. 

12. 应 用 牛 额 法 于 方程 一 a 二 0, 导出 求 立 访 报 让 的 先 代 公式 ,并 讨论 其 
收效 性 . 

解 说 Fr 一 证 一 ay 六 tr 一 3 天 (rr 一 6r 牛顿 法 选 代 公 式 为 


Fr ) Ma_ ?+a 
x 轩 = 下 一 0，1，2，… 
I Fr 3 ” 


当 &a 芭 旧时, 玉 zx) 的 单 根 为 ,此 时 牛顿 法 在 根 = 附近 是 平方 收 化 的 


一 它 和 
Fit] 一 3 


则 x -x 0, 故 对 Fz EE RR 均 收 襄 . 
13， 应 用 牛顿 法 于 方程 (ri = 1 一 忆 二 0, 导出 求 Ya 的 迷 谍 公式 ,并 用 此 


5 
xr 


公 丰 求 wi15 的 值 . 


解 f(r) 一 1 年 ,f(x) 一 各 ,x 天 0, 根据 牛顿 选 代 公式 有 


Tp Ih = rs ( 3 
2 2 
Xx! 


可 知 产 (Cr) = 一 名 一 0 均 联 x (0, va 时 ,迭代 收 筑 . 


对 于 “115, 取 zu 一 9, 迁 代 计 算 ,得 
I1 = ,330 434 78， ra 一 10.702 425 53， zi 一 10.723 741 4 
74 一 10.723 805 29， rs = 10,723 805 29 
页 115 2 10.723 805 29. 


14， 应 用 牛顿 法 于 方程 f(z) 一 x" 一 4 一 上 和 f(x) 一 1 一 所 一 0 分别 性 


出 求 Y9 的 先 代 公式 ,并 求 
Ih Ya 一 fi 
tr Wa Oo) 
解 车 fx) 一 zx"~a, 则 有 (ec) 二 nr", 牛顿 送 代 法 为 


ri = 一 2 | Vat | = 0, 1, 2 
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如 _7—1 

根据 定理 7.4 知 (9) 一 | 
本 (Wa CO— ror) 1 n—1 
一 2 WA 


一 了 ,牛顿 选 代 法 公式 为 


co+D- 全 |， Re 0 1. 2 
人 此 
有 Ca 六 十 
根据 定理 7.4 知 (VO) 一 一 三 ) 
lim NE 一 I 四 1 
Fe Ca — ry? 加 J 
15- 证 明和 办 代 公式 


3a) 
“mt 30 十 已 
是 计算 Ya 的 三 阶 方法 . 候 定 初 值 rs 这 分 靠近 究 工 * 


lim YA 一 it 
kr ta — re )? 
证 朋 ”着 没 PCz) = 三 3 , 则 有 an = ql) 用 PCVS) Ya. 由 


一 Ya' 求 


gz 的 定义 ,有 


(3x 十 Ge 人 zz = rtr + 3a) 
对 上 式 两 端 连续 求 导 二 次 ,得 


Grplr) 十 (324 a (ry 一 3x2 十 3a 
69(z) 十 12zg (tr) t+ tar tay ry = x 
lB C7) 十 38xgp Cr) (3x: + ay (ry) = 6 
代 工 = Ya 依次 人 上 三 式 , 并 利用 gly2) = Ya .得 
ya) =0, ya) = 0, 
所 以 造 代 公 式 是 求 va 的 三 阶 方法 


lim 一 


] 
Er 3 


y= 二 


S| 
二 | 
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16. 用 牛顿 法 解 方程 组 


取 0 ii12) - 


解 车 i(xty) 一 下 十 一 4 f(T 二 一 六 一 1; 则 


1 

， Zz 2Yy ， 4 
下 【3 一 | | [Ft{r; vx}! = 

2 工 ~—2y 1 

4y 


牛顿 千代 法 为 


1AH] > | 
tit 
Rd 


代 初 值 "(1.5,1.2) ,人选 代 计 算 ,得 


[27 


rT 并 ， Fix, 
-re 


yt y 


人) (9 250) 本 ) (人 138 1 ) 
yn .225 000 " yt 1. 224 ?44 9 和 8 


x 1.581 138 830 rx 1. 581 138 $30 
Cs)= ( ). CG.)=( 


wil \1.224744871 y 


) 


) 


1. 224 744 B71, 
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第 8 章 矩阵 特征 值 问题 计算 


1.1 重点 .难点 全 本 


已 知 为 二 (ay js 称 pa) 一 det 一 4 为 是 的 特征 多项式 ,特征 方程 
P44) 二 dett 守 一 存 ) ~ 0 的 根 称 为 4 的 特征 值 , 相 应 于 特征 值 4 的 齐 次 方程 组 
{以 一 和 ?x 二 9 的 非 零 解 x 称 为 撼 取 生 的 对 应 于 i 的 粕 征 向 重 . 

求解 算 阵 4 的 特征 值 及 对 应 的 特征 向 量 的 办 法 有 两 类 :-- 类 是 宪法 与 反 震 
法 , 属 选 代 法 : 另 一 类 是 正 交 相似 变换 的 方法 , 属 变换 法 .如 绰 可 比 法 和 QR 方 

1. 车 法 

笑 法 是 一 神 计算 矩阵 主 特 征 秆 ( 按 横 最 大 的 特征 值 及 对 应 特征 向 量 ) 的 选 
代 方 法 . 

定理 8.1 设 和 ER "有 # 个 线性 无 拓 的 特征 向 量 , 主 特征 信和 1 满足 | 1 i 
人 > 1 社 久 1 这 下 室 1 1 加 对 任意 非 等 初始 向 量 关 -wn. 按 下 述 方 法 村 造 
的 向 量 序 列 {un} im)， 

DD 
二 A] 


pH = ma vy } ‘8.1) 


1 


Bad 


有 ely Hm 一 上 C2) imgs < 


i 
maxt x ) 
其 中 maxtw*) 表示 向 基 加 的 绝对 值 最 大 的 分 量 . . 

将 法 可 用 原点 平移 法 进行 加 速 收 敏 的 计算 , 对 于 实 对 称 矩 阵 ,也 可 几 珊 利 
商 胡 速 法 进行 计算 . 甘于 主 特征 值 是 其 他 情况 的 , 仍 可 用 罕 法 进行 计算 . 

2Z. 世相 法 
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反 天 法 用 来 计算 非 奇 异 实 答 阵 控 机 最 小 的 特征 值 及 其 特征 向 量 , 绪 合 原点 
平移 法 ,也 可 用 来 计算 对 应 于 一 个 给 定 近 似 特征 值 的 特征 向 量 . 
反 短 法 交代 公式 为 : 任 取 初始 向 量 ww 一 车 天 和 , 即 


二 . 
(8, 2) 
1 .| 


maxt ve) 
定理 8.2 设 志 为 非 奇 红 实 奸 阵 , 且 有 # 个 线性 开关 的 特征 向 量 , 其 对 应 的 
特征 值 满 足 
1 
则 对 任何 初始 向 量 ‘as 半 0} ,由 反 舌 法 构造 的 向 量 序列 iy4)， {ws} 满足 


or | 
maxt x } 


定理 4.3 设 AE R*" 有 #t 个 线性 无 美的 特征 向 量 . 的 特征 值 及 对 应 
特征 向 量 分 别 记 为 4 及 xi 一 1， 2， ,1) ,而 记 为 4) 的 近似 值 .C4 一 pp) 人 
在 , 且 :网 一 让 受 | 多 一 P| 0 关门; 则 对 任意 非 等 初始 向 量 二 v(m; 闫 9)， 
由 原点 平移 的 反 老 法 选 代 公 式 


Ve 二 CA ph wl 
8, 3 
th 一 一 全 天 一 1，2， | 


maxft 


《1》 me 一 一 1 2) limmax(ty) 一 一 将: 


产生 的 向 量 序列 {w) 、 iu! 满足 


Ti 


{1) lm max(x)’ 

(27 lmrmaxt ys: ) 一 = ma A tk 2), 

友 雷 法 进 代 公式 (8.2),(8.3) 中 的 是 通过 解 方程 组 vw 一 tI 和 (4 一 
Blw 二 机 -| 而 得 . 

3. 雅 可 比方 法 


雅 可 比方 法 是 求 一 个 实 对 称 和 矩阵 的 全 部 特征 值 和 特征 阿 量 的 方法 ,其 基本 
思想 是 通过 一 系列 正 交 相似 变 换 , 把 实 对 称 矩阵 化 为 对 角 阵 ,该 对 角 阵 的 对 角 
元 就 是 原 矩 降 的 特征 值 , 所 用 的 正 交 变 搞 为 平面 旋转 变换 (也 称 Givens 变换 )， 

记 A: 一 上 各, 雅 太 比方 法 为 

并 二 1 (8.42 
其 中 民 为 平面 旋转 变 挽 上 RCp, 9， 六 ，p， g(tp 二 四 为 声 ,中 按 模 最 大 的 非 对 第 
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元 的 下 标 ,9 的 选取 应 使 < 由 一 0. 当中 ”一 时 


cos8 一 冯 ， sing 一 signlay cos 如 
当 a 11 了 光 a 时 , 今 
ly tl ， 
一 2 二 人 + I= tang sgnCe) 
2am Icl+ ve+t+l 
1 。 
cosd = » Sing = 1eostd 
v1i+er 


这 样 就 确定 了 让 一 RC(p， gq， 用， 并 使 sg 二 ay 二 0 RTRI*…RY 的 各 列 分 别 是 
各 的 特征 值 对 应 的 特征 向 量 . 
雅 可 比方 法 是 收盘 的 . 
4，QR 方法 
QR 方法 是 用 来 求实 短 阵 4 的 全 部 特征 值 的 方法 . 它 主 要 由 商 步 组 成 , 首 
先 通过 反射 变换 (Househotder 变换 ) 把 4 正 变 要 似 变 搞 为 上 海 森 伯 格 阵 : 然 后 
把 所 得 之 上 海 森 伯 略 阵 通 过 QR 方法 变换 为 本 质 收 葡 的 上 三 角 阵 . QR 方法 
如 下 : 
4 一 点 一 外 及 
一 Ri 一 上 下 
其 中 的 QR 分 解 常 用 平面 旋转 变换 完成 . 
关于 QR 方法 的 收 敏 性 用 浸种 为 对 称 纸 阵 时 的 相关 结论 见 教 材 . 为 如 速 收 
化, 也 可 使 用 带 原点 平 称 的 QR 算法 . 当 特 征 值 为 共 固 复 根 时 ,可 用 双 步 原点 平 
称 的 QR 方法， 


k= 2 3 (8.53 


1.2 习题 全 解 
1. 用 团 法 计 贡 下 列 眶 阵 的 主 特征 值 及 对 应 的 特征 向 量 ， 


了 3 一 3 一 4 3 
4 和 if {2) A; 一 区 6 | 
一 2 一 1 3 3 3 1 


当 特 征 全 有 3 位 小 数 稳 定时 造 代 终 止 . 
解 ”套用 大 法 公式 


【1) 站 | 二 
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本 了 0 
max(t RE ) 


取 四 = 《1,1 1 关 0, 将 雪人 民 入 上 式 , 计 算 结 困 见 下 表 ，; 


点 ui maxt th ) 


1 Cl OD.75, OY 8 

2 Cl. 0, B48 648 B40, ~ 0, 297 207 2077} .3 

二 [1, 0. 608 798 347, — 0. 388 839 6817 铀 .38 二 人 R50 
时 Ch, 0. 605 776 832. — 0. 394 120 752) 9.605 429 O02 


? Cl B05 B09 732. C—O. 394 F368 924) 9. PN5 572 O02 


即 A&， 的 主 特征 信 4 = 9.605 572， 特征 向 量 x 有 (1，0.605 610， 
— 0.394 369)7, 
将 4: 代入 因 法 公式 , 取 柚 = (1 1, 1)7 ,计算 结果 见 下 表 : 


上 CE maxtv} 

1 {0, 285 714 286， 0.714 285 714, 1) 7 

2 (0, 162 790 698. 1, 0, 551 162 791) 6, 142 857 143 
5 【一 0 476 667 405, 1. 0.275 116 331) Bb, 400 967 982 
10 (0. 398 164 195, 1. 0.155 903 744} B. B55 264 397 
16 《一 0.604 221 865, 1, 0,150 937 317》 B. S69 534 947 
17 一 站 -自生 288 O82, 1, 0. 130 B81 2917 8. R69 699 412 


帮 A 的 主 特征 信 X41 xs 8.869 699， 主 特征 向 量 为 《一 吕 604 288, 1， 
0, 150 8817T. 
2. 利用 反 堵 法 束 答 阵 


5 2 1 
4 一 |23 1 
1 1 1 


的 最 接近 于 日 的 特征 慎 及 对 应 的 特征 向 量 . 
解 ”4 的 特征 伟 的 近似 值 思 一 6, 将 
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2 1 
B=A—H= 2 —3 1 
1 1 一 5 
进行 三 角 分 解 :可 一 EL, 其 中 
1 oO 0 一 3 1 
1 DD 
1 5 _1 
P=loo01|, LI=|z 1 9 pi 3 5 
] 0 0 4 27 
o 写 1 0 0 全 
根据 tm = 0 ,1，17 解 得 ,tn 一 一 上 一 , 琢 
maxtw ) 
Ly = Pair Ey, = y+ my = 
计算 得 以 下 结果 : 


站 


本 1 


Hs 
了 7 


Pb; 


= (1.618 518 519, 0. 807 407 407, 0,185 185 185)™ 
— (C1, DO. 498 855 835, 0. 114 416 475)7T, A m= 6.617 848 97 
= (0. 498 855 835, — 0. 135 O11 442，1. 108 O09 154)7 


re 0. 742 944 316, 0., 397 A06 559, 0.205 186 8877 


一 (1，0. 534 907 597, 0. 276 180 698?7 ， 2 7.345 995 895 
= (0.334 907 397, 0.008 726 899, 0. 993 O18 48)™ 
= (0.787 588 409, 0. 408 053 844, 0. 183 892 31177 


;3 一 《1，0.518 105 446, 0, 233 487 833)7 , A 7,269 698 727 


= (0.51 105 446, — 0.025 864 89, 1,. 020 451 912)7 

= 60.772 837 002, 0. 405 513 711, 0. 188 972 576)7 

= (1, 0.524 707 939, 0. 244 518 023)T, A a 7, 293 933 905 
= (0.524 707 939, — 0.017 835 946, 1,014 2683 757}™ 

= 00. 777 569 535，0. 406 086 225, 0,187 827 547)T 


i = C1, 0.522 250 689, 0. 24]1 557 235)7， > 7,286 058 616 


= (0. 522 250 689, — 0.019 568 109, 1,015 654 488)T 

= C0.776, 020 139, 0.405 957 918, 0, 188 O84 164)™ 

= 1, 0.523 128 07, 0. 242 370 209}T, A 2 7, 288 626 351 
= (0.523 128 07, —0.019 193 826, 1.015 355 O61)1 

= {0.776 528 141, 0. 405 985 642, 0. 188 028 715)T 
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u; = C1), 0.522 821 544, 0. 242 140 2d457T , A 7. 287 783 336 
观察 结果 得 ,4 的 与 8 最 接近 的 特征 值 约 为 7. 288, 对 应 特征 向 量 为 (1,0. 522 8， 


0. 242 1)T。 
3. 束 卸 阵 
4 人 0 0 
0 1 3 
与 将 征 愤 4 对 应 的 特征 向 量 ， 
解 ”和 宪 阵 有 特征 值 4,4,2. 故 可 用 车 法 求解 ,本 加 = 51， 1, 137 ,利用 筹 法 
计算 得 
真一 (444 和， 划一 《人 1， 1， 1)， maxtn) 一 4 
Ta 《44 4) e 一 (1 tiy maxfwy 一 革 
看 与 特征 鼻 4 对 应 的 特征 向 量 为 (1,1,1)". 
4.《I} 谴 各 是 对 称 姑 附和 ilzrliy > 一 1 是 各 的 一 个 特征 值 玛 相应 的 等 
征 向 醒 . 又 该 产 汶 正 充 阵 ,使 不 一 利 一 (1 0， 0)T, 证 明 如 一 B4BP7T 的 第 
一 行 和 第 一 列 除 了 4 外 其 食 元 素 均 为 零 . 


(2) 对 于 产 隆 
2 10 2 
各 一 5 - 
2 8 11 


1 一 9 是 其 特征 值 ,x 二 ( 拖 ， 直 ,性 】 是 相应 于 9 的 特征 向 醒 , 斌 求 一 初等 反 


射 阵 卫 ,使 Px 二 ei ,并 计算 BB 二 PAP'. 
解 (1;) 量 二 PAP', B87 一 {PAP )' 一 P47TP' 一 Ph4PT, 故 下 为 对 称 搜 阵 . 
又 A 二 江 , 卫 为 正 交 阵 且 Px 二 a, 放 由 PAP Px 一 P4x 得 Be = APr 一 jel， 
妈 4 是 昌 的 特征 秆 ,et 是 县 的 与 4 对 应 的 特征 向 量 . 再 由 Be, 一 je 知 rl 二 b; 
如 二 0b 一 人 根据 对 称 性 知 aas 一 让 一 … 一 和 ,一 0, 妈 吾 的 第 一 行 ， 
第 一 到 的 元 素 除 广 ， 一 4 外 ,其 余 全 为 零 . 
《2)》 取 上 = 工 一 上 作为 反射 罚 面 的 法 向 量 ,可 将 x 变 为 e1. 此 时 上 一 


(一 去 ,十 , 三】 ,反射 给 阵 
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2 1 2 

3 3 3 

1 2 _2 

r= 3 

2 .2 1 

3 3 3 
9 DO 必 
从 而 有 且 一 PAPT 一 。 18 
0 0 ~—9 


5. 利用 初等 反射 贬 将 


正 灾 相似 转化 汰 对 称 三 对 站 阵 ，. 
解 。” 对 向 量 (3,4)' 作 及 射 变换 ,使 其 变 得 与 e 一 (1, 9) 平行 ,此 时 


ov 和 3 一 4 
3 _ 4 
于 二 工人 ss 5 
1 2 十 (,)' 本 了 一 4 3 
5 5 
所 以 ,所 求 的 反射 阵 为 
1 0 0 
3 4 
H- | “3 7 
_4 3 
0 5 5 
且 
1 —5 0 
_so 73 14 
HAH 一 ” 站 25 
14 23 
0 相 一 王 


6. 设 息 -| 是 由 素 斯 审 东 德 方法 得 到 的 超 阵 ,又 设 》 是 1 的 一 个 特征 
催 章 ， 
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《1) 证明 起 降息 对 应 锥 特征 向 重 是 x 二 PP sy 
《2) 册 于 挫 出 的 了 应 如 何 计 算 x? 
证 明 (1) A | 一 中 2 了 下 和 PP- ， 设 对 应 的 特征 值 为 1, 册 4? 
一 好, 所 以 有 
下 PP 一 
ACP) Pi P29) 一 ACP, P,P sy) 
圾 x 一 户 P.…P,:y 是 和 4 的 与 特征 值 4 对 应 的 特征 向 量 ， 
<2) 因 忆 ， P; ，*…,， PP : 痢 是 反射 矩阵 且 是 在 计算 入 . ; 的 过 程 中 得 到 的 , 故 
可 志 采 用 如 下 过 程 计算 x. 今生 一 PP， 
Hi = PH k=2,3,..,n—? 


x= Hiy 
7- 用 带 和 位 移 的 QR 方法 计算 
1 2 站 3 1 0 
tj 二 |2 一 1 1|; 2) 8= |1 2 1 
中 1 3 n 1 1 


的 全 部 特征 值 . 
解 (人 2) 记 4 一 4 到 ”一 al 作为 平移 因子 来 计算 4 的 全 部 特征 值 . 
31 一 全 


.828 427 124 — 4.242 604 686 人 .707 106 781 


PisPrtAl — 1:1} = 有 R= 1. 732 050 S05 — .977 330 268 
0 0 0 408 人 24 245 
--2.0 4. 224 744 名 0 
A: 一 RPTPh+ aI= |1.22474487 1,666 666667 0.235 702 28 
0 0D. 235 702 26 4, 443 333 333 


5 O33 333 333 


3.472 151 ?17 一 1.5606989 0.052 733 495 
有 3 了 4 一 全 和 一 有 一 n 1. 370 88 834  — 0.226 301 
0 0 日. 038 502 921 
~ 2.330 49 345 0. 306 779 526 a 
上 ;一 RP PE 二 :T= NO 306 779 526 11. 978 401 S22 0. O06 7D2 831 
0 DD, O06 792 831 3.372 247 B22 


$3 = 3.372 247 B22 
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6.731 113 823 —0,380 950 572 0.000 363 611 


PiaPitAhy 一 = 有 和 = 0 1. 375 442 B92 0.000 330 107 
0 0 0, 00D 033 499 
— 2,371 041 162 0.073 625 778 0 
A = RPEPE+T r= | 0.073 625 778 1.998 760 14#5 0 
0 0 3. 372 281 32 


二 起 有 一 个 特征 植 4， 二 3.372 281 32， 对 站 ,的 子 策 阵 


2.37] 041 162 0.073 625 7787 
和 | 0.073 625 778 1.998 760 145J 
继续 进行 变换 , 取 s = 1.998 ?60 145, 得 
4.370 421 512 一 0.073 615 329 
0 — 0.001 240 327] 
一 2.372 281 308 — 0.000 020 895 
—0.000 020 895 11,998 760 45 


因此 4 的 另 两 近 估 特征 值 分 别 为 一 2 372 281 308 和 1.998 760 145. (实际 上 妇 特 


pith 一 一 太一 | 


种; = RP -sf = | 


征 值 为 十 士 荆 V 续 和 2, 只 要 再 对 和 变换 一 次 即 可 得 到 准确 值 士 一 2 和 2.) 


{2) 记 B= 二 利 . 扶 5 二 a 的 办 法 选取 平移 因子 , 则 5 一 1. 
2,236 068 1.341 541 吕 . 和 47 214 


了 《有 一 5 和 一 下 一 0 1 ,D095 445 日 355 4 
QO 好 0.816 497 
3.6 0. 489 898 0 了 
一 RPiPi + = |0,489 898 1.733 333 ~—0.745 356 
D — DO.743 356 0, S66 66FJ 


#2 = 0O, PAB B67 
,973 961 0.658 B18 ~—0,122 782] 


PPstB,—s:D) =R= 0 1,224 403 —0.583 259 
D 0 — 0. 447 537J 
3, 708 S43 0.20] 5895 0 
县 RPEPL 十 5 一 10.201 695 1.979 853 0.272 439 
O272 439 0.311 508 


3 .11 B08 
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,402 917 308 0.300 218 995 OD. 01é 147 ?47 


PaPstB Cem)= 让 = 0 1. 675 653 371 0D. 2688 341 036 
0 0 一 0.044 216 968 
-726 337 0.099 318 0 -了 
B, 一 APD PT 十 一 |0.099 318 2.005 887 一 0.007 189 
D 一 站 007 189 Db.267 979 J 
和， 一 站 .267 979 
3.459 784 0.149 160 一 0.000 2067 
PysPs (下 一 iD 一 及 一 0 1.784 156 —0.007 186 
| 0 0 — 0.,000 030J 
.730 619 0.049 ?781 昌 
B: = RPEPE, -3 一 |0.049 781 2.001 435 0 
0 0 0.267 949 


现在 收缩 ,继续 对 8B; 的 子 和 矩阵 
i [2 619 0,049 781] 
0.049 78] 2.001 435 + 
进行 变换 .得 到 
3.732 05I 0 
B; ‘= PtB, — s PE 十 sf = | | 
号 3 
故 求 得 站 的 近似 特征 值 为 
一 3.732051， 避 一 2， A = 0.267 949 


8. 试用 初 等 反射 阵 将 


1 1 1 
站 二 |2 一 1 一 1 
一 4 5 


分 解 成 f 屯 的 形式 ,其 中 也 为 正 变 阵 , 呈 为 上 三 角 阵 
解 ” 第 一 步 ,将 各 的 第 一 列 变 为 与 e 平行 的 向 量 , 取 = (1 十 2: 一 2:) 


一 3,n 一 【1.2， 2) TT 二 +gel 二 《4,2， 2)7 ,A = 方 | ul 有 = 0 = 12, 


因此 ,所 求 反 射 阵 为 
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P: = Bu 一 


一 3 3 -3 
中 站 = | 站 全 "| 
站 一 3 3 


第 二 步 ,将 户 妆 的 第 2 列 中 的 二 维 向量 (0, 一 3)7 变 成 与 志 一 (1，0) 平行 
的 向 量 . 取 # 一 一 3 一 (一 3， 3 入 一 ng 0 
9, 因 此 ,所 求 反 射 阵 为 


1 0 0 一 3 3 一 
P: = | 总 -| ea | 上 3 "| 
bb 一 1 a 口中 3 


[| 2 | 一 
2 ha 


we wl ls 


1 2 2 
3 3 3 一 3 3 —3 
0= P,P. : 3 3 | 
2 2 1 [| 3 
3 3 
9， 课 


一 1 4 
试 确定 入 及 六 1! 煌 征 值 的 界 ， 
解 ”根据 Gerschgorin 图 盘 定 理 知 ,4 的 特征 秆 4; 位 于 圆 盘 
1 2 一生 | 和 2， 一 1，2，--…，m 
之 内 ,如 2 二 和 实时 一 1 22. m,n. 


故 4 的 特征 值 范围 为 二 太志 


请 


2 
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10. 证 点 一 人 2 】 :及 证 4 为 各 1 的 特征 值 ,4, 是 4: 的 特征 值 , 一 
人 A 


3 2 
(ay ez gs)” 为 对 应 于 入 ， 辣 ( 的 特征 向 量 , 二 《有 ， 讨 )T 为 对 应 于 和 ,各 的 
特征 向 章 ， 求证 ; 

1) Ait A 为 息 的 特征 值 ; 

C2) x 二 (grt gs gr 0 0)7 为 对 应 于 和， 各 的 特征 向 量 ,y' 二 (0, 0,0， 
局 , 房 ) 为 对 频 于 ,入 的 特征 向 量 ， 

证 阴 (1) 4 的 特征 方程 为 detCA4 一 1) = 60, 即 det( 和 :1 一审 :sy det( 和 zs 
一 让 2,:) 二 0. 由 于 4, 为 起 1 的 特征 第 ,4, 是 4* 的 特征 值 , 夺 det( 直 1 一 XT.) = 
位 ，dett 六 ;一 AT 二 上 1 所以， 入 均 是 二 的 特征 值 ， 

(2 = x 


ER TT x,! 
,A 0 | ， 
ac = { ) 0 i= 0 =4.)0 |= ix 
OO vA:: 

0 0 D 

| 0 自 
A (™ 0 ) 总 0 总 
一 一 一 A 一 i¥, 
及 0O 4 加 0 ‘Eo 入 

-J 站 2: 3 时 上 


因此 ,zx 是 对 应 于 ii , 各 的 特征 向 量 ,” 屁 对 应 于 4; ,和 丰 的 特征 向 量 ， 
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第 9 章 常 微 分 方程 初 值 问题 数 值 解法 


1.1 重点 、 难 点 全 村 


1.1.1 基本 概念 


本 章 主 这 钱 究 一 阶 常 微分 方程 梓 值 问题 
y = fr Cr ED= i vy) ar 人 hyeER) 
Ym ys xn EE La, bY 
的 数值 求解 方法 . 根据 常 微分 方程 理论 , 当 让 关于 y 满足 李 普 希 花 条 件 时 , 初 什 
问题 适 定 . 
所 滑 数 值 解法 ,就 是 寻求 ytx) 在 系列 离散 节点 
TY] a Ep 
上 的 近 息 信和 相 邻 晓 节 点 的 间距 训 , 一 on 一 ; 称 为 步 
长 ， -一般 取 等 步 长 的 节点 ,此 时 记 步 长 为 六 
莉 值 问 是 49- 1) 的 求解 都 采用 “ 步 进 式 ”". 数值 方法 有 单 步 与 多 步 之 分 ,也 
有 显 式 与 隐 式 之 分 . 


《9.。 12 


1.1.2 简单 的 数值 方法 与 基本 概念 
欧 拉 法 Ne = yh = 1 2 
后 退 欧 拉 法 加 :一 加 十 有 Feng yo)y 二， 1 2 
梯形 方法 wm 一 > 十 也 BE ew yr yl = 0, 1 2 
改进 的 欧 拉 法 
人 


内 - 
Vt] 二 .ya TT FL Wn Er Yt )] 
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上 述 方 法 都 是 单 步 方法 ,其 中 后 退 的 欧 拉 法 和 梯形 方法 是 隐 式 方法 ,需要 
渤 代 求解 ,其 选民 公式 及 收 雍 性 条 件 详 见 教材 . 
问题 (9. 3} 的 单 步 法 的 一 般 形式 为 


Vr] 一 Ja 十 pt nn， er， 页 《9 2) 
定 兴 % 1 设 ytz) 是 初 值 问题 (9. 1) 的 准确 解 , 称 
了 = yr gtr vr) YT) h) 《98. 3) 


为 单 步 法 (9. 2) 的 局 部 截断 误差 . 
定义 .2 设 yix) 是 初 值 问题 9. 1) 的 淮 确 解 , 若 存在 最 大 整数 户 使 显 式 
单 步 法 59. 2) 的 局 部 截断 误差 满足 
To = yr her DC) (9.4) 
则 称 方法 49. 27 具有 皮 阶 精度 ,也 称 方 法 (9.2) 鸭 p 阶 塘 法. 若 将 (3 4) 的 展 
式 写 成 


Tn = fico MATa he + OUR ) 
刚 称 giri， ytza7 37h 站 ! 为 局 部 截断 误差 主 项 . 
容易 证 明 , 欧 拉 潜 和 后 退 的 欧 拉 法 都 是 一 阶 方 法 ;梯形 法 和 改进 的 欧 拉 鞭 
是 二 阶 方法 . 


1.1.3 龙 格 - 库 增 方法 


显 式 龙 格 - 库 塔 方法 的 一 般 形式 为 
Nt 7 ys hrs Va {9.50 
其 中 


r 


pt 一 DPK, 


7 一 ] 


天 1 一 天 ro CH. 6) 
KFK, = fx, AA, y, HAS KD), i 二 2, 3 rr 
算式 (9.5),(9.6) 称 为 + 级 显 式 龙 格 - 库 塔 方法 、 
常用 的 2 级 2 阶 龙 格 - 库 塔 方法 有 改进 的 欧 拉 公式 和 中 点 公式 .2 级 2 阶 显 
式 龙 格 - 库 塔 方法 处 可 能 达到 三 阶 精 庶 ， 
经 典 的 四 阶 龙 格 - 库 塔 公式 为 
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yo KK KK) 


6 
Ki = fe, $) 
Ke = f(z 十 访 和 3 十 广 1) ， 《9.7) 


K; = f(z 十 序 h， 十 羡 兢 ，) 


Ks = fers 十 睛 ,加 十 RE 


1.1.4 单 步 法 的 收 全 性 与 稳定 性 


定 匀 .3 设 yl(x) 是 初 值 问题 (8,1) 的 准确 解 ,y, 是 其 一 数值 方法 5 如 
(9.2)) 的 数值 解 , 若 对 于 固定 的 = 二 十 大, 当 h 一 0 时 有 yy 一 yx): 则 称 
该 数值 方法 是 收 放 的 . 

定理 9.1 假设 旦 式 单 步 法 

Yl {9. 8) 
县 有 卢 阶 精度 , 且 增 量 函 数 pir， yy 5) 关于 x 满足 李 普 希 兹 菜 件 
i gr yr hy glx yo | Lol yi 

及 设 初 秆 mw 是 准确 的 , 即 mm = yz , 则 其 整体 截断 误差 

Vr) yn = Oh?) 

定义 9.4 车 单 步 法 (9.8) 的 增 量 函数 gg 满足 

Perr pr O00 = Flr y) 
则 称 该 单 步 法 与 初 值 问题 (9. 1) 相 容 ， 

单 步 法 i9.2) 收 襄 的 充 要 条 件 是 此 方法 与 初 值 间 题 (9. 1) 相 容 . 

定义 9.5 单 步 法 (9.2) 用 于 解 模型 方法 y 一 ay .Rel2) 之 0, 车 得 到 的 解 
3 一 ECR)yr 满足 | EL 所 7 ,之 二 , 则 称 方法 (9, 2) 是 绝对 稳定 的 .在 yx 二 品 平 
面 上 ;使 | EL} | 之 1 的 变量 转 成 的 范围 , 称 为 绝对 稳定 域 , 它 与 实 轴 的 公称 为 
绝对 稳定 区 间 ， 

t.1.5 线性 多 步 方法 

一 般 形 式 

<s3 


下 
or = Dearyor +h DR (9. 9 


i 
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四 阶 阿 当 姆 斯 显 式 公式 与 隆 式 公式 


十 下 (55 一 3 二 9 


Ms 二 rz 9 二 19fus — 3f Tf 


有 其 米 尔 尼 方法 . 汉 明 方法 .预测 -校正 方法 等 详 抑 教材 ， 
线性 多 步 方 法 的 构造 方法 主 于 是 基于 数值 积分 和 泰勒 展开 这 两 种 途径 ， 
有 美 局 部 截断 误差 等 单 步 法 的 基本 概念 对 于 才 步 法 也 类 似 成 立 . 


1.2 习题 全 解 


1， 用 黑 拉 法 解 初 值 问题 
y= 二 lv wo 一 
取 步 长 太一 1. 计 算 到 工 一 0.3( 保 留 到 小 数 点 后 4 位 ). 
解 ” 欧 接 法 公式 为 
Br yn hr yh IO00y), n=0,1,.2 
代 一 0 人 下 式 , 计 算 结 果 为 
FD = 0D yO 2 ye = OO Ov ytd. BW a ya = 0.005 0] 
2. 用 政 进 的 欧 术 法 和 梯形 法 解 初 值 问题 
yy 十 Ty vd) =0 
取 步 甘 开 二 00,1: 计算 到 一 0.5, 并 与 准确 解 y 二 一 e 7 十 下 一 + 十 ] 相 比 较 ， 
解 改进 的 欧 拉 法 为 


yo + 记 hLfr,. yO flrs yt hf rs wy 
将 Fz 让 二 十 x 一 y 代 大 上 式 , 得 


hi 
yo 一 (一 二 第 )y1 a ha 二 xLsx 


梯形 法 公式 为 


,A 
Sm Zh" Dh 
将 9 二 81 二 0.1 代 信 上 二 式 , 计 算 结 果 见 下 表 ; 


Ert] es) rmi(tl 二 x1)] 
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一 一 一 
x+， |j 朴 进 欧 拉 污 9 i en} — yn | 梯形 法 > | tr) 一 2 


,1 NA. DO S08 局 337 418 036 > 1073 | oO. O05 238 O93 1 0, 755 132 781 x 10°+ 


.86 G48 TB YY L107 


D2 TO. 927 S00 [0D. B58 2353 O78 x IQ0 2 | 0.021 AN 896 
.3 [0.0 dd 388 160,962 BO 182 x 10 © | 0O.049 367 239 700.185 459 $33 10 5 


Da OM RID 71 | 0. 125 O71 B72 x 10 2 | 0. DR ON3 BG2 | 0. 223 738 443 > 1073 


D0 | 0.144 992 237 00.132 291 S68 > 10 °° [0, 143 722 388 0, 253 Nd8 ET Xx Ly 3 


可见 梯形 方法 比 改进 的 殴 拉 法 精确 . 
3, 用 梯形 方法 解 初 值 问 题 
y Fy=0 
iytB)Y = 1 


a 


证 明和 其 近似 郝 为 y， 二 (工人) ,并 证 明 当 及 一 时, 它 收 级 于 原初 使 问 题 的 准 


确 解 Y= 二 
证 明 ”梯形 公式 为 


0 一 十 去 [Far 3 十 天 rn yat] 


将 Fe = 一 和 代 人 得 


解 得 


因为 y = 二 所 以 


2 一 上 请， 
Mi 人 汇 ) 
对 时 0 以 上 为 步 长 经 # 步 运算 可 求 得 ytzx) 的 近似 值 3%; 需 x 二 项， 
n= 艺 ， 代 人 上 式 有 
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; 2h ， 2 工 
0 
加 [05 


4. 利用 欧 拉 方法 计算 积分 | 加 电 在 点 二 一 吕 5.1.1.5, 之 的 近似 值 ， 


解 。 令 ytx) 一 | 中 , 则 有 初 值 问 题 


y= Se WN} 二 9 
应 用 网 拉 法 , 联 廊 = 0.5, 计 算 公 式 为 
ya = pier n= 0,1,2,.3 
由 y(07 二 二 Di 得 
yO a = 0.5, yl.0) a yo = 1.142 
Yl) a yy = 50, yl2.0) a yr =: 7. 245 
5. 取 和 下 二 0.2, 用 四 阶 绿 典 龙 格 - 库 塔 序 法 求解 下 列 初 值 问题 : 


VI yo 站 之 
1) : (2) ， 


1 
ven 一 了 Loy = 1 


解 。 利用 四 阶 绎 上 典 尝 格 - 库 塔 方法 . 对 于 问题 (1) ,f(r 二 工 十 3 对 于 
间 题 (2), f(z, wy) 二 Hy 取 0.24% 二 0) 三 1. 分 别 计算 两 问题 的 近似 


“了 十 工 
解 见 下 崇 : 
La 《ly 的 租 ya i) 的 解 疝 
二 
.2 1.242 SDD N00 1. 727 548 208 
0Q. 4 1. 583 833 920 2. ?4d2 951 299 
.8 2 44 212 913 4. 094 181 355 
0D.B 2.651 041 8532 5. B29 210 728 
ln 3. 436 302 273 了 . 996 O12 143 


6, 证 明 对 任意 短 数 1, 下 列 龙 格 - 库 诺 方法 是 二 阶 的 ， 
yo = 十 本 AKC + K;) 


KI = fr yy,} 
民 ; 二 Fr 十 块 ，Va 十 的 下 1 
RK, = Fr, tt tl hk, wi] OAK.,) 
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证 明 ”根据 截断 误差 定义 ,只 要 证 明 Ti 一 Ci 即 可 .而 
了 = yr yr hgtr, 3， 有) 


gplr, y: 内) 一 二 [f(r 十 内 yy 十 太 y CT 十 rf 十 (1 一 中 hy 十 


Cl shy 《zy 
困 此 只 须 将 ?Cr 二 ji) 和 pz，yy ph) 都 在 处 展开 即 可 得 到 余 项 表达 式 ， 


zz 十 读 ，y 工 二 (2 一 fly, 二 而 弛 cv， 人 十 
靖 呈 信人 ar vy) Oh 
Fr (ih yt hy (tr)) = 
Fa 二 一 0 人 zz， DTU Dh SE, OCs) 


所 以 
了 rr (二 (CD 十 二 乱 六 (0) 十 于 如 一) 一 


.1 La po | ] 3 
| A y+ 下 二 二 (7 + hy br) Em tr yD) |= OR) 


故 对 任意 参数 1. 此 方法 是 二 阶 的 、 
?7. 证 明 中 点 公式 
op = hk 
琴 ， 一 Fr Yn) 
KK; = fx 十 羡 ， Ys 十 SK ) 


是 二 阶 的 ,并 求 其 如 对 稳定 区 间 . 
解 ”根据 裁 断 误差 定 义 


Tr = Mr ~ yr, ) hf (xz, t 1 YT 上 


Cr 一 


ym) tt hy (x, ) 二 全 (ee 十 亲生 Cr TD) 一 


hh af (rss yr) ， 
ye — hk ferss YC)) 上 这 EE 
Gyr,) Bf lrm pir) | 


IFFY Fr, yr)) 
下 | (名) + 
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hi, fr TY 

BY [| aray 
(yr)) Se | Oths) | 一 
下 wy Ah: 2 ， 他 上 本 ， ys Ee 
EI 【了 医 yr) ay 1 {Cy ry 这 | yn t 


Oy = CRY 
因此 ,中 点 公式 是 二 阶 的 . 
车 yy = AytReta) < 天 0 , 则 
ww 一 [1 十 贡 十 去 Oh |> 


即 当 | 1 十 矣 十 本 (有 < 1, 一 2 过 衣 世 9 时 ,中 点 公式 绝对 稳定 . 
8， 对 于 初 值 问题 
y lO00ty— ri)+ or y(0)=1 
{1) 用 默 拉 法 求解 .上 汕 长 五 取 计 么 范围 的 值 ,才能 使 计算 稳定 ?” 
{2) 著 用 中 阶 龙 格 - 库 塔 方法 计算 ,还 长 月 如 何 选取 7? 
(3) 车 用 实 形 公 云 计算 ,还 长 及 有 无 限制 ? 
解 (1) 欧 拉 法 的 绝对 稳定 区 间 为 
.1 十 (一 1008) 所 1] 
即 Dh 0, 02 
(2) 四 了 脐 龙 格 - 库 塔 法 的 绝对 稳定 区 间 为 


+ 1 2 1 3 了 1 | 
1 证 逢 十 到 ( 守 ) 十 有 上 | 


二 2. 了 785 


A 
{3 梯形 法 的 络 对 稳定 区 间 为 疫 一 一 1007 EE (一 o0.0), 即 


0 


即 一 2.785 委 市 日 < 直 : = 0, 027 835. 


A 


故 对 步 长 无 限制 
9, 分别 用 二 阶 亚 趟 阿 当 姆 斯 方法 和 二 阶 隐 起 阿 当 姆 其 词法 解 下 列 初 值 问 
题 : 
y= 1—y, v0=0 
取 玉 二 0.2, ys 二 0 wn 二 0.181, 计 划 ytl,0) 并 与 准确 解 y 一 1 一 e。 ” 相 比 较 . 


解 ”二 阶 阿 当 姆 斯 绍 式 和 隐 式 方法 分 别 为 
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Jr 这 el 十 二 (37 一 访 )》 
三 
Yr 二 ft 二 让】 


将 了 = 1 一 > 代入 上 二 式 , 化 简 得 


是 式 方法 一 
_ 2 一 2 
愧 式 方法 Jrl 一 DCA" 十 了 二 


取 卡 二 0.2, 二 0 0.181. 计 算 结果 见 下 表 : 


LT 显 式 y， | fr) CO— yn ! 凰 式 入 Cr 一 | 
-1 上 

Dd N326 7 297 995 3 1 3 个. 了 和 六 D. 229 135 x 1073 
DF 站 .和 ?79 站 139 B36 S$ x 10 3 站 ,451 了 :43 Bol 和 555 二 证 7 x 了 和 一 
[| ,545 423 DD. S24 BO 3 x IO-? Lr A 0. 753 F100 x 10™ 
了 .已 | oO.B2f 7351 总 .上 4 5143 8 10°: 0 a2 HB3 2 D. BH G61 x 107? 


可 见 , 隐 式 方法 比 明 式 方法 精确 . 
10. 证 明 解 y 二 六 zw) 的 直列 差 分 公式 


ye Cy 1 


1 
Mot] = ty 


是 二 阶 的 ,并 求 出 局 部 截断 误差 的 主 项 
证 明 ”由 局 部 截断 误差 的 定义 知 


Ta = yix, 十 在》 一 Cw)) 十 CE 一 由 7 一 


Cy Cr 二 外 yr ay nA) 


可 
i 


yr Ti hy Cr 于 和 Cr) 十 


sw 1 ve， 
En ¥Y Cre} TOth 5 YT) 


本 [ye) 一 人 Cr) 十 二 由 ta 一 二 让 om) 十 De ;|- 
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站 | fy rn) thy Cr) + yO )— yr) + 
4 Li 2 


3 fa 一 Ar (Cao 十 于 是 37Ca) 二 + Oh) )]= 


全 一 字 一 二 )ytr+( 计 一 二 1+ 宇 Ym) 十 


1 1 3， _ 
( 言 二 让 一 去 一 时 Jiasrteo) 十 DO) = 
— 于 由 Gan) + Oh!) 


故 方法 是 二 阶 的 ,局 部 规 断 误差 的 主 项 为 一 hz 
11. 试 证 明 线性 二 步 法 
yes CB Dy by, = DL + 3) fo (C3841)7.] 


当 所 天 一 ] 时 ,方法 为 二 阶 , 当 上 56 二 一 1 时 ,方法 为 三 阶 . 
解 ”由 局 部 截断 误差 定 习 知 
Tnz = yrs 2A) TD yr A) yr ) CO— 


Lt OY Cr 2h) + 30+ Dy (Cr)] = 

7) 二 2hy 《zn 二 地 (2 Gr) 十 吉 (2h) ?yr (xz) 十 
TA PCr) OY TF BD EYE) + hy Ces) + 
Epy ryt yt yr Oy | 
2 J 
ey) Bt) | ye) + hy Cr) + y+ 
上 3 1 + 让 和 ” 

译 C2AY YY Cr 十 DG ;| 0+ Dy tr,) 一 


(1 ol—A) ytr)+ 


|2+6 一 1 (二 3) (36 : D |ay’ tx + 


| 
L 
1 1 1] sw 
2 二 到 地 一 1 一 可 示 十 3) | 证 3 十 
I 
L 


1 1 
可 二 百人 季 一 1 一 可 他 十 3) | yr Cr,) 二 


第 9 章 ” 常 微 分 方程 初 值 问题 数值 解法 135 


| 3 十 去 G 一 1 一 也 (6 十 3) we (xn) + Oa!) = 


一 地 人 二 DY) 一 ( 言 一 玖 Phy Cr) 二 OCR) 


当 上 了 关 一 1 时 ， 
1 


村 (十 Day tr,) 十 (CR 


Ts 过 一 


该 方法 为 二 阶 : 当 户 一 一 1 时 ， 
? 


3 ， : 
Te 一 一 《 计 一 去 人 hy te) + O05) 


该 方法 为 三 阶 ， 
12. 未 方程 组 
”一 一 10 dv 
| = 10u— 1iv 
前 责 性 比 , 用 四 阶 尽 -发 广 活 求 解 时 ,最 大 步 长 能 取 和 多少? 


一 10 9 
解 根据 刚性 比 的 定义 ,车 方 各 组 的 扼 降 4 = | 1。 ,1 | 的 特征 值 
满足 条 忻 Ret4y) < 之 007 一 1,， 23, 则 刚性 比 为 
max | Reta,) | 
与 一 Ls 
jn 上 Recau ny | 
和 的 两 个 特征 值 为 
A 一 一 1， Ja 一 一 如 
所 以 与 一 20 
利用 四 阶 恨 - 玫 方 法 求解 时 , 御 对 稳定 区 和 间 在 站 皇 [ 一 2.78，07 ,南昌 二 环 魏 
一 2 98 


一 人 132， 


一 20 
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附 录 


由 测试 题 ( 一 ) 
1, 填空 题 (20 分 ) 
(1) 设 x" 一 2.40315 是 真 值 了 一 2.401 8 的 近似 值 , 则 了 ”有 位 
有 效 数 字 . 
(21 设 所 Xx) 二 工 十 7 一 1, 则 差 商 ( 均 差 YY[0.1,2,33 二 .f/f[0,1， 
213:4] 一 
(37y YT” 的 相对 误差 约 县 x” 的 相对 误差 的 异 . 


td4) 所 会 三 点 Co 了 TD)， (Cryrormy， (rssrxr)) 的 水 平 直 斌 是 


“一 


(5 用 单 节点 的 高 斯 - 勒 让 德 求 积 公式 计算 积分 | 科 二 dx 的 值 为 


n 


2.《410 分 } 用 二 次 拉 格 朗 日 插值 多项式 tx) 计算 sin0. 34 的 值 ,插值 节点 
和 相应 的 函数 秆 是 (0,07,{0. 30,0, 295 5),(0.40,0. 389 4), 

3. tl0 分 ) 求 隙 数 fx) = e' 在 区 间 L0,]1] 上 的 二 次 最 佳 平声 遂 近 密 项 起 及 
相应 的 平方 逼近 误差 :61 :( 小 数 点 后 至 少 保 留 五 位 }. 

4. 《£10 分) 用 二 分 法 求 方程 六 zx = 二 茹 一 x 一 1 二 0 在 [1.0,1,5] 区 间 内 的 
一 个 要 ,误差 腿 e = 107. 

5.(12 分 ) 选取 常数 2,6, 使 得 max | er 一 ar 一 占 达到 最 小 - 

6. (12 分 ) 用 矩阵 三 角 分 解法 解 方程 组 


1 1 2 | 3 
DO 2 ] 2| | 1 
1 —l1 2 211r; 3 
2 2 5 9 了 
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7. 012 分 } 给 定 一 阶 方程 组 初 值 问题 
[> = 2 十 3 


x = ly 
| 一 zt0) 一 1 
{1} 写 出 用 欧 拉 预 估 - 校 正法 求解 此 初 值 问题 的 计算 公式 : 
《2) 取 步 长 中 一 0.1, 计 算 y(0.1) 和 xt0. 1 的 近似 值 . 


8. (14 分 ?证 明 区 间 [a, 门 上 带 权 上 zz) 的 正 交 多 项 式 PCrCa= 1， 


个 根 者 是 单 根 . 且 位 于 区 间 (a,8) 内 . 

自 测 试题 {一 】 参考 答案 

1. 填 上 是 

60103 《4210 《3) 六 y= Df) C50 


2. 0. 493 36 
3, 1,0]2 99 ~ 0, 85] 12x .0.839 18ri; |||: = 0.008 S89 1 


4, N= 6 = 1.25x: 一 1.375 = 1.312 5.7, = 1.343 75, 


TI: 一 1]1.328125ra = 1.320 312 5 
Ee te— 1}lnte—1) 


5, UU 二 Ee—1,p= 


2 
站 ， 工 ] Ts Os Ei 1 
?01 了 略 2) yt011) a ]. 245 ,ct0,1) xr 1. 335 
8. 上 略 
自 测 试题 (二 ) 


1. 填空 是 (20 分 ) 
(1) 襄 rs 一 0 231 是 真 秆 rr 二 00.229 的 近 想 年 , 则 zz, 有 


数字 . 
(2) 求 积 公式 | f(ydz ~ 和/( 当 ) TF 的 乒 数 精确 诬 为 ， 


1.2, 


(3) 设 方程 站) = 训 十 x 一 1, 则 差 商 ( 均 差 }f[0.1,2.3] = 


3+4] 一 


3 


位 有 效 


‘fi 0, 
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(入 4 [1 出 xD 二 
4 2 ”一 一 


Ee DIT] 
(5) 车 Stx) 二 41 , 
诗人 1 志和 3 


是 三 次 样 条 还 数 , 则 = 二 号 一 vc 一 
2.《12 分 ) 已 知 单调 连续 画 数 y 一 f(z) 的 如 下 数据 : 


车 用 插值 法 计算 ,x 约 为 多 少时 Fr = 1{ 小 数 点 后 保留 5 位 ). 
3. (12 分 ) 用 短 阵 三 角 分 解法 解 线性 方程 组 : 


1 2 3 4 rr 
1 3 4 5| |r: 3 
2 1 4diz| 

3 2 5 Lr 


4. (10 分 ) 选取 第 数 a, 使 得 | | e 一 a ! dx 达到 极 小 ,并 求 极 小 值 . 
5.010 分 ) 建立 Gauss 型 求 积 公式 |， dr~ Ay fir) 下 Ar 


86.512 他) 设 & 六 0,z > 了 .证明 


Tri + Bn) 


是 二 01],2 
rta 


是 求 Va 的 三 阶 选 代 方法 

7. (10 分 ; 用 复 化 幸 浦 生 公 式 求 积分 | zlgrdz 的 近似 值 ,要 求 取 7 个 等 
节点 (包括 端点 1 和 ?) ,小 数 点 所 保留 四 位 以 上 . 

8. 《14 分} 用 插值 法 求 在 + 二 0 与 cosx 相 灯 ,在 x 二 也 与 cosx 相交 的 二 次 
多 项 式 Petr) .并 推导 插值 余 项 的 估计 式 


| Pi: (zz) Co— cosr | 把 证 (zx 也) 


身 测试 题 (二 ) 参考 答案 


1. 填空 题 

《12 (232 ALO (C4)3 (SYa 一 太一 3 一 1 
2. x a 1.321 479 

3. xz 一 [1,1, 一 1 了 


4. 4a 二 二 ,最 小 值 为 0.420 838 


En 


TL = D045 336 36] ,7 = 0. B42 193 058 
一 |， 


有 035 301 293,A: = 0, 964 698 706 
6， 略 
7, 80.116 2 


8, P: 《zy 一 一 十 1 
克 


自 测 试题 (三 ) 
1. 填空 题 (20 分 ) 
(1) 要 使 得 V20 的 相对 误 兰 不 超过 0.1%% ,应 取 位 有 效 数字 . 
(29 设 Im = | VTTEdr. 则 用 2 点 高 斯 公式 所 得 近似 值 为 (第 


果 保 留 ? 位 有 效 数字 ). 
(3) 求 方 程 z= fx) 根 的 牛顿 适 代 公式 为 
100 


14) 设 年 汪 [ 
99 


997 
9 上 计算 辣 的 条 件数 tondt 一 


人 5) 已 知 序数 f(x) 一 | se dt, 取 所 一 6.5, 利 用 欧 拉 法 计算 近似 入 fC0. 5) 


2.(10 分 } 已 知 关 x) 的 函数 值 0) = 一 5 7 = 一 3,f( 一 1) 二 一 15, 六 2 
一 一 9, 求 Newton 艾 荆 桂 值 客 项 式 Nytr) 及 站 1,5) 的 近 灿 得. 

3. 515 分} 求 函数 ft.r} = arctanr 在 0;1] 上 的 一 次 最 佳 平方 逼近 钉 项 式 
{小数点 后 至 少 保 留 5 位 }。 
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4. (15 分 ) 给 定数 据 
+ | i130 1. 32 1 34 ] .36 1. 38 


Tr) 3. 02 10 3.903 30 #4, 255 60 673 中 53.177 44 


用 复 化 Simpson 公式 计算 了 = | (zydz 的 近似 值 ,并 估计 误差 . 


de 0 fn 1 | 
5. (10 分 ) 给 定 线 件 方程 组 . p al 四 一 | 其 中 Bed sds sls 
了 
均 为 已 知 常数 .是 cbr 天 0.011 写 出 Gauss-Seidel 选 失 格式 :27 分析 该 先 代 格 式 
的 收 侣 性 ， 


BG aq ee Ti 


r2 2 1 
6. C10 分 ) 求 三 阶 答 阵 |0 2 :| 的 QR 分 解 ,其 中 如 是 正 交 答 阵 ,中 是 上 一 
L2 1 2 


角 矩 阵 ， 
7. (010 分) 证 朋 方 程 fr 一 7 一 27--5 二 0 在 区 阿 (2,3) 内 有 唯一 根 户 


用 区 间 分 半 法 计算 p 的 近似 值 .x, 时 , 试 确定 迁 代 次 数 使 | zu 一 所 | 二 二 义 
10 (不 要 求 计算 mw) 


i = Fioriy) 


8 (10 分 ) 对 常 微分 方程 初 值 问题 | ，， ,所 了 过 刀 使 用 预测 校正 
Val = 
轩 人 — vhf y,) 
A 
公式 3 十 七 [5Fza ,区 了 十 辣 站 一 ri 洒 " 求 其 局 部 截 


断 误差 ,并 指出 该 公式 是 一 个 几 阶 的 公 
自 测 试题 (三) 参考 答案 


1, 填空 题 
. . : te fir,) js _ _ 
Cd F201.642 O42 3 一 一 TPO 9 206 HID.3 


2. 一 4 十 56x 一 3, 一 5 
3. y= 0.0d2 B09 二 O79 Balx 


内 录 id41 
4, 0.343 984 和- —0.27X10° 
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